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Foreword 


This is a time in which the deployment of HV and EHV insulated cable systems is 
exhibiting significant growth and, most importantly, when serious consideration is 
being given to the adoption of such systems in a much greater variety of cases than 
in the past. 

What this implies is that insulated cables are currently looked upon by users as 
an enabling network solution that can help them address effectively, rapidly and 
reliably a variety of complex problems, for example through the deployment of 
mixed overhead and underground lines, rather than, as was often the case in the 
past, a niche solution to be employed only when absolutely inevitable. 

This trend, driven by environmental concerns which have made overhead lines 
less attractive, and facilitated by technological progress which has made insulated 
cables more competitive, requires an ongoing objective evaluation of the possibili- 
ties offered by modern cable and mixed systems. 

Against this background this book represents a timely and useful contribution 
towards a better understanding of the factors that need to be taken into consideration 
when introducing insulated cable systems in the network. 

The book’s thorough coverage of many of the associated technical aspects, cou- 
pled with its analysis of the broader economic, environmental and social factors 
involved, will certainly be of great interest for cable system users, cable engineers 
and students of electrical engineering. 


Milan, May 2009 Aldo Bolza 


Preface 


We shall not cease from exploration 
And the end of all our exploring 

Will be to arrive where we started 
And know the place for the first time. 


T.S. Eliot, Four Quartets 


In the development of the AC power systems, long transmission lines and overhead 
lines (OHL) have been considered for many years an indissoluble binomial, whereas 
the use of HV and EHV insulated cables was devoted mostly to DC submarine links; 
in the last decades, a strengthened sensibility towards the environment, with conse- 
quent big hindrances to OHL installations, and an increased reliability of high qual- 
ity extruded insulations (cross-linked polyethylene — XLPE) of cables have induced 
the transmission system planners and grid owners to install numerous AC HV and 
EHV cable lines: they have unquestionable qualities of adaptability to the problems 
given by the territory when realizing modifications, widening or reinforcement of 
the transmission network. Nowadays, in the jargon of electrical engineering, words 
such as undergrounding, UGC (Underground Cables) and mixed lines (a cascade 
composition of OHL and UGC) have a clear meaning. 

The planning of any new line (OHL, UGC or mixed line) in the electrical grid 
must be always validated by extensive network simulations, among which power 
flow, short circuit, transient analyses; generally those studies ought to be preceded 
by a propedeutical investigation, in order to assist the system planner in the configu- 
ration, comparison and choice of the more suitable solutions, chiefly when the field 
of choice can be very wide as in the insulated cable market. In fact, a target of this 
book was the development of this propedeutical phase regarding the electric energy 
transmission standpoint with particular attention to the operating requirements. 

Starting from the classical equations of transmission lines, developed more than 
a century ago, and by rearranging them with original procedures, the authors have 
created some novel capability charts (for UGCs in Chapter 3, for mixed lines in 
Chapter 4) which immediately visualize the best possible performances within the 
limits of current and voltage compatible with the line lifetime and power quality. 
These capability charts (which take into account the steady state ampacity and volt- 
age levels, subtransient voltages, shunt compensation degree and other important 
parameters) are of primary importance (in the same way of power circle diagrams): 
they offer an original and efficient guide extremely useful for both the planning of 
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UGC or mixed line systems and a full compatibility of their insertion in the power 
network. 

Through the book, there are also several occasions of cultural widening in the 
field of the circuit theory: one of the most meaningful contributions is the recovery 
of Ossanna’s theory, completely developed in 1926 (and in the book newly presented 
in Section 3.9.2). It gives a simple but extremely elegant and powerful analytical di- 
rect computation of the possible voltage regimes obtainable in any AC power system 
by fixing the complex power at one of its ports, so avoiding the iterative procedures 
up till now largely used. This method, particularly useful to complete the line perfor- 
mance analysis, ought to be also generally appreciated as one of the more important 
tool in the circuit theory. 

Also, the authors’ attempt in Chapter 6 of framing the economical comparison 
between UGC and OHL on an engineering fair basis takes aim at helping the in- 
volved stakeholders to clear the field of subjective impressions. 

Briefly, the book offers: 


¢ a wide panorama of the large EHV cable installations which has been enhanced 
by means of some detailed reports of paradigmatic installations (Chapter 1); 

¢ an outline of positive sequence modelling to analyse the steady state regimes of 
typical OHLs, of cross-bonded (with phase transpositions) cable systems and of 
gas insulated lines (Chapter 2); 

¢ call to the constant use of some international rules and national standards; 

¢ novel capability charts (with a lot of examples both for UGCs and mixed lines) 
which offer a throughout visualization of transmission line performances com- 
patible with chosen current and voltage constraints; 

¢ criteria regarding the energization and de-energization phenomena of no-load 
transmission lines (UGC and mixed lines); 

¢ criteria for the dimensioning and check of the shunt reactive compensation of 
UGC (uniformly or lumped); 

¢* comments of educational purpose on the physical meaning of some formal ana- 
lytical expressions; 

¢ introduction to the multiconductor analysis of undergrounding with self-made 
power frequency matrix procedures (Chapter 5), which allows investigating the 
electrical behaviours of all cable conductors (phases and sheaths) with some 
case studies and offers the possibility of comparing and validating simplified 
approaches; 

¢ a technical and economical comparative procedure between overhead and ca- 
ble lines which can contribute to the Environmental Impact Assessment recom- 
mended by the European community (Chapter 6). 


The authors hope to have given the electric engineers a modern and sound tool to 
face the challenges of the future electrical grids where the undergrounding will play 
a key role. 


Padova, June 2009 Roberto Benato and Antonio Paolucci 
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Chapter 1 


HV Cable World Statistics 
and some Large Installations 


1.1 Introduction 


The National Academy of Engineering has rightly stated that the electrical grid is 
the greatest engineering achievement of the twentieth century [1]. 

The invisible presence of electricity and its grid in the industrialized societies is 
absolutely taken for granted in spite of its huge extension and possible vulnerability 
(blackouts!). When thinking of the invisibility of electricity, engineers turn their 
thoughts to cable undergrounding. 

From the onset of the transmission and distribution of electrical power under- 
grounding has always played a key role. 

It is suggestive to remember that the first electrical line supplying Pearl Street in 
New York (1882) constituted of DC underground insulated conductors laid in ducts 
under the street. Moreover, the first known public supply of electricity in Great 
Britain began in 1881 when a small river was harnessed to provide electrical power 
for lighting the streets of Godalming in the county of Surrey. At that time, there were 
no legal powers to enable would-be suppliers of electricity to bury cables under 
streets or footpaths. In consequence, the Godalming cables were laid in the street 
gutters [2]. 

Worldwide, the need for undergrounding installations became more and more 
pronounced with each passing year owing to an ever-growing concern for safety and 
amenities in densely inhabited urban and suburban areas as well as the preservation 
of aesthetic values in many localities. 

The necessity of enforcing and making sturdier the electrical transmission net- 
work and the insurmountable difficulties in erecting new overhead lines force more 
and more towards the use in the power networks of undergrounding links. 

At such level of maturity of the electrical grid the importance of cable systems is 
destined to strongly increase in the medium and long term. 

Many issues [3] can be addressed to justify this asseveration: 


e The ever-growing energy consumptions in conurbations of industrialized coun- 
tries can be only meet by using underground cables; 


R. Benato, A. Paolucci, EHV AC Undergrounding Electrical Power 1 
© Springer 2010 


2 1 HV Cable World Statistics and some Large Installations 


¢ Great difficulties of erecting new overhead lines (even outside the large urban 
areas) due to a strong environmental awareness of the public; 

e In some countries, the exposure limit of the magnetic fields will force the grid 
owners or operators to substitute parts of existing overhead lines with under- 
ground cables; 

e Parts of the existing cable grid constituted of oil-filled paper insulation cables 
(installed up to the end of 1950s) must be replaced and re-investment in under- 
grounding will be needed in the next few years; 

e In many developing countries (e.g. several regions of Turkey), there is a consid- 
erable pent-up demand in respect of reliable power supply systems; 

e The Third World countries will surely expand their electrical grid by predomi- 
nantly using cables. 


In many of the above-mentioned situations, the electrical line is not more monotype 
but is a combination of overhead and cable line. This book will carefully deal with 
this new type of line: the mixed line. 

Another issue deserves a very careful consideration. 

The authors, when facing the question of land constraints for electric infrastruc- 
tures, have often begun with this outstanding passage taken from an EPRI book [4] 
that postulates the necessity to develop common corridors in order to satisfy the 
ever-growing exchanges between areas of big industrialized systems: 


“Common Corridors 


For both ecological and economic reasons future expansion of transmission sys- 
tems especially near urban areas will have to consider sharing their rights-of-way 
with other land users. Independent use of land by single utilities can lead to inef- 
ficient use of natural resources. Development of new rights-of-way in future years 
will have to consider the location of suitable highway networks, gas lines, sewage 
systems, nearby utility transmission lines, telephone lines, underground transmis- 
sion systems, and recreational areas. Such development will require a coordinated 
effort by the government and industry. A current study in [5] outlines the approach 
that the electric utilities may follow to encourage full use of their transmission line 
corridors.” 


This EPRI idea regarding the USA, even if with wide territory resources, is 
strongly recommendable to many other countries. 

This topic is so felt and argued in Europe that the International Council on Large 
Electric Systems (Cigré) has established a Working Group B1.08 entitled Cable 
systems in multipurpose or shared structures [6] and also a Joint Working Group 
B3/B1.09 Application of long high capacity Gas Insulated Lines in structures [7] to 
study the feasibility and the characteristics of these technical solutions. 

The rationalization of fundamental services in the land, such as electricity and 
transport both via rail and via highway, represents a higher degree of conscious- 
ness in conciliating environmental issues and technology. Of course it is also a big 
technical challenge because it is not trivial that there is a full compatibility between 
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possible different infrastructures. Moreover, there are a lot of projects of new gal- 
leries between neighbouring countries (separated by mountain chains) and these 
possibilities represent absolutely unique opportunities for integrating (railway or 
highway) transport and energy transmission within the same corridor or within the 
same structure [8-13]. 

This could strengthen the scarce transnational network between States whose 
electrical networks are responsible for giving high hindrances to an effective inter- 
national electric market. 

Sharing different nature structures and power transmission line within the same 
corridor represents an optimization of service rights-of-way: conciliating technolog- 
ical society needs respecting nature and territory seems to be the challenge of the 
twenty-first century. 

Once again, multi-purpose structures can cause cable technologies to further take 
off. 


1.2 Statistics of Cable Installed Lengths 


Data on the lengths of AC underground cable and overhead line currently installed 
are important in order to have an idea of the consistency and the possible future 
developments of undergrounding in the power grid. 

The following statistics are based on a Cigré Brochure [14] which has collected 
data by means of a questionnaire sent to the Study Committee members. 

The statistics have been divided into five voltage ranges chosen in order to group 
together similar design and operational principles. 

The chosen voltage ranges have been: 


50-109 kV; 
¢ 110-219kV; 
¢ 220-314kV; 
e 315-500kV; and 
¢ 501-764kvV. 


The Brochure highlights that in some cases data capture proved difficult, particularly 
for countries with a multitude of small independent utilities (for example the USA 
and Germany). 

In addition, the national systems for maintaining such data have been discon- 
tinued in some countries, since utilities have been released from state control and 
experienced reorganization, merger and acquisition. 

However, the total length of underground cable circuit expressed as a percentage 
of the total circuit length is shown in Figure 1.1. 

The results show that the large majority of circuits are overhead and under- 
grounding has still a great gap. The proportion of circuits that are underground 
ranges from 6.6% (50-109kV) down to 0.5% (315—500kV).! 


' According with English literature, the decimal numbers are denoted with full stops. 
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50-109 kV 110-219 kV 220-314 kV 315-500 kV 501-764 kV 


Figure 1.1 Percentage lengths of AC underground cable circuits on the total line length for each 
of the 5 voltage levels 
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Figure 1.2 Percentage lengths of AC underground cable circuits on the total line length at HV 
(220-314 kV) for different countries 


Figure 1.2 shows the percentage of underground circuit lengths installed in each 
country in the 220 to 314kV voltage level whereas Figure 1.3 shows it in the 315 to 
500 kV. 

In order to simplify data collection, little technical detail was included in the 
questionnaire. However, where possible, respondents were asked to split the lengths 
of underground cable into those using lapped paper technologies and those using 
extruded polymeric insulation. The results are shown in Figure 1.4. 

The decreasing proportion of extruded insulation used at the higher voltages re- 
flects the incremental development of these cables. 
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Figure 1.3 Percentage lengths of AC underground cable circuits on the total line length at EHV 
(315-500 kV) for different countries 
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Figure 1.4 Percentage of the AC underground extruded polymeric-insulated cable 


Lower voltage, low stress cables were developed first and, as the technology im- 
proved, extruded insulation was applied to higher voltage and used at higher stress. 

Since the early 1960s 50 kV extruded cable has been in use, whereas 400 kV and 
500 kV transmission circuits using extruded insulation were not introduced until 
around 2000. 

The data on installed lengths clearly show that utilities have a strong preference 
for overhead lines rather than underground cables. 

For the 50kV to 109 kV range, 93% of the network is overhead. 

In the voltage range 501—764kV there are not cable installations. 
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Table 1.1 Major EHV Underground Projects (at power frequency) [15] 


Country 


Denmark 
(Copenhagen 
Southern cable 
route) 


Denmark 
(Copenhagen 
Northern cable 
route) 


Germany 
(Berlin/BEWAG 
Mitte- 
Friedrichshain) 
Germany 
(Berlin/BEWAG 
Friedrichshain- 
Marzahn) 

Japan 

(Tokyo (3)) 
United Arab 
Emirates 

(Abu Dhabi) 
Spain (Madrid) 
Barajas 


Denmark (Jutland) 


United Kingdom 
(London) 


The Netherlands 
(Rotterdam) 
Austria 
(Wienstrom) 


Italy (Milan) (6) 


Rated 
(phase- 
to- 
phase) 
voltage 


(kV) 


400 


400 


400 


400 


500 


400 


400 
400 
400 
400 
380 


380 


Type and 
number 
of joints 


dd) 


CPFJ 
72 


PMJ 
42 


CPFJ+PMJ 
48 


CPFJ+PMJ 
30 


EMJ 
264 


PMJ 
12 


CPFJ+PMJ 
96 


PMJ 


Type of 
installation 


(2) 


DB 


DB 


D&M 


DB&D 


DB&D 
DB&T&M 


DB&D 


Route Number 


length of 
circuits 
(km) 
22 1 
12 1 
6.3 2 
5.5 2 
39.8 2 
1.3 4 
(5) 
12.8 2 
14.5 2 
20 1 
2.25 1 
52 2 
8.4 2; 


Conductor 
cross-section 


Winter 
transmission 
total capacity 


(mm?) 
(MVA) 


1600 Cu 
975 


1600 Cu 
800 


1600 Cu 
1100 


1600 Cu 
1100 


2500 Cu 
2400 (4) 


800 Cu 
not available 


2500 Cu 
3440 

1200 Al 
1200 


2500 Cu 
1600 


1600 Cu 
1000 


1200 Cu 
1400 


2000 Cu 
2100 


Com- 
mis- 
sion- 

ng 

year 


mo 


1997 


1999 


1998 


2000 


2000 


2000 


2004 


2004 


2005 


2005 


2005 


2006 


(1) CPFJ — Composite Prefabricated joint, PMJ — Premoulded Joint, 


EM J — Extruded Moulded Joint 


(2) T — Tunnel, DB — Directly Buried, D — Ducts, D&M — Ducts and Manhole 
(3) Cable system prequalified following Japanese Specification 
(4) 1200 MVA/circuit with forced cooling in the future, 900 MVA/circuit now 
(5) 15 core kms/4 circuits x 3 phases = 1.3 km 
(6) See mixed line Turbigo—Rho in Section 1.6 
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Table 1.2 Design criteria adopted for 400 kV Projects [15] 


Project Cable — Conductor Electric Metallic Outer sheath Joints 
circuit stresses screen (1) 
IN/OUT 
(kV/mm) 
Berlin A 1600 mm? 11.5/5.4 Cu wires + PEwithFRV 39 CPFJ 
Cu 5 segments Al laminated 
B 1600 mm? 12.5/6.2 Cu wires + PE with FRV 15 PMJ 
Cu 5 segments Al laminated 
Cc 1600 mm? 12.5/6.2 Corrugated PE FRV 24 PMJ 
Cu 6 segments Al 
Copenhagen 1600 mm? 11.5/4.9 Extruded  PEwithSCL 72 CPFJ 
Cu keystone Lead 42 PMJ 
Madrid A 2500 mm? 11.6/6.5 Cu wires+ PEwithFRL 48 CPFJ 
(Barajas) Cu 6 segments Al laminated 
B 2500 mm? 12.5/7.2 Al welded FR PE 48 PMJ 
Cu 6 segments laminated 
London 2500 mm? 11.6/6.5 Cu wires + PEwithFRL 60 CPFJ 
Cu 6 segments Al laminated 
Jutland 1200 mm? 12.6/6 Al wires + PEwithSCL 96 PMJ 
Al stranded Al laminated 
Rotterdam 1600 mm? 11.8/5.9 Extruded lead PE 6 PMJ 
Cu 
Vienna 1200 mm? 12.1/5.7 Cu wires + PE 30 PMJ 
Cu 5 segments Al welded 
Mixed line 2000 mm? 11.8/6.4 Welded Al PE 66 PMJ 


Turbigo—Rho 


Cu 6 segments 


(1) FRV — Flame retardant varnish; FRL — Flame retardant layer; SCL — semi conductive layer 


1.3 Large Installations of EHV Cable Systems 


In order to show the great development of EHV undergrounding in the last few 


decades it seems particularly meaningful to present the Tables 1.1 and 1.2 [15]. 


A lot of EHV cable installations are part of a longer overhead line (e.g. Turbigo— 


Rho, the Barajas project): once again the importance of mixed line arises. 


1.4 Land and Submarine 150 kV AC Cable Link 
Sardinia—Corsica: SAR.CO 


This installation can be considered as submarine HV AC reference installation [16]. 
Even if the islands of Sardinia and Corsica were partly connected through the 
200 kV DC cable installed in the early sixties (called SA.CO.I. — Sardinia—Corsica— 
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Italy) [17], the need to improve power interchange between the two islands has 
required a new connection to satisfy the highly increased energy demands of these 
two tourist sites. The choice for the connection has been a 150kV XLPE AC cable 
with 400 mm? copper conductor (as shown in Figure 1.6). The total length of the 
connection is 31 km and is composed of 10.5 km of land cable on Corsica, 15 km of 
submarine cable and 5.5 km on land in Sardinia (Figure 1.5). The nominal continu- 
ous rating is I50 MVA. 

For the submarine part a three-core armoured cable has been selected (see Fig- 
ures 1.6 and 1.7), while for the land part three single-core cables (Figure 1.7). 

The maximum sea depth of the line is 75 m. The submarine cable (for construc- 
tional details see Figure 1.7) was buried in the sea bed for the major part of the 
route (about 70%) and most of all in the proximity of the seashore in order to pre- 
vent mechanical damage; additional protection has been applied in some particular 
cases. 

The land cable was laid flat at 200 mm spacing and directly buried in trench at 
a depth of 1.4 m (see Figure 1.8). The metallic screens were cross-bonded. 

The fibre-optics cable runs in a duct of the same trench: successively, it goes on 
embedded in the three-core submarine cable. It is worth noting that the design of 
the land and the submarine cable core is similar: this facilitates the transition joints 
between the single-core land cable and the three-core submarine one. 

Before the installation of such a link, the 50 Hz AC Corsica grid was not syn- 
chronous with that of the Sardinia being linked to the SA.CO.I. system (200 kV DC) 
by means of the intermediate conversion station of Lucciana (North-East Corsica 
coast). 


Bonifacio cs 90 kV 


ae ae 


10.5 km 


15 km 


S. Teresa 


Figure 1.5 Sardinia—Corsica 
(SAR.CO) cable link SARDINIA 
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Submarine cable ¢=207 mm Land cable o=84 mm 


Figure 1.6 AC three-core submarine extruded polymeric-insulated cable and single-core land one 
(SAR.CO) (Courtesy of Prysmian) 


DAnfkwW bd 


LEGEND 
1) Copper conductor (400 mm?) 

2) Semi conductive screen 

3) XLPE insulation 

4) Insulation screen (semi conductive water swellable tapes) 
5) Lead sheath 

6) PE Sheath 

7) Fillers 

8) Fibre-optic cable 

9) Reinforcement and bedding tape 
10) Galvanized steel wires armour 
11) Polypropylene yarn (¢=207 mm) 


Figure 1.7 Constructive details of three-core submarine cable (SAR.CO) 
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Figure 1.8 Installation details of land cable (@ = 83.9 mm) laying (SAR.CO) (Courtesy of Terna) 


Table 1.3 Power cable characteristics (SAR.CO) 


Land Submarine 
Type Single-core Three-core 
Cross-section (copper) 400 mm? 400 mm? 
Copper compacted watertight @ 25.1mm 
Semi conductive layer 
Cross-linked polyethylene (XLPE) insulation ¢ 67.8mm 
Water barrier (hygroscopic tapes) 
Al tape longitudinally welded 
PE with graphite coating ¢ 83.9mm ¢ 207 mm 
Cable mass 8.6 [kg/m] 75 [kg/m] 
Maximum voltage for the equipment 170kV 170kV 
Operating phase-to-earth voltage 87kV 87kV 
Operating phase-to-phase voltage 150kV 150kV 
Insulation voltage level (lightning) T50kV 750 kV 
Short-circuit withstand current 20kA/0.5s 20kA/0.5 5 
Route length 5.5 + 10.5km 15km 
Ampacity 580A S80 A 
Transmission capacity at 1SOkV 150 MVA 150 MVA 
Phase conductor DC resistance at 20 °C 0.047 Q/km 0.047 Q/km 
Sheath DC resistance at 20°C 0.15 Q/km 0.35 Q/km 
Shunt capacitance 0.18 uF/km 0.15 uF/km 


Loss factor at maximum operating temperature 0.001 0.001 
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Table 1.3 collects a lot of data about the design of two cable types used in the 
SAR.CO installation. 

The submarine cable route has taken into account the existing 200kV DC 
SA.CO.I cable route: starting from San Teresa di Gallura in Sardinia and arriving 
in Corsica at Cala Sciumara, the 150kV AC cable crosses the SA.CO.I 200kV DC 
cables only one time and in the remaining part of the route is laid at a distance of 
about 400 m. 


1.5 The Madrid “Barajas” Airport Project (Spain) 


This project can be considered as a reference tunnel installation [18, 19]. 

AENA, the Spanish airport authority, had to extend the international Barajas Air- 
port of Madrid. 

This extension included the construction of two new runways. 

Since the existing 400 kV overhead transmission lines operated by REE would 
have obstructed the incoming airport planes, they had to be dismantled and replaced 
by two cable (each with 2500 mm? Copper conductor) circuits installed in a tunnel 
crossing under the new runways. 

Figure 1.9 shows (inside the circle) the double-circuit cable link (12.8 km) in the 
electrical map of a part of the Spanish grid. 


‘os Angel o 7 soil Cabariiss 
CN PI fr i 
os. Tramesa | 
§ t — Corcedifia 
1 Aero — Cristo oa SAS ce Moratzara! pegos 
= oe SSS Colmenar Viejo. ¥oin)oo32 oo ussin Guadalajara. 
sh sbAPAT Aa ety ~ Guadghjar 
ry. 5 ae Habid —8 SS efi “ NZ iy if 
Navalperal ~ en N27 - Bf of Movers 
Ge Pinares vAldemron oa J gas Minerale 
over urns anal 8 SAN SE OPARACUT 
inte Nuevo > —_VILLANUEVA DE weil foo A 5 \ Og ant 
eet drsas tA CAIADA S/n ANGHUELO 
P Felson ; > - , 
0 Fugrtenovill 
JOSE CA 
Atpoo 


© Alabe Foret 


Villarrubia 
RNE  deSantiago 


Tarancol 
Figure 1.9 Electrical map 
with the indication of the 
double-circuit cable (Barajas: ba Puebla de /Porsa Ata 
12.8km) (Courtesy of RED } / | 
Eléctrica de ESPANA) +—100 km — 
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1.5.1 Milestones of the Barajas Project 


2001 Economical agreement between REE and AENA 
June 2002 Presentation of the project for authorization procedure 
2003 Beginning of the cable laying with the tunnel still un-concluded 


(cable laying concluded in January 2004); 
February 2004 In service 


In order to respect the foreseen timetable, the yard could count on up to 200 
skilled worker group simultaneously, with five Pirelli jointer teams and four ABB 
ones present together on site. 


1.5.2 Technical Characteristics of the Link 


In a part of the route the two single-circuit overhead lines, San Sebastian de los 
Reyes to Moorata (35 km long) and San Sebastian de los Reyes to Loeches (25 km 
long), became a double-circuit. In this part, the overhead lines have been substituted 
by a double-circuit cable i.e. by an underground link 12.8 km long. 

The six single-core cables of the double-circuit are installed in the same tunnel 
with the possibility of forced air cooling. The double-circuit is not electrically in 
parallel. 

Each cable circuit contains 48 pre-fabricated joints, sectionalized, straight and 
earthed, and six outdoor terminations plus all required bonding cable and Surge 
Voltage Limiters (SVLs) for the cross bonding arrangement. The ABB joints are 
equipped with capacitive sensors in order to perform partial discharge measure- 
ments. 

Additional current transformers and surge arresters for differential protection of 
the total cable system are also installed. 

Table 1.4 sums up the chief characteristics of the installation. 

The bonding of sheaths consists (see Figure 1.10) of a combination of 5 x3 minor 
cross-bonding sections (810m long) and two end single-point bonding sections of 
300 m and 400 m long. These two sections have allowed to absorb the possible route 
variations (with respect to the project) and to order 90 drums of equal length. 

The cross-bonding is performed without phase transpositions. 

Even if the winter ampacity reported in Table 1.4 is rather high (i.e. 2482 A) the 
link is at present operated at an average current of about 600 A per circuit. 


1.5.3 Tunnel and Earthing System Characteristics 


The tunnel is buried at about 2m depth and in a part of the route runs parallel 
to a river so that some pumps have been installed (251/s) for the possible water 
extraction in emergency conditions. 
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Table 1.4 Barajas installation characteristics 


Maximum voltage for the equipment On 420 kV 
Phase-to-phase operating voltage Un 400 kV 
Phase-to-earth operating voltage Uo 230kV 
Withstand short circuit 50kA/0.5 s 
Voltage insulation levels Waveform 250/2500 ps 1050 kV 
Waveform 1.2/50 1425kV 
Route length 12.8km 
Ampacity (Winter) 2 X 2482 A 
Transmission capacity at 400 kV (Winter) 2 x 1720 MVA 
Joule power losses in the double-circuit (Winter) 448.6 W/m 
Ampacity (Summer) 2 x 2006 A 
Transmission capacity at 400 kV (Summer) 2 x 1390 MVA 
Joule power losses in the double-circuit (Summer) 294.7 W/m 
Inlet air max. temperature (Summer) +42°C 
(Winter) +25°C 
Max temperature of the air in the tunnel +50°C 
Maximum forced cooling air speed 5m/s 


Single point bonding Insulated sheaths 


Major sections 


300m 2430 m 2430 m 2430 m 2430 m 2430 m 400 m 


Figure 1.10 Sketch of the sheath bonding and earthing (Barajas) 


In a location, the route has required a crossing of the river by means of a circular 
section viewable tunnel. The overall digging volume has been equal to 650000 m°. 
The tunnel has been constructed by 6500 pre-cast reinforced concrete units with 
inner dimensions 2 x 2.25 m (length x height) and thickness of 250 mm. 

Some photographs during tunnel construction are shown in Figure 1.11. 

For the earthing grid, the steel reinforcement of the gallery has not been used; 
it has been preferred a formation of four copper conductors (6 = 120mm? and 
¢ = 35 mm°’) located at the tunnel vertexes; a metallic ring bonds the four earthing 
conductors at regular intervals (400 m). 

The tunnel has several technological equipments e.g. closed-loop camera, smoke 
detector system and telephone repeaters (see Figure 1.12) for maintenance and in- 
spections. 
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Figure 1.11 Tunnel under construction (Barajas) (Courtesy of RED Eléctrica de ESPANA) 


1.5.4 Power EHV Cables 


The Barajas power cables have a cross section of 2500 mm? and their character- 
istics are reported in Table 1.5. The two manufacturers ABB and Pirelli have de- 
livered a circuit each. Some photographs of Barajas samples are shown in Fig- 
ure 1.13. 
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Figure 1.12 Detail of the 
telephone system (Barajas) 
(Courtesy of RED Eléctrica 
de ESPANA) 


Table 1.5 Power cable characteristics (Barajas) 


Manufactures ABB PIRELLI 
Cross-section of single-core Copper cable 2500 mm? 2500 mm? 
Milliken 6 sectors conductor ¢ 65mm ¢ 65mm 
Semi conductive layer ¢@ 70mm @ 71.6mm 
Cross-linked polyethylene (XLPE) insulation @ 125.9 mm ¢ 122mm 
Semi conductive layer ¢ 129.7 mm ¢ 126mm 
Water swelling tape - ¢ 128.3mm 
Aluminium sheath - ¢@ 138.8mm 
Copper wires + aluminium foil ¢ 130.9 mm = 

PE outer jacket ¢ 148mm ¢@ 142.5mm 
Cable mass 37 [kg/m] 40 [kg/m] 
Drum length 810m 

Numbers of drums 90 

Drum mass & 35 x 10° kg 


Drum dimensions @ = 4.4m; width = 3m 


Figure 1.13 Two samples of the “Barajas” Pirelli EHV cable (Courtesy of Prysmian) 
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1.5.5 Cable Laying in the Tunnel 


Figure 1.14 shows a suggestive sight of the double-circuit in tunnel whereas Fig- 
ure 1.15 offers a sketch of the same sight with the detailed indications of the spac- 
ings and tunnel dimensions: the cable spacing is 0.5m whilst the circuit spacing is 
1.44m. 

The cables are fixed at the tunnel lateral walls every 6 m with a snaking of 0.25 m: 
cables rest on saddles and brackets. 

Due to the mechanical withstand to the short circuit electrodynamic forces, every 
3m the cable are tied together by means of three phase spacers (see Figure 1.16). 

The four inlet forced cooling stations (a Pirelli project) are located along the 
tunnel at a distance of 2480 m (see Figure 1.17): each station is equipped with three 


Figure 1.14 Double-circuit 
cable in tunnel (Barajas) 
(Courtesy of RED Eléctrica 
de ESPANA) 


2.25m 


Figure 1.15 Sketch of the 
double-circuit cable in tunnel 
with indication of the spacings 
(Barajas) 


2.00 m >| 
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Figure 1.16 Cable snaking 3m 3m 
(Barajas) IK >I >| 


Saddles and brackets ] 
Phase spacers 


Figure 1.17 Forced cooling 
system (Barajas) 


fans (38.3 kW rated) with a maximum air speed of 5 m/s. The outlet cooling stations 
are three. 

The temperature control is driven by the RTTR (Real-Time Thermal Rating) sys- 
tem. 


1.5.6 The Transition Compounds and Protection Schemes 


The two transition compounds cover an area of 60 x 40m? and are not equipped 
with disconnectors or circuit breakers hosting only the oil outdoor terminations and 
the surge arresters (see Figure 1.18). The outdoor terminals are supplied by the two 
manufacturers, whereas the six surge arresters are from ABB. The current trans- 
formers for the differential relays are toroidal-type. 

The transition compound with a unique earthing grid (designed for 50 kA short 
circuit current) hosts a civic building (8 x 2.5 x 2.75 m*) and an emergency diesel 
generator 250 kVA rated (other five emergency diesel generators have been installed 
in the five sites of forced air adduction). 
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Figure 1.18 Transition compound (Barajas) (Courtesy of RED Eléctrica de ESPANA) 
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LEGEND 
87: Differential relay 


21: Distance relay 
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Figure 1.19 Relay architecture (shown only for one circuit) (Barajas) 


In the utility building there are a 25kV/400V transformer, service switch- 
board, telecommunication and cable protection panels, tunnel monitoring terminals 
(closed-loop cameras, fibre-optics terminal, Data Thermal System and Real Time 
Thermal Rating). 

The protection architecture is shown in Figure 1.19. 

With regard to spare units, in one area given by the Airport authorities three 
drums and five joints for each circuit have been stored. In the ventilation station, 
a spare fan is also present. 
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Successively to the conclusion of the project, in the electrical stations of San 
Sebastian and Morata it has been necessary to install two three-phase shunt reactors 
of 150 Mvar in order to perform a good operation of the link. 


1.6 380kV Double-Circuit Cable of Mixed Line Turbigo—Rho 
(Italy) 


This underground line can be considered as a directly buried reference 400 kV in- 
stallation [20, 21]. 

The need of a 380 kV transmission line between the important network-node of 
Turbigo (reinforced by a 800 MW power plant) and the substations of Ospiate and 
Bovisio (which play a fundamental role in the power delivery of Milan City net- 
work) had been felt from several years, aiming at avoiding congestions subsequent 
to the high power flow in that area and reinforcing the West-East transmission axis 
in the Northern Italy. 

The Turbigo—Rho link can be considered as a mixed line composed of a 20 km 
overhead line from Turbigo to Pogliano Milanese and 8.4km double-circuit cable 
line from Pogliano Milanese to Rho: the undergrounding was needed in order to 
cross both strongly populated areas and natural reserves. 

In order to fulfil the thermal power rating of the overhead line, it was necessary 
to install an underground double-circuit cable with a 2000 mm? copper conductor. 

The two circuits were directly buried in two separate trenches (circuit spacing 
equal to 6 m) at a depth of 1.5 m, flat-type with 0.35 m cable spacing: the 6 m circuit 
spacing allows giving a good independence of the two circuits both in the thermal 
regimes and for maintenance purpose. 

Outdoor terminations were installed in the overhead/underground transition com- 
pound. Continuous temperature monitoring with a DTS system and continuous par- 
tial discharge monitoring at joints were also installed. 

A special magnetic field shielding has been provided for some parts of the route 
in accordance with the Italian laws [22, 23]. 

Figure 1.20 shows the electrical map with the indication of the mixed line route. 

The project construction started in July 2005 and has been in service since 30 
June 2006. 


1.6.1 Milestones of the Turbigo—Rho Project 


September 2004 — Authorization to the construction 
July 2005 Opening of the yard 
30 June 2006 In service 


It is worth noting that the short time for the realization (only 12 months) has been 
driven by a very precise schedule of the digging times. 
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Figure 1.20 Electrical map of Turbigo—Rho 380 kV mixed line 


Table 1.6 Turbigo—Rho installation characteristics 


Maximum voltage for the equipment Un 420 kV 
Phase-to-phase operating voltage Un 380kV 
Phase-to-earth operating voltage Up 220 kV 
Withstand short circuit 50kA/0.5s 
Voltage insulation levels Waveform 250/2500 ps 1050 kV 
Waveform 1.2/50 1425kV 
Route length 8.4km 
Ampacity 2x 1600 A 
Transmission capacity at 380 kV 2 x 1053 MVA 
Joule power losses in the double-circuit 245 W/m 
Laying depth (to bottom of the trench) 1.5m 
Laying arrangement Flat 
Phase spacing 0.35m 
Circuit spacing 6m 


1.6.2 The Undergrounding Link of the Turbigo—Rho Mixed Line 


For the underground part, the choice of the cable cross-section has been driven by 
the necessity of avoiding, in any possible operating condition, a bottleneck from 
a power rating standpoint. Since the overhead line can carry at thermal limit (in 
the cold months) a current of 3200 A (in accordance with the Italian Standard CEI 
11-60 [24]), a double-circuit cable with 2000 mm? cross-section has been chosen. 

Table 1.6 reports all the chief installation characteristics at a glance. 

Another important feature of the installation regards the partial discharges (PD) 
monitoring. All installed accessories have an embedded capacitive sensor and are 
permanently connected through a special cable to the monitoring system. Today 
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there are no standards for PD measurement on-site, mainly due to external and back- 
ground noise level that can be neither foreseen nor limited below a given threshold. 
All accessories are connected to several local acquisition units, and all the local 
units are connected to an Ethernet LAN in S-ring configuration. Such configura- 
tion allows continuity in PD monitoring activity even if one or more channels of 
one or more acquisition units should experience any problem. In order to continu- 
ously monitor the performances and the behaviour of such an important connection, 
the system is provided with a variety of monitoring systems, e.g. a DTS system 
for thermal monitoring which collects all the information about the cables physical 
variables by means of suitable sensors (in this case, fibre-optic cables laid on top 
of the hottest phase). However, the most innovative of them is the above-mentioned 
permanent monitoring of partial discharges. 


1.6.3 Power EHV Cables 


The solution chosen for the underground part is a single-core XLPE-insulated cable 
designed for the voltage levels of Table 1.6. 

The cable (see Figure 1.21) has a Milliken type (with six segments) copper con- 
ductor, with a cross section of 2000mm7, longitudinally watertight by means of 
hygroscopic tapes. 


NYDN fF WN 


LEGEND 
1) Cooper Milliken conductor (water-blocked) 
2) Semi conductive screen 
3) XLPE insulation 
4) Semi conductive screen 
5) Semi conductive water swellable tapes 
6) Welded Aluminum Sheath 
7) PE outer sheath 


Figure 1.21 A sketch and a sample of the underground part of the “Turbigo-Rho” EHV mixed 
line (Courtesy of Prysmian) 
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Insulation consists of extruded cross-linked polyethylene (XLPE) suitable for 
operation at conductor temperature equal to 90°C. 

Under the metallic sheathing, a longitudinal water barrier is applied in order to 
limit the water penetration along the power core in case of cable damage. 

The cable is provided with a longitudinally welded aluminium sheath, thick 
enough to withstand the rated phase-to-ground short circuit current (i.e. 50kA per 
0.5s). 

The metallic sheath acts also as a radial water barrier. The anti-corrosion protec- 
tion consists of an extruded PE sheath. The cable outer sheath has an overall thin 
layer of graphite to permit the after laying voltage test of the sheath. 

Table 1.7 reports other characteristics. 

Metallic sheaths of the cables are electrically connected in a cross bonding sys- 
tem. The cross-bonding is performed with phase transpositions (see Figure 1.22). 

Each phase consists of 12 cable spans with four major cross bonding sections 
composed of 710m long minor sections. 

The minor sections are connected using sectionalized joints. At these points the 
cable sheaths are cross-bonded via suitable link boxes, provided with SVLs. 


Table 1.7 Power cable characteristics (Turbigo—Rho mixed line) 


Manufacture PRYSMIAN 
Cross-section of the single core copper cable 2000 mm? 
Milliken 6 segments conductor ¢ 57mm 
Semi conductive layer ¢ 61.2mm 
Cross-linked polyethylene (XLPE) insulation ¢ 113.8mm 
Water swelling tape ¢ 119.3 mm 
Aluminium sheath ¢ 121.3mm 
PE outer jacket ¢ 132.3 mm 
Cable mass 31 [kg/m] 
Drum length 710m 
Numbers of drums 72 

Drum mass ~ 23.7 x 10° kg 
Drum dimensions @ = 4.2m; width = 2.5m 


Figure 1.22 Sketch of the sheath bonding and earthing: cross-bonding with phase transpositions 
(Turbigo—Rho mixed line) 
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Straight joints are used to connect the major sections. At these points the cable 
sheaths are solidly bonded and earthed through three-core disconnecting link boxes. 

At the terminations, the cable sheaths are connected to earth using disconnecting 
single-core link boxes. 

Joints are of premoulded type, with the electrical part consisting of a single piece 
sleeve. As already mentioned, all accessories have been provided with a capacitive 
sensor in order to allow partial discharge monitoring. 


1.7 Cable Laying 


The cables are chiefly running parallel to roads of local authorities property (Prov- 
ince of Milan, Municipalities) with many crossings (roads, rivers, gas pipes, medium 
voltage cables, railways, etc.). 

As already mentioned, the cables are installed in two separate trenches located 
on the opposite sides of the roads, 6 m spaced out. 

With regard to the main road installation configuration, the cables are laid in 
flat formation in a sand backfill at 1.6m depth protected by concrete slabs as in 
Figure 1.23. 

With regard to the road crossings, the “open trench” technique has been adopted: 
the trench is similar to that of Figure 1.23 but the cables are laid inside polyvinyl 
chloride (PVC) pipes blocked in flat formation by concrete. 

For river and road crossings the directional drilling technique has been used: the 
cables are laid inside PE pipes. For the crossing of rivers Olona and Lura the pipes 
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Figure 1.23 Standard cable laying in trench (Turbigo—Rho mixed line) 
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Figure 1.24 (Turbigo—Rho mixed line) Cable laying for the reduction of magnetic field (screening 
raceway) (Courtesy of TERNA) 


Figure 1.26 (Turbigo—Rho mixed line) Pulling winch and cable laying (Courtesy of TERNA) 
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Figure 1.27 (Turbigo—-Rho mixed line) Joint location and trench (concrete slabs are also visible) 
during installation phase (Courtesy of TERNA) 


are filled with water: in order to improve the heat dispersion, these pipes are suitably 
fitted with expansion tanks. 

Another important constraint that the installation has overcome is the limit level 
of 3uT imposed to Electromagnetic Field (EMF) by the Italian law [23] for new 
transmission lines at sensible receivers (i.e. where human presence lasts more than 
four hours per day). A special screening solution has been adopted in order to mini- 
mize the magnetic field in the sensitive areas: the cables are installed in a trefoil con- 
figuration inside single plastic ducts which are embedded in a special high magnetic 
permeability steel raceway suitably protected against corrosion (see Figure 1.24). 
This ferromagnetic raceway is composed of different elements: the hull, the cover 
and the closing clips. The design of the described open raceway combines the closed 
perimeter shielding efficiency with an open shape and the absence of welding. The 
particular shape allows also the raceways to follow the curves of the trench and lat- 
eral and vertical variations of direction. It is worth underlining that the dimensions 
of the raceway are tailored for this project. Indeed they have to reach the optimum 
compromise between the magnetic field screening effect and the minimization of 
the derating factor for the power cables. 

Figures 1.25—1.27 show different phases during installation. 


1.7.1 The Transition Compounds and Protection Schemes 


Figure 1.28 shows the transition compound (which realizes the scheme of Fig- 
ure 1.29) where there are not circuit-breakers but only disconnecting devices. In 
the same Figure 1.28 it is possible to see the surge arresters. 
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Figure 1.28 (Turbigo—Rho mixed line) Transition compound (Courtesy of TERNA) 


Transition Transition 
compound compound 


a lt 
= | —Hibre-ontios _ _ 


TURBIGO STATION 
OSPIATE STATION 


I 

I 

I 

i 

I 

I 

Fibre-optics. | 
. =f | 

== pe Fibre-optics_ "77 La 


Fae See een Sener 


1bre-op! ee es 
Peete tgs 8 I - ‘ I —— J 
ee ee ee Fibre- a sieeeeed -- 
| J ae, te, pees 
Foy ae eee een) |___Fibre-optics Block 
79R} Block Block LZ9R 
79L 7a. 
l Dekada kod peabs ad veces asain einen sae J 
Overhead line Double-circuit cable Overhead line 


LEGEND 
87: Differential relay 


21: Distance relay 
79: Reclosure system 


Figure 1.29 (Turbigo—Rho mixed line) Relay architecture 
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Another requirement has been to assure the highest level of reliability, also con- 
sidering the protection systems currently in use for 380kV. In particular the se- 
lectivity of the protection devices for rapid and slow re-closure cycles had to be 
maintained: in case of permanent faults (on the cable section) the line shall not be 
re-closed, but for faults on the overhead section (generally transient events) the cycle 
has to be completed. 

All these features are present in the relay architecture scheme of Figure 1.29. 

The different relays can detect whether or not the fault is on the cable section. If 
so, they block the (slow and fast) re-closure cycles. 
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Chapter 2 


The Positive Sequence Model 
of Symmetrical Lines 


2.1 Introduction 


In the power frequency steady state regimes of a power line, the positive sequence 
model is often assumed: this model is almost always acceptable for studies of bal- 
anced regimes of overhead lines (even if they are asymmetrical and with earth wires) 
whereas, in the case of single-core, unarmoured UGCs, it can be accepted only if 
cross-bonding with phase transpositions is performed in order to both zero the in- 
duced current circulation in the sheaths and symmetrize the system. 

On the contrary for the Gas Insulated transmission Line the solid-bonding con- 
nection is required: the lowest value of the enclosure resistance allows the induced 
current in each enclosure to be almost equal in magnitude and opposite in angle 
to the corresponding phase current, so that a very good approximation consists in 
considering each phase as decoupled from the others (not only capacitively but also 
inductively), and hence in modelling the line at balanced regimes by means of the 
positive-sequence circuit. 

After these remarks, it appears licit to turn to the classical treatment of the AC 
uniformly distributed parameter transmission lines. 

However, it is worth remembering that also in the aforementioned cases, the 
determination of equivalent parameters to be introduced in the positive sequence 
model hold a given degree of uncertainty: the multiconductor analysis together with 
experimental results can be helpful to deepen, with more accuracy, the knowledge 
of the phenomena. 


2.2 The Transmission Matrix of a Uniform Line 


Let us consider a line at power frequency steady state regime as an uniform sym- 
metrical transmission line: Figure 2.1 shows the single-phase model with distributed 
parameters at positive sequence, where each infinitesimal element dx is qualified by 
the four parameters (also called primary constants) r, £, g, and c. Their values can 
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Figure 2.1 Single-phase model at positive sequence of a line 


be computed after Sections 2.3.1—2.3.4 in case of cross-bonded cables with phase 
transpositions and after Sections 2.4. 1—2.4.4 in case of solid-bonded GILs. For OHL 
the formulae are well-known. 

The voltage and current references at the sending- and receiving-ends (both of 
the entire line and of the infinitesimal element) are also defined in Figure 2.1. 

In the steady state regime each element of infinitesimal length dx has, at w an- 
gular frequency, a longitudinal impedance and a shunt admittance given by: 


zdx = (r+ jol)dx , 


2.1 
ydx = (g +jmc)dx ; ey 


where Zz is the per unit length (or kilometric) longitudinal impedance expressed in 
Q2/km and y the per unit length (or kilometric) shunt admittance in S/km. 

The voltage increase dU,, between the two ports of the infinitesimal element 
dx (located at distance x [km] from the receiving-end R) is the voltage drop on the 
longitudinal impedance zdx so that 


dUy, = 1,zdx . (2.2) 


The current variation d/_, between the two ports of the infinitesimal element is equal 
to the current derived from the infinitesimal shunt admittance 


dl, — ox an dU y,] * y -dx. (2.3) 


By considering that in (2.3) the term of superior order dU 9, - ydx can be neglected, 
(2.2) and (2.3) can be rewritten in the following form: 


dU 

ae ae oe (2.4) 
aly = yU (2.5) 
dx So " ; 


It is worth remembering that the steady state phasors Uy, and J, rotating with 
constant angular frequency w determine in the time domain (with their projections 
on a reference axis depending upon their phase displacement) the sinusoidal in- 
stantaneous values. Consequently it remains only to solve the system of differential 
equations composed of (2.4) and (2.5) as if the current and voltage phasors were 
complex functions of the sole spatial variable x. 


2.2 The Transmission Matrix of a Uniform Line 31 


By £-transforming (C = Laplace) the complex functions Uy, = U(x) and 
I, = I(x) of the x real variable and remembering the well-known transform of 
a function’s derivative i.e. 

d 
F f(x) 
dx 


= s£L6f(o)} — £O) 


the following system arises 


SL{U 9. $ — Upp =2-L{1,} , 
SEAL ds _ y- LLU ox} , 


where (with the abscissa origin at receiving end) the obvious relations U 9(,=9) = 
Gor and I (<9) = Lp have been applied. 
Rearranging and reordering, it has 


SLLU),} — ZL{L,} = Uop, 
—yL{U,} + SLL = Ir ’ 


from which it is immediate to obtain (by means of Cramer’s rule) 


, (2.6) 


SIpt+ WUor 
Srey 


Lihat = (2.7) 


If the new quantities k and Z,, are defined 


so that: 


——— (2.8) 


I ae 
LU} = —— =. (2.9) 
S _— 
By implementing the inverse Laplace transform of them, the current and voltage 


phasors as a function of x are expressed by 


Un, =A, Upp + By Ip: (2.10) 
iS C.6U Doe ks.. (2.11) 
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where 
A, = cosh(kx) ; B, = Z,-sinh(kx) ; (2.12) 
C, = [sinh(kx)]/Z,; Dy = cosh(kx) 
x [km] distance to the receiving end of the line; 
k= /z-y [km] propagation factor; 
L.= 5 [Q] characteristic or natural impedance. 


If the hyperbolic functions are expressed as exponential ones (by using the Euler— 
Lambert relations) it is possible to obtain the transmission formulae (solution of the 
telegraphist’s equations) which were first written in 1876 by O. Heaviside [1] and 
then used, rearranged and divulged by Arnold [2], Réssler [3], Steinmetz [4] et al. 


If x = d, it yields 


Voxaad)=Uos and Log =Ls 


Lys 
so that 
Ups =A:Unrt+ B-Lr, (2.13) 
Is =C-Ugrp+D-Ipr., (2.14) 
where 


= Z,-sinh(kd) : 


A B 2.15 
C = [sinh(kd)]/Z,; D = cosh(kd) . ac 


The single-phase model of the line (of length d) can therefore be considered as 
a two-port network which is characterized by means of A, B, C, D parameters, also 
called as transmission hybrid parameters since they have different units (and differ- 
ent circuital meanings) as highlighted in any power system analysis book [5, 6]. 


Since the line has uniformly distributed parameters, the two-port network results 
symmetric so that: 


A=D. 


Equations 2.13 and 2.14 can be rearranged in the concise elegant matrix form, where 
M is often named transmission matrix, 


Uos |_| AB] | Vor 
Is 7 le 5 IR a) 
—— —_—- 

Ws M WR 


which allows obtaining the hybrid vector wy at sending-end once fixed the hybrid 
vector W p at receiving-end. 
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The relations (2.15) allow verifying that 


det M = AD — BC = [cosh(kd)]* — [sinh(kd)]? = 1, (2.17) 


so confirming that the two-port network constituted by the line obviously satisfies 
the reciprocity principle. 


2.3 Computation of Single-Core Cable Kilometric Parameters 


As already stated and remarked in Sections 2.1 and 2.2, the computation of AC 
primary parameters is referred to a cross-bonded (with phase transpositions) single- 
core unarmoured cable. 


2.3.1 Computation of r (Cable) 


It is important to reaffirm that with the hypothesized perfect cross-bonding the in- 
duced sheath currents are zeroed and that, as is well-known, the alternating current 
Joule losses in the phase conductor are (due to the skin and proximity effects) higher 
than those in direct current. 

This fact brings to consider for the phase conductor a kilometric resistance rac 
greater than that in direct current rgcyo as shown in [7](IEC 60287): 


lac = Tdcd° (1+ ys + yp) [82 /km] ’ (2.18) 
where: 
ys = skin effect factor; 
yp = proximity effect factor; 
Tacge => DC resistance computed at the maximum operating conductor 


temperature 1° 


Tac => kilometric cable resistance at maximum operating temperature 0°. 


Table 2.1 Permissible conductor temperatures for HV and EHV cables [8] 


Cable insulating material Maximum operating Short circuit 
temperature [°C] temperature [°C] 
Impregnated paper 85-90 160-180 
Low-density polyethylene — LDPE 70 150 
High-density polyethylene - HDPE 80 180 
Cross-linked polyethylene — XLPE 90 250 


Ethylene Propylen Rubber — EPR 90 250 
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Table 2.2 Electrical resistivity at 20°C and temperature coefficients of materials used in cable 
construction [7] 


Resistivity Temperature coefficient 

Material used in P22 Q m at 20°C a 
Conductors 

Copper 1.7241 x 1078 3.93 x 1073 

Aluminium 2.8264 x 1078 4.03 x 1073 
Sheaths and armours 

Lead or lead alloy 21.4 x 1078 4.0 x 10-3 

Steel 13.8 x 1078 4.5 x 1073 

Bronze 3.5 x 1078 3.0 x 1073 

Stainless steel 70 x 1078 Negligible 

Aluminium 2.84 x 1078 4.03 x 1073 


Usually it is conservative to assume that 3° is the maximum operating conductor 
temperature [8] compatible with the insulating materials (see Table 2.1). It should 
be noted that XLPE and EPR are the most interesting dielectrics. 

For the computation of rg.yo the following well-known formula can be used start- 
ing from the kilometric resistance rz9° at 20°C 


race = F290 - [1 +a- (0° —20°)] [Q/km] , (2.19) 
where: 


a => temperature coefficient 1/°C; 


v° => maximum operating conductor temperature. 


With regard to the resistivities and temperature coefficients of the metals used in 
cable construction, Table 2.2 sums up the values [7]. 

The computation of y, and yp, must be performed after the detailed formulae of 
IEC 60287 [7]. 

With regard to the large cross section conductors, it is worth remembering a seg- 
mental conductor design (Milliken H. Electrical Cable. Patent Publication 1933, US 
1904162), which is commonly referred to as a Milliken conductor: it has insulated 
segments consisted of layers of wires (more or less insulated each other) stranded 
around the segment “centre” aiming at reducing the magnitude of skin and proxim- 
ity effects [9, 10]. 

Nowadays, cross sections up to 2500 mm? are customarily used and this topic 
continues to be reviewed [11] and developed so that any formula must be used with 
due caution (possibly supported by experimental validations of manufacturer). 


2.3.2 Computation of £ (Cable) 


In the hypothesis that the current distribution is uniform in the whole cross-section 
of the phase conductor and that the sheath currents are zeroed by the cross-bonding, 
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Figure 2.2 Single-core cable @) 
laying arrangements 


Open trefoil arrangement 
flat arrangement with axial spacing s 


the kilometric inductance (to be given to each phase) can be computed by means of 


K-s 
£=0.2In (ca) [mH/km] , (2.20) 


where: 


S => axial spacing (mm); (see Figure 2.2); 
K = 1 for trefoil arrangement; 
K => 1.26 for flat laying (this value can be easily confirmed by means of 
the Geometrical Mean Distance among conductor centres 
GMD = ¥/3)2* 523° S13 = A/s +s +25 = 1.26-5 when the cross- 
bonding is performed with phase transpositions); 
GMR => geometrical mean radius (mm) of the phase conductor. 


Usually in EHV XLPE-cables the conductor is hollow Milliken type; for GMR, 
see [12], Vol. II, page 155. 

It is worth noting that the axial spacing can be reduced to a minimum (i.e. touch- 
ing trefoil) but with a decreased ampacity. 


2.3.3 Computation of c (Cable) 


The values of the single-core cable kilometric capacitance c (between the phase 
conductor and the metallic sheath) can be computed by considering the cylindrical 
symmetry, so that (with ¢9 = 8.8542 - 107!7F/m) 


27 E0Ey 


(8) 


(c) [F/m] > c = —“*—__ [yF/km], (2.21) 
181 
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where: 
€; = relative dielectric constant of the insulating material; 
do = outer diameter of insulating medium excluding semi conductive screen 


or layer; 


d; => conductor diameter including semi conductive screen or layer, if any. 


2.3.4 Computation of g (Cable) 


The shunt conductance g (S/km) is the parameter taking into account the active 
power “dielectric” losses g - UB (per kilometre) given for each phase by the electric 
field in the insulating material [10, 13]. 

It is worth highlighting that the circuital element given to it (conductive shunting 
of capacity as in Figure 2.1) allows evaluating correctly the dielectric losses, but it 
not constitutes a realistic modelling in order to qualify the complex physical phe- 
nomena which are the fundamental cause and are referable to a kind of hysteresis 
in the dielectric polarization given by the alternating electric field at angular fre- 
quency w: only in a very small part the losses are referable to a proper Joule effect 
due to imperfect insulation, also present in DC. 

The shunt conductance can be computed with 


g=b-tand, (2.22) 
where 


tanéd = loss factor of insulating material (see Table 2.3), 


b=ac [S/km] = _ the per unit length capacitive susceptance of the cable, 


so highlighting active power losses proportional to capacitive power absorbed along 
the cable and hence proportional to the frequency and to the square of the voltage. 
The total (in the three phases) dielectric losses are given by 


waror = 3-g-Ug =g-U, = [W/km] 


where U,, [V] is the rated phase-to-phase voltage (see Table 6.4). 

The dielectric losses of HV and EHV cables ought to be always accounted for 
whilst they can be neglected for low and medium voltage since they are very small 
in comparison with the Joule losses in the longitudinal phase resistance. 

In any case, the IEC standard 60287-1-1 states that the dielectric losses should be 
taken into account for values of the phase-to-ground voltage Up equal to or greater 
than those of the Table 2.4. 
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Table 2.3 Values of relative permittivity €, and loss factor tan 5 at 50-60 Hz (see Section 2.3.4), 
for different insulating materials after [7] 


Type of cable Ey tand 


Cables insulated with impregnated paper 
Solid-type, fully impregnated, preimpregnated 4 0.01 
or mass-impregnated nondraining 
Oil-filled, self-contained 


up to Up = 36kV 3.6 0.0035 
up to Up = 87kV 3.6 0.0033 
up to Up = 160kV 3.5 0.0030 
up to Up = 220kV 3.5 0.0028 
Oil-pressure, pipe-type 3.7 0.0045 
Internal gas-pressure 3.4 0.0045 
External gas-pressure 3.6 0.0040 
Cables with other types of insulation 
Butyl rubber 4 0.050 
EPR — up to and including 18/30kV 3 0.020 
EPR — above 18/30kV 3 0.005 
PVC 8 0.1 
PE (HD and LD) 2.3 0.001 
XLPE up to and including 18/30(36) kV — unfilled 25 0.004 
XLPE above 18/30(36) kV — unfilled 2.5 0.001 
XLPE above 18/30(36) kV — filled 3 0.005 
Paper—polypropylene—-paper (PPL) > 63/110 kV 2.8 0.0014 


Table 2.4 Values of phase-to-ground voltage above (or equal to) which the dielectric losses should 
be considered [7] 


Type of cable Up [kV] 
Cables insulated with impregnated paper 
Solid-type 38 
Oil-filled and gas-pressure 63.5 
Cables with other types of insulation 
Butyl rubber 18 
EPR 63.5 
PVC 6 
PE (HD and LD) 127 
XLPE (unfilled) 127 
XLPE (filled) 63.5 


2.4 Computation of GIL Kilometric Parameters 


It has been already premised in Section 2.1 that the single-phase model at the pos- 
itive sequence of a GIL with solid-bonding connection represents a good approxi- 
mation: if the enclosures of a GIL (which are always very low resistive) are solid- 
bonded, they allow the circulation of induced currents almost equal (in magnitude) 
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(pitrendx td 


yn = kilometric resistance of the phase conductor 
gadx cdx Yen = kilometric resistance of the enclosure 


¢ = kilometric inductance of the coaxial circuit 
c => kilometric capacitance of the coaxial circuit 
g => kilometric conductance of the coaxial circuit 


a 


Figure 2.3 Single-phase model of a GIL element dx at positive sequence 


to those of phases and opposite in angle [14-17], so that they can be considered 
as return paths of each current in simple single phase circuits electromagnetically 
decoupled, having parameters computable on the basis of typical coaxial configura- 
tions (see Figure 2.3). 


2.4.1 Computation of GIL Apparent Kilometric Resistance r 


It is plain that, with the aforementioned hypothesis for the phase and enclosure 
current phasors, the total Joule losses must be computed by fixing (in the model of 
Figure 2.3) for the phase conductor an apparent kilometric resistance r equal to the 
sum of phase resistance rp, and enclosure resistance ren 


r=fontren [Q/km] . (2.23) 


It is worth mentioning that the thicknesses of the pipes and the used spacing render 
negligible the skin and proximity effects at power frequency [16]: so that rp, and ren 
can be evaluated in DC. The value of phase resistance rp is given by: 

Pph 


Dh ee [2/km] , (2.24) 


where pn is expressed in [Q mm/?/km] and the cross-sectional area is 
Sph = -(R>— Ri) [mm’] (2.25) 


with R» and R; the outer and inner radii of phase conductor expressed in millimetres 
(see Figure 2.4). 

Data concerning the involved materials and the resistivities can be obtained by 
Tables 3.4 and 3.5 or by more precise “manufacturer’s” information. 

Analogously for the computation of ren: 
Pen 
os [Q/km] , (2.26) 
Sen = 0 -(R{— Rj) [mm], (2.27) 


II 


Ten 


where R, and R3 are the outer and inner radii of the enclosure expressed in mm. 
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Figure 2.4 One phase and 
corresponding enclosure of 
a GIL 


The resistivities Ppn and Pen of material (Q mm? /km) must be evaluated at the 
operating temperature: the maximum allowable ones are stated in IEC 61640 [18]. 


2.4.2 Computation of GIL Kilometric Inductance £ 


By considering the typical coaxial structure of a GIL (see sketch in Figure 2.4) and 
by adopting the above-mentioned solid-bonding technique, if it is assumed that op- 
posite currents with magnitude / flow in the phase conductor and in the enclosure 
with uniform distribution (no skin and proximity effects), it is possible to compute 
the kilometric inductance as the sum of three terms £ = €, + €y + £- (the first rep- 
resents the phase internal inductance, the second the inductance due to the magnetic 
field between phase and enclosure, the third the enclosure internal inductance); it is 
also simple to ascertain that the value of €) is prevalent among the others and that, 
for the assumed hypothesis, the magnetic field external to the enclosure is zeroed. 

Since (see Figure 2.4) R,, Ro are the phase inner and outer radii respectively and 
R3, Raq the enclosure inner and outer radii respectively, it is easy, by integrating the 
magnetic field between R2 and R3, to express the inductance /, through 


Ho, R3 
£, = —In— fH : 22D 
ae a [H/m] , (2.28) 


where: 
}io => magnetic permeability of the free space = 4a x 10~-’ H/m. 


In the usual installations, it is In (#) = | in order to obtain, as is well-known, an 


optimal configuration of the electric field. 
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In order to obtain the expression of £, it is convenient to premise that (in accor- 
dance with electromagnetic theory) the following relation can be written 


en, = / [(u0H7) /2]dv = £,17/2; (2.29) 


it correlates the magnetic energy ena, expressed by the integral computed in the 
whole volume of the phase conductor (of 1 m length), with the corresponding quan- 
tity £,17/2. 

Once computed eng, it yields immediately 


£, = 2en,/I? . (2.30) 


It is sufficient to observe that the magnitude of H, in the integration interval Rj—R2 
can be written (since the circular symmetry and the hypothesis of uniform current 
density op, in the phase) by means of 


Opn «1+ (x7 — RP) I 


i= Qnx sd ca (R3 — Ri) 


(Ri <x <R). 


Analogously, the computation of £, can be performed by writing for the magnitude 
of H,, in the integration interval R3—R,4 the two relations (note that the current —/ 
flows in the enclosure) 


I en’ 4 2_ R?2 I 
4 Cent (x 5) 


5 On =——T CCR SX Ra). 
m (Ri — R3) 


~ Onx 2x 


By developing the above-mentioned procedure, the following formulae (see also [19, 
20]) can be achieved 


= | Ry (2) + Fe = (H/m) , (2.31) 
’) 


den | (RR AR 


_ Ho Ri (®) (R3 —3R} 
See =e Be) a@geey | e ) 


In order to have some orders of magnitude, it seems useful this example. Once the 
following values are assumed: 


R,; =80mm, Ro=90mm, R3=250mm, Ry = 260mm, 
by means of (2.31), (2.28) and (2.32) and by converting from H/m to mH/km the 


following values yield: 
£, = 0.0074, ¢, =0.20433, = 0.0026 mH/km, 
L=l,+ 4+. = 0.2144 mH/km , 


where it appears that £, is the most important one and, as first approach, it always 
gives sufficient accuracy. 
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2.4.3 The Computation of GIL Kilometric Capacitance c 


The value of capacitance can be easily computed considering the cylindrical sym- 
metry and that the permittivity of insulating gas mixture is ¢9 = 8.8542 107!? F/m 
(as for the free space): 


(re. Ra St ee a (2.33) 


in (#2) 18 -In ( #) 


2.4.4 Computation of GIL Kilometric Shunt Conductance g 


It is worth noting that the dielectric used in the GIL (pure pressurized SF¢ or a mix- 
ture of it with No, at about 7 bar for U = 380 kV) has always a very low loss factor 
and gives losses that result absolutely negligible also if compared with the small 
Joule longitudinal losses due to the low apparent resistance r expressed by (2.23). 
Also in precise economic evaluations of losses, it is always assumed g = 0. 


2.5 Some Other Matrix Relations Deriving 
from the Fundamental One 


In the transmission equations (2.13) and (2.14), there are four variables: the voltage 
and current phasors at sending and receiving ends. 

When two of the four variables are fixed, the remaining ones are univocally de- 
termined. 

Consequently, in this paragraph, different ways of using the transmission equa- 
tions are highlighted in order to both give the reader a test-bed for getting some 
practice and have suitable formulae for Chapters 3 and 4. 

As a first example (2.34), as much typical as (2.16), can be derived by rearrang- 
ing (2.13) and (2.14): it expresses the vector wp (hybrid), composed of the voltage 
and current at receiving-end, as a function of the vector wy, at sending-end: 


Gor |_| D -B] | Yos |. Sagal 
tlle a | ee ( Ma=M'. ven 
KY a ey 

WR M, Ws 


It is easy to demonstrate that M_, is the inverse of M (or that the product of M and 
M, is the identity matrix), by remembering that AD — BC = 1 in virtue of the 
reciprocity principle (see (2.17)). 

As a second example, if the sending-end voltage phasor and the receiving-end 
current phasor are fixed, it is easy to obtain the matrix relation (with D = A) 


U, 1/A -B/A] [U, 
oa = Ee 1/A | rt (2.35) 
a 


M, 


which allows determining the other two phasors Upp and J. 
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As it will be used in Chapter 3, in order to know the variables along the trans- 
mission line, it is rather immediate to achieve from (2.34): 


U D, —B U 

OX | aX =x |,}| =0S 
le 

Wy M,. Ws 


which allows calculating the vector w, along the cable (at a distance x from S) 
having A, .... D* obvious meaning (once known the vector w,). 

It will be often sufficient to compute only the values of voltage Uy, and current 
I, at line midpoint 7 with (2.37) 


Z| Dr | ae | 
= : 2.37 
Ir -Cr Ar Is eon 
KS SS 

wr Mr Ws 


by giving x the distance of T from S. 

In Chapter 4 other relations, besides those of this chapter, will be developed 
bringing again reasons to appreciate the matrix algorithms [21, 22]. 

Furthermore, it is worth remembering the matrix relation: 


Ups |__ [ A/C. -1/C.] FL 
[es ]- [42 we] [E] i 
SE 


Zz 


which involves an impedance matrix Z (where for the current J p at port R as well 
the same phasor references of Figure 2.1 are held). 

In order to deduce the matrix relation in (2.38), (2.16) can be made explicit so 
that, by setting J p = 0: 
UR_! 


Us _A =Z d = —=Z 
=== an = : 
I Ga Se Ig (Cr 


By making (2.34) explicit, with J, = 0, it easily yields: 


Up D Us 1 
i= CE ~ Err and Te G7 ESR: 


2.6 Cascade Connections of Two Port Networks (TPN) 


Let us suppose to have the cascade connections of Figures 2.5 and 2.6. 
As is well-known, the resulting transmission matrix is given by the product of 
the transmission matrices. 
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Figure 2.5 Cascade of 2 TPN S ®@ H @ K 
A B, A, B, 
ee eel. 
S ® H @ K ©) R 
A; B, A; B, A; B; 
C D, C> D) peas C3 D; 


Figure 2.6 Cascade of 3 TPN 


For each section ©, @ and @ the corresponding transmission matrices can be 
calculated as already thoroughly explained in the previous sections: 


A; |B, A>| By A3|B3 
Ci|Dil? |G |Dal’ |Cs|Ds|- (2.39) 
M, M, M, 


The equivalent cascade matrices are given (with evident symbols) 


Ap (Bo Ajp3 | Bias 
M,, = M,M,= >; M);, = M,M,M, = - (2.40) 
Cy | Di C173 |D123 


For example, by considering the ports in Figure 2.6, the following expressions ap- 
pear immediately understandable: 


‘ —1 ; 
ws =M,- wy; Wr= Mz -we: 
ws =M,-weK: wy = M,3-WR: 

= . _ yy-l 
ws = M3: wR: wr=M),-ws. 


It is worth noting that generally even if D,; = A,, D, = A>, D; = A, (see (2.15)), 
the same is not necessarily valid in a cascade matrix where for example it is (but in 
particular cases) D,, # Aj, and Dj, A A493. 

In any case, as long as non-reciprocal two port circuits are introduced (e.g. 
a phase shifter two-port), relations such as A153 - D493 — By23-C 493 = 1 remain still 
valid. 


2.7 Parallel Connection of Equal Two Port Circuits Thermally 
and Electrically Decoupled 


In case of a transmission link of a double circuit of equal cables in electric parallel 
(a situation that can be useful to gain a great power capability and enhance the 
reliability at the same time), the equivalent matrix can be easily deduced on the 
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za 
2 


Figure 2.7 The parallel of two, perfectly decoupled, equal circuits 


basis of the obvious equations in (2.41) valid for each of the two circuits (which 
are supposed to be thermally and electrically decoupled) with equal sharing of the 
currents shown in the scheme of Figure 2.7. Of course, it is true if the port voltages 
have the same reference (thick links in the sketch) as it is always valid in the positive 
and negative sequence models. 

Therefore it is possible to determine the equations in (2.42) and hence the equiv- 
alent transmission matrix M_, of the double circuit, once computed the single circuit 
matrix M, 


Uns = (A)- Ung + (B)-, 2.41 
po | — 
Uos = (A) Vor + (4) -Le. | (2.42) 
Tg = (2C)-Uor + (D)- LR. 


IQ | is 
IS | le 
NO 
a 
Is 


M M 
It seems interesting to note that the same result can be achieved (following a more 
intuitive procedure) by giving the double circuit the kilometric parameters z/2 and 


2y (since z and y are those of the single circuit) so having the same propagation 
factor and the halving of the characteristic impedance. 
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2.8 The Shunt Reactive Compensation 


It is known that the kilometric cable capacitive susceptance has often a so high 
value to determine (chiefly with great lengths, voltages and cross-sections) steady 
state and transient regimes which could become very severe or even unacceptable 
(as it will be shown in Chapters 3 and 4), unless a suitable shunt compensation is 
adopted. 

Usually the cable link modelling must be completed with the insertion of the nec- 
essary single-phase shunt reactors at the two cable ends and in case (depending upon 
the link length) also in some other intermediate locations (lumped compensation). 

By considering for example the circuital scheme of Figure 2.8, the steady state 
regime at angular frequency w of a short link compensated only at the ends can be 
studied by means of the matrix cascade M,-M-M z, being that M is the cable trans- 
mission matrix and M_, is the transmission matrix which qualifies the behaviour of 
a shunt reactor with inductive susceptance equal to Y, (as it will be shown in Sec- 
tion 2.8.2); it remains a useful basis also for the switching transients studies, since 
the subtransient regime at the same angular frequency w (e.g. Sections 3.7 and 4.7) 
has a precise “indicial” role for the transient evaluations. 


Figure 2.8 Lumped shunt 
compensation of short UGC 
(d < 20-25 km) 


By considering that the lumped shunt reactors (located at short mutual distances 
< 20-25 km) can be numerous and can render the diagrams and the analytical re- 
lations heavy, it appears to be more efficient and meaningful, as a first approach, to 
imagine alternatively a uniformly distributed shunt compensation (see Section 2.8.1) 
in a way consistent with the uniformly distribution of the cable capacitance: it al- 
lows easily computing (and with sufficient accuracy) the necessary shunt compen- 
sation degree and having an equivalent transmission matrix which, even if ideal, is 
extremely useful. 


2.8.1 The Uniformly Distributed Compensation 


The compensation of the uniformly distributed cable capacitive susceptance would 

theoretically require the insertion of uniformly distributed inductive susceptance. 
Since this is hard to imagine constructively, it is possible to think to insert, in 

each phase and at each kilometre, a reactor of inductive susceptance of finite value 
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Figure 2.9 Kilometric shunt 
branch in the hypothesis of 
uniformly distributed shunt 
compensation 


CABLE | SHUNT REACTOR 


b tand bésn Pp 


b= ac [S/km] 


equal to —j &nhwc = —j &nb (see Figure 2.9), which compensates at each kilometre 
a fraction &, of the cable capacitive susceptance, bringing to compute the resulting 
kilometric susceptance y after (2.43) 


y=b-tand+b-&n-pt+jb-U—§n) ; (2.43) 


where the reactor “loss factor” p (which allows considering the reactor power 
losses, always very low) has been introduced. 
In case of absence of compensation it has &, = 0, and (2.43) returns the usual 
one: 
y =b-(tanéd +j). (2.44) 


Therefore, once fixed the shunt compensation degree &, and evaluated the suitable 
kilometric admittance after (2.43), it is possible to continue directly, by means of the 
already detailed relations, with the computation of all the new transmission param- 
eters (k, Z,, and so on) up to the uniformly compensated link transmission matrix. 


2.8.2 The Lumped Compensation 


The transmission matrix obtained for a cable link uniformly compensated (as de- 
scribed in Section 2.8.1) allows developing (as it will be shown) quick procedures 
for a precise choice of the compensation degree &: once fixed such value the link 
planning must be completed with an accurate configuration of the corresponding 
lumped compensation and with the check of it by means of a specific circuital anal- 
ysis. 

In the example of Figure 2.10, it is evident that the analysis regards the cascade 
of the three matrices M, - M_- M;: it will be sufficient to form each matrix M¢ 
(regarding a reactor of inductive admittance Y;) after (2.45): 


M, = 1/0) (2.45) 


AG 


which has an immediate confirmation in the circuit and in the obvious relations of 
Figure 2.11. 
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Cable link 


(with d<20-25 km) Ye Yerj we bon 


IS 


Ms 
Figure 2.10 Cable link compensated at both ends 


U,=1- U, +0-I, 
1,=Ye U,+1-L 


Figure 2.11 Shunt reactor 
circuit 


If the cable to be compensated has a length greater than 20-25 km it is suitable 
to install shunt reactors also in m intermediate sections, respecting the limits shown 
in Figure 2.12. 

The admittance value are subsequently modified in Y’. = 2Y; =-j océnd/n; 
it is almost trivial to observe that each cable section of length d/n must be modelled 
by its transmission matrix without compensation. 

Once formed the cascade of the different two port circuits, the regime at angular 
frequency w at each port can be analysed by making suitable use of the matricial 
algorithms already described in this chapter and of the procedures which will be 
exposed in the Chapters 3 and 4. 

Hence it will be possible not only to ascertain the effectiveness of the hypothesis 
of distributed compensation, but also to evaluate in detail the effects, more or less 
positive, due to possible alternative schemes. 

Usually the reactor percentage losses range between 0.1% and 0.3%. 


a 
n 


< 20-25 km 


Figure 2.12 Lumped compensation of a cable link 
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Chapter 3 


Operating Capability of Long AC EHV Power 
Cables 


3.1 Introduction 


The operating constraints, regarding current and voltage regimes, involved in power 
transmission of long AC cable links are presented and here throughout developed. 
The possibility of transmitting high power rating is examined and the transmission 
length limits are evaluated. The procedure, which considers positive sequence mod- 
els (as in Chapter 2), is essentially based on the classical transmission equations and 
their pertaining diagrams. 

This chapter has been particularly devoted to XLPE cables and gas insulated 
lines, but the procedure is worth applying to any distributed-parameter transmission 
line included overhead lines. The capability charts are a useful tool to highlight 
the operating characteristics of a cable link and a guide to evaluate the results of 
different degrees of reactive compensation. The results can offer the transmission 
system operators very effective means in order to evaluate the existing and future 
underground cable links to achieve high performances and a long lifetime [1]. 


3.2 The Basic Constraints 


It is well known that the lifetime of an AC cable with solid insulation strictly de- 
pends upon the story of the thermal (temperatures and their durations) and dielectric 
stresses suffered by the insulation. Therefore, a long service life can be reached with 
these constraints in any point along the cable: 


¢ avoiding that the phase currents overcome the cable steady state ampacity /, 
(even if, in many cases, the cyclic and emergency current regimes occur) meant 
as the continuous constant current just sufficient to produce asymptotically the 
maximum conductor temperature with acceptable dielectric degradation; 

¢ avoiding that the phase-to-earth steady state voltages |U)| overcome the highest 
RMS Voltage U//3 (see the Insulation Co-ordination International Standard 
which also reports other limits for overvoltages of different duration) for cable 
both on-load and no-load conditions. 
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It is evident that foreseeing any conditions (steady state, transient and cyclic 
regimes) which the cable is subjected to during the operation of the network is 
rather difficult, if not after exhaustive power flow studies and transient simulations; 
notwithstanding, in order to frame in a preliminary way the possible performances, 
reasonable conservative conditions are generally: 


¢ the power frequency currents at both ends (S for sending-end and R for receiving- 
end), during operation, must not exceed a current magnitude J, < J,; 

e the phase-to-earth voltage level in one of the two ends must be fixed (in this book 
S will be always chosen). 


Therefore, the following “basic constraints” are fixed: 


Frl<k., (3.1) 
Isl <ke, (3.2) 
[os | = Voc - (3.3) 


In the following investigations, where U, = 420KkV, the level Upp = 230kV = 
95% Um/+/3 has been always assumed. 

A thorough analysis of the cable operating characteristics with the aforemen- 
tioned constraints has been presented in [2], which constitutes (chiefly by means 
of a general abacus) one of the first elegant (dated 1986) approaches: it has been 
a strong suggestion for the authors. 

The cable steady state ampacities in undergrounded installations are evaluated 
according to the thermal study affected by: 


e Joule losses in the core conductor and the sheath (and armour if any); 

e dielectric losses; 

e cable spacing and burial depth; 

e thermal equivalent resistances of the cable (e.g. between conductor and sheath 
and between sheath and armour); 

e thermal resistance of the surrounding medium depending upon the soil resistivity 
and mutual thermal resistances between the cables; 

e ambient (soil) temperature; 

¢ allowable temperature of the dielectric. 


Table 3.1 Data of EHV single-core cables (380 kV, 50 Hz) 


Cable #a #b #c #d 
Cross-sectional area mm 1600 Cu 2500 Cu 3250 Cu 630 Cu 
Conductor diameter mm 50.7 63.4 AZ3 30.8 
(Milliken type) 

XLPE insulation diameter mm 103.7 119.9 130.4 92.8 
Metallic sheath diameter mm 114.1 Al 130.1 Al 140.4 Al 101.8 Lead 
Sheath cross-sectional area mm = 500 = 500 = 500 = 1000 
PE jacket diameter mm 125.7 141.7 152.0 117.0 


Total mass kg/m 26 37 45 28 
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The Standard IEC 60287 [3], regarding the calculation of the ampacity in differ- 
ent conditions, is now completely acknowledged by almost all the National Rules; 
numerous scientific books are devoted to UGC constructional and installation (elec- 
trical and thermal) features [4-6]. These contributions show very clearly that in tran- 
sient or cyclic operation the thermal capacities (of the conductor, of the dielectric 
and of the surrounding soil) can temporarily increase the ampacity. 

A systematic analysis of the steady state regimes compatible with the aforemen- 
tioned constraints (3.1) and (3.2) can be developed by means of the transmission 
equations (2.13) and (2.14) of the “positive sequence model” (where D = A): 


Ups =A-Unrt B-Lp (3.4) 
Le= Cs Upp td Aol nc (3.5) 


Tables 3.1 and 3.2 report the main data for the four different EHV (380kV AC) #a, 
#b, #c and #d cables. For a wide survey of the installations and a general frame on 
the EHV cable technology refer to Chapter | and [7]. 

In order to show the application of the procedures, examples regarding chiefly 
the use of cable #b will be presented. 

The per unit-length capacitive susceptance of UCG (see examples of Table 3.2) 
is about 15—20 times that of a typical 380 kV OHL and can require (depending upon 
length d of the cable) strong capacitive power: the reactive susceptance affects sen- 
sibly the power flows in steady state regimes and plays furthermore a key role on the 
transient regimes of the line at no-load. A rough estimation brings to compute, for 
UGC #b, the heavy steady state capacitive power O, © (3UZwc) = 11 Mvar/km 
at 50 Hz and 13 Mvar/km at 60 Hz: it becomes foreseeable that shunt compensation 
reactors are almost always necessary to achieve good effects both in the steady state 
and transients regimes (see Section 3.7). It is worth noting that at f = 60Hz the 
ampacities are lower than at f = 50 Hz due to the greater power losses (both Joule 
and dielectric ones). 

It is worth remembering that the hypothesis (made in Chapter 2) of uniformly 
distributed shunt compensation can be adopted to easily compute the suitable trans- 
mission parameters A, B and C to use in (3.4), (3.5) and in the other equations 
involved: in fact, it has been throughout verified that both the steady state and sub- 
transient regimes remain almost unchanged by applying (with the same &,,) lumped 
shunt compensations on length intervals not greater than 15—20km. 


3.3 First Analysis: U,, (5), 1, Constrained 


A first analysis on the possible operating conditions can be performed if in (3.4) 
and (3.5) both the current phasor Jp at R and the voltage phasor Up, at S are 
constrained as follows, in accordance with (3.1) and (3.3), 


Ip = (1.20) on the real axis, (3.6) 
Uos = (Ue 26), 5 = 0-27, (3.7) 
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so that the remaining variables Upp and J, are univocally individuated (see Sec- 
tion 2.5) by means of: 


1 B 1 B 
Vor = A Yos— GER > VYor= ee a 
(3.8) 
C 1 C 1 
ts= | Yost Gtr eu eT Mc a 
(3.9) 


In order to develop the analysis of the J , regimes determined by (3.9), the phasorial 
diagram of Figure 3.1 is useful (where only for graphical purpose, point N’ has been 
placed sensibly far from point NV): it shows, once fixed on the real axis 1p = I, = 
1.78kA (if cable #b is considered), how the angle 6 of Uys = (230kVZ6) can 
range the whole interval 0-27 giving rise to two different groups: 


i) any regime with |J.,| < J, = || compatible with the cable ampacity ac- 
cording to the constraint (3.2) (e.g. see the phasors OL), OL2, OL3); 

ii) any regime with |/.,| > J, not allowable (e.g. see the phasors OH), OH2, 
OH). 


In particular, the phasor OK; (when 6 = 6,) and the phasor OK2 (when 6 = 43) 
correspond to |J.,| = J, = |Lp|, and concern two very meaningful regimes for 
the power flows between S and R (see Section 3.6). It is evident that for 6 ranging 
between 4; and 52 the voltage phasors Uy, (5) determines in this first analysis all 
the current phasors J, (e.g. OL), OL2, OL3) having magnitudes within the limit J, 
according to (3.2). 


(Briggs) 


Figure 3.1 Phasorial diagram of first analysis 
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The values 6, and 52 can be determined by the following relations 


65 =a-y-f,, (3.10) 
by = 6; + 281, (3.11) 


which are evident by observing the Figure 3.1, where the necessary Brigg’s formu- 
lae are also reported. It is worth including in the analysis procedures a “warning 
message” to display possible incompatibilities of the magnitudes |a|, |b|, |c| in the 
formation of the corresponding triangle. 

The calculation of Upp by means of (3.8) closes the first analysis: after this, a first 
set of steady state regimes compatible with the constraints (3.1)—(3.3) is completely 
individuated in the interval 5;—6>. 

For these regimes (where |/ p| = J, and |1.5| < Jc), it is of paramount interest 
to compute the complex three-phase powers at ports S and R. 


3.4 Second Analysis: U,, (3), 1, Constrained 


A second analysis on the possible operating conditions can be performed if in (3.4) 
and (3.5) both the current phasor J , and the voltage phasor Uy, at S are constrained 
as follows, in accordance with (3.2) and (3.3) 


Is = (40) on the real axis , (3.12) 
Uos = (Ue 240), 8 = 0-27, (3.13) 


so that now the remaining variables Upp and J p are univocally individuated (see 
(2.34) in Section 2.5) by means of: 


Upp =D-Ugg-B-ly > Upp =D-(230KVZ0) — B-(I-Z0), 
(3.14) 


[pe eC Up tay <1 pS Cr C0294 AE): 
(3.15) 


Equation 3.15 can be visualized by means of the phasorial diagram of Figure 3.2. 

This diagram shows how the angle } of Ups = (UpcZ0), with respect to 
Is = 1, on the real axis, can range the whole interval 0-27 giving rise to two 
different groups: 


iii) any regimes with |J.p| < J, compatible with the cable ampacity (e.g. OL’, 
OL!" OL"); 

iv) any regimes with |J p| > J, not allowable (e.g. OH’, OH”, OH"). 

In particular, the phasor OK’ (when 3 = 0’) and OK” (when 3 = ”) individuate 


again two regimes with |J.p| = J, = |J.g5| and with the same relations between 
current angles (as OK, OK; in the first analysis). It is evident that for ? ranging 
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C=C; @-C=C Znty ; 
c =-CU,s (8") 
(mt+7+ O')+ B- a= 


(Briggs) 


Figure 3.2 Phasorial diagram of the second analysis 


between %” and v” the voltage phasors Uy, (v) determines in this second analysis 
all the current phasors J p (e.g. OL’, OL”, OL") having magnitudes within the limit 
I, according to (3.1). 

The values 0” and #’ can be determined by the following relations 


0” =-By-yta, (3.16) 
oY = 9" 426, (3.17) 


which are evident by observing the Figure 3.2. 

The calculation of Upp by means of (3.14) closes the second analysis: conse- 
quently a second set of steady state regimes compatible with the constraints (3.1)— 
(3.3) is completely individuated. 

Also for these regimes (where |J.5| = J, and |Ip| < IJ), it is of paramount 
interest to compute the complex three-phase powers at ports S and R. 


3.5 Voltages and Currents Along the Cable 


The analysis of the currents and voltages at sending-end S and receiving-end R 
does not assure “‘a priori” a satisfying electric behaviour along the cable: in fact, it is 
important to examine the currents and voltages along the cable in order to investigate 
possible points exceeding the fixed J, or the voltage limits; to this aim, the use 
of (2.36) could be suitable. 

For instance in Figures 3.3—3.6 the voltage and current behaviours are shown 
along the 60 km cable #b in the cases &, = 0 and & = 0.608 (distributed compen- 
sation); the regimes of the “first analysis” are pointed out in Figure 3.1 by means 
of the points K;, L;, L2, L3, Kz (equidistant between themselves on the circum- 
ference arc) and the regimes of “second analysis” are pointed out in Figure 3.2 by 
means of the points K”, L’, L”, L'”, K’. 
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Figure 3.3 Current magnitudes along the cable #b (€,, = 0.0) 
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Figure 3.4 Voltage magnitudes along the cable #b (&,, = 0.0) 
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Figure 3.5 Current magnitudes along the cable #b (&,, = 0.608) 
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Figure 3.6 Voltage magnitudes along the cable #b (€, = 0.608) 


It can be seen how the current levels are always within the ampacity limit 
(~ 1.8kA) and how possible voltage levels (steady state) meaningfully high can 
be detected at the end R. 
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To this aim, it is worth noting, as a typical example, the regime L”, which would 
give rise to (see Figure 3.6) the high level Upp ~ 247kV: however it would have 
(see Figure 3.5) currents |J.5| = 1.8kA and |Z p| + 1.4kA and a corresponding (see 
Figure 3.31) strongly capacitive power flows (~ —1000 Mvar) of scarce interest in 
the network operation. In general, strong variations of current magnitude along the 
cable are justifiable when the quadrature components are prevailing. 

A more complete and elegant method for the analysis along the cable is described 
in Section 3.9.2. 


3.6 Power Values Compatible with Basic Constraints 
and with Voltage Levels at the Receiving-End 


The procedures shown for the first and second analysis allow individuating all the 
regimes at power frequency compatible with the basic constraints (3.1)—(3.3): by 
using the couple of phasors Ups, J; and Upp, Ip in each steady state regime the 
computation of the three-phase complex power Sig and Sp (at ports S and R) is 
immediate for each given line. 

The corresponding graphical representation is given in the following Figures 3.8— 
3.11 which refer to cable #b without reactive compensation (& = 0). 

By using the criteria of first analysis, the regimes (in the interval 6;—-d2) compat- 
ible with the constraints |J.p| = J, and |I.5| < J, have been evaluated and their 
receiving-end complex power have been reported forming the curve (a), whilst the 
corresponding sending-end complex power forming the curve (b). 

Analogously, by analysing (in the interval 3”-0’) the regimes |J5| = J. 
and |Ip| < J,, according to the second analysis, the curves (c) and (d) have 
been drawn; the curve (c) represents the receiving-end complex power and (d) 
the sending-end complex power. Therefore, once fixed J, = 1.78kA = I, = 
steady state cable ampacity, the contours (a) and (c) set the boundaries of allow- 
able “receiving-end power area” (see, for instance, grey area in Figure 3.9) whereas 
the contours (b) and (d) of allowable “sending-end power area” (clearly reported in 
the examples of Figures 3.8—3.11). 

In particular, the K points, which highlight “maximum” values for Pr and Ps 
correspond to the homologous regimes of Figures 3.1 and 3.2. Any complex power 
Sp within the receiving-area gives, at ports S and R, current magnitudes < J,; the 
same is valid for any S 5 within the sending-area. 

Since the voltage levels at node R play a key role in the network operation, it 
seems suitable to give them a direct visualization in the “receiving-end power area’, 
by means of phase voltage curves parameterized with |Uyp| constant. It is sufficient 
to implement the well-known expression of receiving-end power 


Uys —A-U. 7 
Sn =3-Usg | SF Bow) a Es) ; (3.18) 
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Figure 3.7 Power balance of a UGC transmission line 


Figure 3.8 Power flows 
compatible with the basic 
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once fixed Ujy = 230kVZo and by setting Upp Zp (e.g. with magnitudes 220, 
225, 230, 235, 240 kV) it is possible to give suitable o—p values such to determine 
complex power S'p in the “receiving-end area”. 

In such a way, it can be clearly singled out the regimes which are not acceptable 
owing to excessively high or low voltage levels at R (even if compatible with the 
basic constraints): for example, heavy reactive power flows or active ones on long 
runs. In particular the voltage curve intersecting the origin of the axes also points 
out the steady state voltage in receiving-end at no-load regime (Pr +jQr = 0), 
when |Ups| = 230kV. 

By the relation (3.19) of the cable power balance (see Figure 3.7) it can be ob- 
served that each receiving-end power Pr differs slightly from the corresponding 
sending-end power Ps, since the active power losses p are always low enough, also 
for UGC long runs. 

Equation 3.20 shows that each receiving-end reactive power QR results clearly 
higher than the corresponding Os, since the uncompensated cable itself “generates” 
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Figure 3.9 Power flows 
compatible with the basic 
constraints (without shunt 
compensation) 
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Figure 3.10 Power flows 
compatible with the basic 
constraints (without shunt 
compensation) 
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Figure 3.11 Power flows 
compatible with the basic 
constraints (without shunt 
compensation) 
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high inductive power q, increasing with the length d. That helps to explain how, by 
increasing the length d, the area contours have progressive translations along the 
imaginary axis of the complex plane. By considering the cases from d = 30km to 
120 km (without compensation), it can be ascertained that the receiving and send- 
ing areas are progressively “moving away” from the x-axis (beyond a perceptible 
“shrinkage” along the y-axis): this highlights that in the power flows a higher and 
higher reactive power (rather than an active one) must be accepted, so opposing the 
chief aim of the EHV transmission links. 

By this important point of view, the length of 25-30km (see Figure 3.8) can 
be considered as a limit for the cable #b unless the shunt reactive compensation is 
performed with a given degree (see Section 3.7). 

Itis worth adding that, beyond the aforementioned drawbacks, in order to have suit- 
able power flows, other heavy problems, regarding the system security itself, are the 
regimes (steady state and transient) of no-load cable energization and de-energization 
besides the load-rejection, which almost always force the use (see Section 3.7) of 
a suitable shunt compensation and give the installation its final arrangement. 


3.7 No-Load Energization and De-Energization 


It is worth remembering that in the network operation, very frequently recurring 
events of greatest concern are the energization (and de-energization) of a no-load 
line since it is always necessary to prepare the operating structure of the grid. 
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Figure 3.12 No-load energization at port S 


The sketch of Figure 3.12 shows the switch-on of the circuit breaker in S in order 
to energize the cable and serves as a first, not complete but meaningful, approach 
devoted to highlight the limit conditions for line energization at no-load. In a generic 
network, the power supply in S can be modelled as an equivalent generator which 
is characterized by its electromotive force Up (supposed to be equal to 230kV) 
and the short-circuit subtransient impedance Z” (for simplicity purely inductive 
Zz =j KY, 

The magnitude of no-load power frequency sub-transient voltage Ujp at R, due 
to the closing at S, constitutes a reference (“indicial” value) to foresee the peak 
value of switching surges due to subsequent evolution of the phenomena and is 
completely defined through the following formulae: 


Uo A 
"” e= 
Ujs =(——z_ | - > (3.21) 
(= hg Cc 
Os 
0, 
U 
—_ : (3.23) 


Yee AE ZG 


with A/C the impedance (almost completely capacitive) as seen from S with R at 
no-load, and expressing with (3.22) the Ferranti’s effect. 

With regard to Z” = j X” evaluation, it is possible to refer to the subtransient 
impedance Up/I.” (from network studies): since the values of subtransient short 
circuit current 1/” (three-phase at S) in EHV networks can be foreseen in the range 
10-50kA, X” corresponds to 23-4.6 Q. 

In order to respect the standard switching levels (e.g. 1050 kV) with a conserva- 
tive margin, it seems advisable that the phasor Uj, does not exceed the magnitude 
Um/ V3 = 242.5 kV: it is easy to ascertain that such target can be reached introduc- 
ing by attempts in (3.23) a suitable compensation degree which modifies suitably 
the parameters A and C. 

In the hypothesis that, extinguished the transient phenomena, the voltage regula- 
tion restores again at port S' the rated value Up = 230kYV, it is possible to compute 

Uo 
ha. = GG (3.24) 
which expresses the no-load steady state (almost entirely capacitive) current Jy) 
(which the circuit breaker b must interrupt in case of de-energization of the no-load 
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cable): for this current the Standards on the Circuit Breaker (Section 4.107 of [8]) 
suggest the limit value of 400 A. 

After all, both the following additional constraints must be fulfilled (by assigning 
by attempts a suitable compensation degree &) 


[Ul n| < 242.5kV, (3.25) 
|Zn_| < 400A. (3.26) 


The curves of Figures 3.13 and 3.14 show that the constraint (3.26) is almost al- 
ways decisive for the determination of &, unless there is an agreement with the 
manufacturer for a circuit breaker with higher /yp. 

Figure 3.15 clearly defines the values of shunt reactive compensation degree &., 
in order to fulfil (3.25) or (3.26) showing which one is the more limiting criterion as 
a function of the line length. 

It appears almost trivial to note that analogous considerations in order to verify 
the additional constraints must be made also for the energization by the port R 
(besides by the port S'), by introducing for the parameters Up and X” suitable values 
in consideration of the supply grid linked to R: the detection and the choice of the 
“best-end switching” [9] are mentioned by CEGB as effective experienced practices 
in network operations. 

Moreover, the cases where the steady state capacitive power absorbed by the 
cables (QnL = 3Uo/NL) exceeds the ability of synchronous generators (dangerous 
self-excitation conditions), located in close proximity of the cable installation, must 
be avoided [10-12]. 

In any case it is absolutely advisable that the Transmission System Operator 
(TSO), when planning a new link, performs both detailed power flows and network 
simulations. 


sunanm due to (3.25) 


mms uc t0 (3.26) 
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Figure 3.13 Limit lengths due to the constraints (3.25) and (3.26) for cable #b 
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Figure 3.14 Limit lengths due to the constraints (3.25) and (3.26) for cable #d 
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Figure 3.15 Reactive compensation degree as a function of cable length due to the con- 
straints (3.25) and (3.26) for cable #b and X” = 20Q 


3.8 Power Capability Charts 


Once the degree of shunt reactive compensation has been chosen (in accordance 
with Section 3.7) in order to fulfil the additional constraints (3.25) and (3.26), it is 
possible to draw the final capability charts compatible with the basic constraints, 
showing the voltage levels Ur at R for a good line operation, reporting above 
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a summing-up table with the values of the following parameters which fully charac- 
terize the installation and with some noteworthy regimes: 


d cable length; 

I. chosen current < /,; 

x" 50 Hz subtransient reactance (seen at port S); 

Eon shunt compensation degree according to Section 3.7 (unif. distributed one); 


UsR RMS value of the subtransient voltage in the no-load cable energization 
from the port S, with a given degree &.; 

Ixy RMS value of the current absorbed in S by the cable at no-load in R, in the 
steady state regime with |U) | = 230kV and with the given degree &,; 

Qn. Capacitive reactive power absorbed in S by the cable at no-load in R in the 
steady state regime with |Uy;,| = 230kV and with the given &; 

Pri receiving-power with cosg@ = 1 (indicated with U in the figures 3.16— 
3.27); 

AP, active power losses at PR}; 

AQ, _ whole reactive power required by the cable at PR. 


The so-depicted “capability charts” (see the following examples) are an extremely 
useful tool to quantify and compare, in a visually immediate way, different planning 
solutions. 

The capability charts of Figure 3.16 show how for the cable #b at 50 Hz the shunt 
reactive compensation can be avoided only for lengths inferior to 23.6km, if the 
capacitive no-load current Jy, must not exceed 400 A (with Ups = 230kV); Fig- 
ure 3.18 highlights (in the table above) that with d = 30km the condition &, = 0 
would maintain |U{'p |' within satisfying limits, but would give rise to Jy, = 509 A. 

Figure 3.19 shows how the choice of &, = 0.214 can give the link the final set- 
tlement: it is of note (beyond the limitation of /N_ to the suitable value 400 A) also 
a slight increasing of Pr; and a good subtransient voltage level, which leaves to 
foresee for the allowable limit subtransient reactance values up to 27 Q, in accor- 
dance with (3.23) and (3.25). Figure 3.17 shows the capability charts of the cable #b 
at 60 Hz without shunt compensation: it is worth noting the worsening of the limit 
length (d = 19.7km) and of the ampacity. In general, the frequency (60 vs 50 Hz) 
plays its important role in the UGC analyses. 

If other lengths of the cable (see examples with 60 km in Figures 3.20 and 3.21; 
90 km; 120km) are considered, it can be observed, with immediate visual compar- 
isons, that the high benefits given by the shunt reactive compensation (after Sec- 
tion 3.7) involve not only the energization — de-energization phenomena, but also 
the power flows. 

It is worth underlining the case with d = 90 km, (Figures 3.22 and 3.23) where 
the compensation degree &, = 0.74 (which gives additional costs) modifies dras- 
tically the no-load phenomena Uj, and yi, and requalifies dramatically the al- 
lowable power flows: it can be noted the power Pr, increase from 652.5 MW to 


' Tn all the figures, for the sake of brevity, |U/ ORI will be written as Usp. 
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Figure 3.16 Capability charts 
without shunt compensation 
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1220 MW (87%), the reduction of AQ, and, in general, improvements of the power 
flows in the four quadrants. 

In the case of d = 120km (always cable #b) the comparison of Figures 3.25 
and 3.24 appears very impressive at first sight, because there is an almost com- 
plete “overlapping” (beyond a sensible widening and a good symmetrization in the 
complex plane) of the “receiving” and “sending” areas, by means of the very high 
Es = 0.81, for which (differently from the foregoing cases) the constraint (3.25) 
is more influencing rather than (3.26). The shunt compensation degree &, = 0.810 
is mandatory for having the subtransient voltage at R from the value 313 kV (with 
Es) = 0) to the value 242.5kV after (3.25) and is also sufficient to carry the no- 
load current to the allowable value 390 A. For the computation of the capital costs, 
it is worth remembering that such a compensation implies the installation of shunt 
reactors in at least five sections along the line for a total rated power of roughly 
0.81 -11-120 ~ 1070 Mvar. The capability chart of Figure 3.24 (with &, = 0) is 
rather interesting for the steady state regimes, since both the power areas, which are 
sensibly far from the real axis, show power flows with strongly reactive components. 

It is also of note that receiving area (with & = 0) does not include the axes 
origin, so showing that the regime Pr +j Or = 0 (or steady state regime at no load 
at R) implies the exceeding in the cable of the set current limit, as clearly shown in 
Figure 3.24 from In, = 2170A > 1788 A= i. 

The analyses up to here have been performed for cable #b over different lengths: 
if cable #d is considered, the features become slightly worse. For instance, the com- 
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Figure 3.17 Capability charts 
without shunt compensation 
at f = 60Hz 
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Figure 3.18 Capability charts 
without shunt compensation 
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Figure 3.19 Capability charts 
with shunt reactive compen- 
sation 


Figure 3.20 Capability charts 
without shunt reactive com- 
pensation 
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Figure 3.21 Capability charts 
with shunt reactive compen- cable # b Pri=1214.0 MW |U"n=240.5 kV 
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Figure 3.22 Capability charts 
without compensation 
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Figure 3.23 Capability charts 
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Figure 3.24 Capability charts 
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parison in Figure 3.26, performed for d = 60 km, shows that cable #d suffers in the 
maximum transmissible active power a stronger reduction than in the ampacity and 
a worsening in the corresponding power factor. 

Figure 3.27 shows the capability chart with a compensation degree &, = 0.254. 


3.8.1 Theoretical Limits of the Length d 


The progressive “shrinkages” and the “shifting” in the complex plane of the sending 
and receiving areas, which can be easily detected in the cases with &, = 0, leave 
to guess that (for a given cable) there is a limit length beyond which the “basic 
constraints” cannot be satisfied: for example the Figure 3.28 shows the shrinkages 
and the locations with lengths 100, 150, 180, 190, 200 km. 

For a precise determination, Figures 3.1 and 3.2 are very useful: it can be verified 
that the points of the regimes 1) and iii) vanish if 


1 
(Lala!) «Wash 
la Ie 


for instance, the limit length of 201.87 km can be computed for cable #b, once set 
the values Up, = 230kV and J, = 1788 A. 

For this length, the sending and the receiving areas degenerate in the points S7, 
and Ry respectively which highlight a steady state regime interesting only from 
a speculative standpoint and deserve a brief comment: from the coordinates of the 
points S;7 and Rz it is possible to infer that each port S and R introduces symmet- 
rically in the line an almost capacitive power equal to that required from the cable 
of half length at no load. 

The analysis along the cable (by means of the algorithms shown in Section 3.9.2) 
shows that at S there is an almost capacitive current equal to the ampacity 1788 A; 
at half the line there is an almost null current; again at R there is a current equal to 
1788 A. 


3.9 Steady State Regimes Within Power Areas 


As seen in the previous paragraphs, the capability charts are extremely useful means 
in order to quantify and compare accurately the cable possible performance for dif- 
ferent types and lengths. 

It is anyway important to describe some additional procedures (see Sections 3.9.1 
and 3.9.2) that allow visualizing in detail the possible regimes within the receiving 
area. 
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Figure 3.25 Capability charts 
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Figure 3.27 Capability charts 5 z 
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Figure 3.28 Worsening of 

the capability charts of cable 1500 
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3.9.1 Enhanced Capability Charts 


It could be useful for transmission system operators to enhance the capability charts 
(see for example Figure 3.29 for the receiving-end area in the case of cable #b with 


3.9 Steady State Regimes Within Power Areas 73 


1500 


-500 pe 


-1500 


-1500 -1000 -500 0 MW _ 500 1000 1500 


Figure 3.29 Enhanced capability charts at the receiving-end for cable #b; d = 60 km; &, = 0 


a length of 60km and with & = 0) adding other curves referring to some current 
values J, smaller than cable ampacity, so that for each transmissible complex power 
chosen within the area the margins are immediately evident. 

For the computation of these curves the methods of the first and second analysis 
can be suitably used by fixing on the real axis current phasors (Ip Z0) and (1§ 20) 
with magnitudes smaller than the limit J. 

It is sufficient to specify that in the first analysis 6 is not constrained when 


5 C 
<F| + | Uos| < te (3.27) 
and in the second analysis ? when 
|AIS| + |C Ugs| < Ic. (3.28) 


Another procedure can be based on (3.15) inferable by (2.34) which has been 
premised in Section 3.4 and here written again: 


Equation 3.15 (valid for a complete analysis) can be interpreted by the construction 
of the two triangles (ON’N and ON’N) of Figure 3.30, which shows how, once 
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Figure 3.30 Triangles for 
Briggs formulae application 
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Figure 3.31 Enhanced capability charts at the receiving-end for cable #b; d = 60km; & = 0.608 


fixed for the magnitudes in S and R a pertaining couple of values 1%, I, lower 
than /, (and holding |Up;| = 230kV) it is possible (respecting the geometrical 
compatibility of the triangle sides and by using the Briggs’ formulae) to individuate 
the angles ~w — t anda + t and hence to know the two possible regimes (the former 
composed of the phasors J and J p,, the latter of J§ and J p,): once obtained, with 
a twofold application of (2.38), the phasors Upp, and Upp, in the two different 
regimes, it is simple to calculate the two corresponding complex powers at R. 
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The capability chart of Figure 3.29 shows at a glance the two steady state regimes 
at no load (as seen at port S, since |J.,5| & 1 kA and Jp = 0 or alternatively at 
the port R, since |p| & 1kA and J, = 0). In order to have no-load currents 
< 400A, it is necessary to adopt a shunt compensation &, = 0.608 as already 
shown in Figure 3.21; in this case, the enhanced capability chart is represented by 
Figure 3.31. 


3.9.2 Application of Ossanna’s Method 


Starting from the single-phase model of a three-phase network seen at port R, the 
Ossanna’s method [13] gives a simple (but very powerful and elegant) algorithm 
to determinate at the port itself the possible voltage regimes (and consequently the 
current ones as well) due to a fixed complex power (for each phase) Sop = Sore!?: 
the analytical direct computation avoids the iterative procedures up till now largely 
used. It constitutes a powerful tool to further investigate the steady state regimes 
inside and outside the receiving areas of the capability charts and is able furthermore 
to determine the limit conditions of physical feasibility. 

Firstly, it appears clear that the cable line system can be studied by means of the 
typical scheme as seen by the port R of Figure 3.32, by considering that in such case 


Uos 
Yorn. = (3.29) 
Uorn. _ Uos/A_ B 
a I Re Us /B A ( 


being /Rsc = short circuit current at R so that (3.31) is valid 


U -U S * 
=ORNL ~OR = de _ (= ) ; (3.31) 
Z Uor 
Jy 
Ze R Ir 
U sm us| Sore! ?. single-phase complex power 


Uoo=|Uornt] 


u = Uor/ Uop 


Figure 3.32 Network model and phasors u (per unit) 


76 3 Operating Capability of Long AC EHV Power Cables 


If the “per unit” method is adopted (exactly as J. Ossanna made at the beginning of 

the past century), and by assuming as basis quantities 
2 

Up 


=—, Z,=|Z|, 
|Z| 


Uw = |Uornil» Soo 


and by posing Uypny;, on the real axis, it is possible to obtain the following p.u. 
relations: 


—F sig = > wW-wWastd® & uli -W asd 
av u* 
=> w—ue?+siX =0, (3.32) 
where 
U U 
ORNL — yo = 1; u=—2- y=|u); u=ue? 
Uop Uon 
(see phasor diagram) ; 
Z ; S i 
=="; s= SR: sasel?: p-p=x. 
Zp - Sob = 


By explicating in (3.32) the real and the imaginary parts, a system of equations is 
reached 


ua +us,—ux+scosy=0, (3.33) 
—uy +ssnxy=0, (3.34) 


so that (3.34) gives immediately 
Uy =S-+SINY = Noss (Noss = Ossanna’s parameter) (3.35) 
and allows transforming (3.33) in 


ux + (s-sin x)* —ux +s-cosy =0 


=> w—uxyt [s-cosy + (s-siny)’] =0. (3.36) 


From (3.36) it yields: 


_ 1+ Vi-4lscos x + (ssin x7] _ | 


1 
+ —scos y—(esin x) 


2 2 
(3.37) 
Ux = 1/2 + Loss ; (3.38) 
1 
Soss = ye —scosy—(ssiny)? (os; = Ossanna’s parameter) . (3.39) 


It is also to highlight that the necessary location of Coss in the real part gives the 
conditions of physical feasibility (described in [13]). 
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So it is possible to determine analytically the two phasors u, and uw, (p.u.) possi- 
ble at the port R, each of them is compatible with the same fixed power s (p.u.): 


u, = [C/2) + Soss] —j Noss] . (3.40) 
uy = (1/2) — foss] — j [Noss] . (3.41) 


being that (3.42) is the constraint of physical feasibility so that it is advisable to fore- 
see a warning message for €o9,, becoming a complex number. The absolute quantities 
are given by U, = u; Upp and U, = uz Upp. 

These relations have constituted a fundamental premise for the wide and interest- 
ing investigations developed by Ossanna [13, 14] on the performance of AC over- 
head lines foreseeable in those years, but now they give this book a meaningful 
contribution in order to enhance the analysis of capability charts. 

In our case, it can be ascertained that the phasors u,, by (3.41), present always 
very small magnitudes (with consequently highest currents) and give (beyond the 
obvious theoretical interest) a precise reference chiefly for the possible investiga- 
tions of voltage collapse: therefore, in the capability chart study (devoted to allow- 
able regimes) only (3.40) must be considered. 

From (3.40) and (3.41) it is immediate to derive [13] the elegant relation u,; + 
(u,)* = 1 = uy + (u,)* which clearly highlights how if uw, is close to | the same 
cannot be for up. 

Therefore it becomes easy to implement few and simple files on the basis of 
(3.40) in order to automatically illustrate in detail the steady state regimes which 
exist not only at R and S ports, but also (by implementing the obvious matrix tech- 
niques described in Section 2.5) in other sections along the line once any complex 
power Sp = 3S pp is set: in this way an exhaustive investigation can be performed 
not only inside but also outside the receiving area, so that it will be also possible 
to individuate the physical feasibility independently from the constraints (3.1) and 
(3.2). 

For example (in case of cable #b, d = 60km, & = 0.608), as a completion 
of the analysis along the cable, with Ossanna’s method, it is possible to obtain, by 


Table 3.3 Three-phase receiving power S p = 3SopR 


#Regime Three-phase complex power S p 
MW + j Mvar 


1000 + j0 
0+j0 
—1000 + j0 
500 + j 500 
0 + 5500 
—500 + j 500 
750 — j 250 
0 —j 250 
—750 — j 250 


OANDNHWNK 
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Figure 3.33 Voltage and current magnitudes along the cable #b, d = 60 km (&, = 0.608), in the 
regimes 1-9 of Table 3.3 


hypothesizing the three-phase power in Table 3.3, the curves (along the cable) of 
Figure 3.33. 

It is trivial to observe that, by using the suitable references, the fixed power can 
be located in any quadrant. 

The possibility of knowing the magnitudes of the currents at different locations of 
the cable line becomes extremely useful in order to refine some temperature checks 
in those parts of the route having high thermal resistivity. 


3.10 Cables with Gas Insulation (GILs) 


The proposed procedure has been suitably applied to another kind of cable with 
gas mixture insulation namely Gas Insulated Transmission Line: Tables 3.4 and 3.5 
report the constructive data of a typical 380kV GIL and the positive sequence pa- 
rameters, respectively. 

As already mentioned in Chapter 2, the enclosures of gas insulated lines have 
a very low per unit length longitudinal resistance (see ren in Table 3.5) due to 
a highest cross-sectional area: subsequently they are always solid-bonded so that 
an induced current almost equal to that of phase (but with opposite angle) can flow 
in it giving a strongest screening effect with very low additional power losses. 

In this way there is an almost perfect decoupling between the phases so that 
a very good approximation consists in the study of GIL as a symmetric three-phase 
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Table 3.4 Typical data of 380 kV (50 Hz) GIL #5341 Al 


79 


Gas insulated line 


Cross-sectional area of phase (pipe) (Al IACS = 61%) 
Cross-sectional area of enclosure (pipe) (Al alloy IACS = 52.57%) 


Phase outer diameter 
Phase thickness 


Enclosure inner diameter 


Enclosure thickness 
Insulating gas SF¢/N> 
Pressure 


% 
bar 


5341 
16022 
180 
10 
500 
10 
20/80 
a 


* IACS = International Annealed Copper Standard (IACS = 100% => Copper Standard 


Conductivity = 58.108 m/(Q mm7)) 


Table 3.5 Positive sequence (50 Hz) parameters of 380 kV GIL #5341 Al 


Installation 


DC resistance 

of phase at 60°C 
DC resistance 

of enclosure at 50°C 
Apparent resistance 


at positive sequence (50 Hz) 
Phase-enclosure inductance 
Phase-enclosure conductance 
Phase-enclosure capacitance 
Longitudinal impedance (50 Hz) 


50 Hz shunt admittance 


Characteristic impedance 


at 50 Hz 
Propagation constant 


Natural power at 400 kV 


Capacitive current related 


to Uy = 400kV//3 
Steady state ampacity 


Gas insulated line 


Directly buried flat with axial spacing 1.3m 


Toh 
Ten 
r=Ppnt+Ten 
£ 


c 
Z=Pfpnt+len + jot 
y=stjoc 


mQ/km 
mQ/km 
mQ/km 
mH/km 
nS/km 
uF/km 


mQ/km 
mS/km 


Q 


1/km 


A/km 


A 


6.286 
2.330 
8.616 


0.204 
negligible 
0.0545 
8.6 + j 64 
0+ 50.017 


61.46 Z—0.07 rad 


0.0001 + j 0.001 


MVA Zrad 2603 Z—0.07 rad 


3.95 


2400 


line with decoupled phase conductors having a very low apparent resistance equal 
to r =Prph + Fen (about the DC values). 
Moreover, from Table 3.5 it is of note the low per unit length capacitance c = 
54.5 nF/km and the high steady state ampacity (2400 A). 
The capability charts of Figure 3.34 highlight immediately how a 100 km length 
GIL can have satisfying power flow performance with low power losses, good volt- 
age levels at R, a satisfying behaviour for line energization (see U;'p), the low Ini 
no-load current in the line de-energization, even without shunt compensation. 
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Figure 3.34 Capability chart for GIL; d = 100km 


Moreover with the procedure in Section 3.9.2 it is possible to verify that, along 
the entire GIL length, the different regimes appear wholly satisfying. 


3.11 Regimes with Ups 4 230kV 


In the previous cases it has been always set |Ups| = U, = 230kV according 
to (3.3). Other capability charts with Ups < or > 230kV could be helpful to trans- 
mission system operators in order to optimize the network regimes in daily load 
variations. 


3.12 ‘“‘Receiving Area” and ‘Sending Area’’ as Set Intersection 


The curvilinear lattices which embody the enhanced capability charts can be drawn 
by means of another type of procedure, which presents unquestionably the advan- 
tage of a noteworthy formal simplicity but has the drawback of not requiring a care- 
ful critical study on the possible interactions between the involved quantities of the 
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network (it is usually very helpful to deeply understand the knowledge of the phe- 
nomena). 

The abovementioned analyses which individuate and discriminate the steady 
state regimes compatible with the basic constraints can be also developed by means 
of the following attractive procedure (it is less circumstantial but equally valid) 
which is based on the set theory. 


3.12.1 The Determination of the Receiving Area as Set Intersection 


By utilizing (3.6)—-(3.8) written in Section 3.3 for the first analysis and by not ex- 
cluding the cases where the phasor Uy, determines the regimes with magnitudes J's 
not allowable, it is immediate that the relation 


U B.-I 
Sp = 3: (= - ait) ‘TR (3.43) 


(with 6 = 0-27 and Ip = I, Z0) represents the enclosed contour which delimits 
the set of the receiving-end complex power compatible with the two constraints (3.1) 
and (3.3). 

Then by using (3.12)—(3.15), written in Section 3.4 for the second analysis (and 
by not excluding the cases where the magnitudes of Jr become not allowable), the 
following relation 


Sin = 3A Ug = Bal) CC Ue PA (3.44) 


(with & = 0-2” and J, = J/,Z0) represents the enclosed contour which de- 
limits the set of the receiving-end complex powers compatible with the two con- 
straints (3.2) and (3.3). 
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Figure 3.35 Receiving power 
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It can be easily understood how the contour delimiting the “receiving-end power 
area” of the regimes compatible simultaneously with the three “basic constraints” 
(3.1)-(3.3) can be obtained by means of an intersection procedure of the two sets 
above delimited as the example of Figure 3.35 clearly shows (in whole accordance 
with the capability charts of Figure 3.24). 

It is important to highlight that the procedure allows identifying in automatic way 
also the contours of subsets corresponding to Jr and Js inferior to J, (e.g. for the 
drawing of the enhanced chart). 

The structure of the given relations shows that the contours Sp, are circumfer- 
ences and the contours S p, are ellipses. 


3.12.2 The Determination of the Sending Area as Set Intersection 


By following some procedures completely analogous to those of the preceding para- 
graph, it is possible to achieve the following relation: 


-—C-U Ip\* 
SS BU vg | ae as ER 3.45 
951 = 3¥os ( A 25 A (3.45) 
with 6 = 0-27 and Jp = J, Z0 and the trivial one: 
Ss = 3Uosls (3.46) 


with } = 0-2n andJ, = [,Z0. 
The intersection of the two sets gives the “sending power area”. 


3.13 The Analysis Along the Cable with Lumped Shunt 
Compensation 


Once the investigations by means of the capability charts have been made (based 
on the hypothesis of distributed compensation), it is advisable to perform a further 
analysis by realistically considering a lumped shunt compensation (according with 
the criteria in Section 2.8.2) based on the same compensation degree &, already 
evaluated with the idealized hypothesis; also in this new configuration it is possible 
to check that 


i) the energization subtransient voltage Ujp and the steady state current Jy, are 
practically within the limits already fixed; 

ii) the receiving powers compatible with the ampacity constraints and the voltages 
along the whole link from the port S to the port R have slight variations. 


Let us consider the example in Figure 3.21 (cable #b, d = 60 km, &, = 0.608 here 
repeated in Figure 3.36) where a lumped compensation is performed with shunt 
reactors installed every 20 km after the general rule recommended in Section 2.8.2, 
with the same compensation degree &, = 0.608. 
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cable #b Pei=1214.0MW 1U",=240.5kV 
d=60 km AP = 7.99MW In, = 400 A 
AQ=-172.9 Mvar | Oy, = -276 Mvar 


X"=202 Ey=0.608 


MW+; Mvar 


1214+) 0 
1000 + j 625 
| #3 | -1214-j0 


-1044 - | 474 

[#5 [ 300+) 750_| 
[#6 |  O+jo | 
| #8 | 1000+ 250_| 


1000 =: 1500 


-1000 F--- 


-1500 ——innd —4 
-1500 -1000 -500 0MW 500 


Figure 3.36 Capability charts (cable #b; 60 km; &, = 0.608) 


For the checks ii) it will be necessary (in correspondence to different meaningful 
complex powers, located along the border of the receiving area or inside it, as those 
indicated by the points 1-8 in the capability charts of Figure 3.36) to compute the 
current and voltage phasors, (or the hybrid vector w,,_¢_;9) in different ports along 
the lumped-compensated cable from R to S, and to visualize the diagrams of the 
magnitudes in Figures 3.37 and 3.38. 

It can be realized by means of the matricial expression: 


w, = Myr wre. (3.47) 


x 


where each matrix M_, p (transmission matrix between port x and port R), is given 
by the ordered product of the matrices of the system elements between x and R (see 
Figure 3.39). 

For example (see Figure 3.39) it is: 


Miyr=M,-M,:M,:Myy-M,-M,-M,-M,-M,-My., (3.48) 


where: 


M,, is the transmission matrix of the cable circuit (uncompensated) of length 
equal to one-fifth the distance between a reactor and the following one; 
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Figure 3.37 Steady state current magnitudes along the cable (cable #b; 60km; &, = 0.608, 
Uos = 230kV) with the lumped shunt compensation of Figure 3.39 in the regimes #1-8 of 
Figure 3.36 


My and Myy are the shunt compensation matrices equal to those M, and M;, 
already presented in Section 2.8.2. 


It is worth remembering that the hybrid vector wp is computable directly by 
means of Ossanna’s method (already presented and used in Section 3.9.2) once fixed 
each receiving power from #1 to #8, the voltage level Uy, and the matrix M gp. 

The checks i) of page 82, regarding the subtransient regimes of energization at 
no-load and the steady state current J, , can be performed by extracting the param- 
eters A and B from the matrix M gp and by applying successively (3.21)-(3.24): in 
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Figure 3.38 Steady state voltage magnitudes along the cable (cable #b; 60km; &, = 0.608, 
Ups = 230kV) with the lumped shunt compensation of Figure 3.39 in the regimes #1-8 of 
Figure 3.36 
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Figure 3.39 Composition of matrices for the analysis of lumped shunt compensation (cable #b; 
60 km; &, = 0.608) 
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this way it yields (with & = 0.608) the values Uj‘, = 240.55 kV (with X” = 20) 
and Jy, = 401 A (almost equal to the results given by distributed shunt compensa- 
tion). 

By means of other examples regarding steady state regimes #1-8, it can be as- 
certained that generally the hypothesis of ideal distributed shunt compensation is 
compatible with the real installation of lumped compensation: some differences in 
the receiving capability (e.g. with Sp #2 or Sp #4) gives some slight operation 
limitations. 

The described procedure is obviously applicable to any configuration also very 
irregular for distances between reactors with different rated power. 


3.14 Conclusions 


This chapter gives efficient means in order to evaluate the power capability of EHV 
AC cables. The power capability charts represent an immediate and precise outlook 
on the possibilities of active and reactive power flows along the link. In the same 
charts, the regimes that give rise to receiving-end voltages too high or too low are 
highlighted. The method is based on the classic transmission formulae and has been 
fully developed in matrix form. 

The procedure for the computation of the voltage and current magnitudes in any 
point along the cable has also been described. In many cases, the cable limit lengths 
can be overcome by means of shunt compensation: the capability charts show the 
different benefits depending upon the compensation degree. 

The procedure has been also successfully applied to gas insulated transmission 
lines showing promising performance in power transmission. 


References 


1. R. Benato, A. Paolucci: Operating Capability of Long AC EHV Transmission Cables, Electric 
Power Systems Research, Vol. 75/1, July 2005, pp. 17-27 

2. R. Arrighi in the name of WG 21.13: Operating Characteristics of long links of ac high voltage 
insulated cables, Cigré 1986, Paper 21-13, pp. 1-10 

3. IEC 60287: Electric cables — Calculation of the current rating, (in 8 parts: 1.1, 1.2, 1.3, 2.1, 
2:23.13 3:2;.3.3) 

4. G.J. Anders: Rating of electric power cables: Ampacity computations for transmission, distri- 
bution, and Industrial Applications, IEEE Press, New York, 1997 

5. G.J. Anders: Rating of electric power cables in unfavourable thermal environment, IEEE Press, 
Wiley-Interscience, New Jersey, 2005 

6. S.Y. King, N.A. Halfter: Underground power cables, Longman, London and New York, 1982 

7. E. Peschke, E. von Olshausen: Cable systems for high and extra-high voltage, 1999, Publicis 
MDC Verlag, Germany 

8. IEC 62271-100 High voltage alternating-current circuit-breakers, 2001 


References 87 


9. British Electricity International: Modern Power Station Practice, EHV Transmission, Vol. K, 
Chapter 9, Sect. 3.3: Control of switching surges, Pergamon Press (Third Edition, 1991, Ox- 
ford) 

10. W. Diesendorf: Insulation Coordination in High Voltage Electric Power Systems, Butter- 
worths, London, 1974 

11. C. Concordia: Synchronous Machines, J. Wiley & Sons, New York, 1951 

12. E.W. Kimbark: Power system stability, J. Wiley & Sons, New York, 1948 

13. J. Ossanna: Neue Arbeitsdiagramme tiber die Spannungsdnderung in Wechselstromnetzen, 
EuM, 1926, pp. 113-126 

14. J. Ossanna: Ferniibertragungsmoglichkeiten grofer Energiemengen, ETZ, 1922, pp. 1025- 
1029, 1061-1063 


Chapter 4 


Operating Capability of AC EHV Mixed Lines 
with Overhead and Cables Links 


4.1 Introduction 


The previous chapter has been devoted to the operating capability and constraints 
involved in power transmission of long AC cable links. 

In this chapter the approach has been extended in order to show the transmission 
performance of some typical mixed configurations with overhead lines (OHL) and 
underground cable lines (UGC). 

In order to erect a new long link, (but also, often, to fit an existing one due to 
new territory necessities) it could be necessary to have a mix of transmission tech- 
nologies such as a cascade connection of OHL-UGC-—OHL. In fact, the overhead 
line can be unsuitable in certain parts of the route for tackling some environmental 
and orographical problems, which can be easily overcome with (as already stated in 
Chapter 3) UGCs. 

For example, an UGC section allows a HV/EHV link to pass inside galleries 
or through areas too wide for OHL span such as large lakes or arms of the sea: 
it is sometimes called a “siphon”. This configuration may also solve some critical 
issues, due to strong local oppositions or some laws, since it drastically reduces 
the environmental and magnetic impacts and hence it permits the UGC line to pass 
through or near a protected site (sensible places such as a school or natural park) 
or an urbanized area also with economical convenience (see Chapter 6). It is worth 
remembering that some important installations of mixed OHL—UGC links have been 
realized for several years to allow great power flows towards town centres. It must 
be emphasized that the use of UGCs gives usually the economic burden of the shunt 
compensation reactors. 

For the aforementioned reasons, the transmission system operators will have to 
face more and more the analysis of mixed lines aiming at achieving their possible 
performances without jeopardizing the system safety and the power quality. 
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4.2 Mixed Lines: OHL—UGC-OHL 


A rigorous approach to the issue requires to fix, for each OH and UG section, the 
specified constraints regarding the safety and the lifetime of section itself and to 
research successively the maximum compatible performances of the whole link. 
A simpler approach [1] can be adopted when the situation is that of Figure 4.1 
(© OHL-@ UGC-© OHL, with lengths d), d2, and d3), where the intermediate 
section is constituted of UGC: in this (very usual) case the steady state constrains 
expressed by (4.1) and (4.2), analogous to (3.1) and (3.2), appear sufficient to satisfy 
the ampacity not only of the cable but also of OHL © and ®: 


Lal<te, (4.1) 
[Peds tes (4.2) 
[Wos| = Ure (4.3) 


the value of J, represents a current level which must not exceed, with a conservative 
margin, the steady state ampacity of cable section @: it is sufficient that 7, (imposed 
at H and K) is little lower than the ampacity of OHL © and © to obtain (as it 
will have to be systematically verified) allowable current levels also at S and R. 
Ultimately, (4.3) fixes the voltage at S', similarly to the settlement already adopted 
in Chapter 3. Also in the following treatment, where U,, = 420 kV, the level Up. = 
230kV = 95% Unifa/3 has been considered. 

For each line section a suitable “positive sequence model” has been assumed: 
this requires phase transpositions for OHL and for cable line also the cross-bonding 
arrangement and allows the use of classical transmission equations based on positive 
sequence parameters (see Table 4.1 for OHL and Table 4.2 for single-circuit UGC) 
and of relevant transmission matrices (see Chapter 2). 

This is a useful guide also for little unbalanced systems provided that a slightly 
conservative choice of J, and Uo, is made. Successively, a deep insight can be per- 
formed by means of multiconductor study similarly to Chapter 5 [2-4]. 

In the case hereafter considered, for OHL a typical 380kV line of Italian grid 
(equipped with bundles of three sub-conductors) has been assumed and for UGC @ 
almost always a 380kV double-circuit consisting of single-pole 2500 mm? copper 
conductor XLPE cables (in the same configuration as in Figure 6.2 of Chapter 6): by 
considering the wide circuit spacing (5m), the two single-circuits can be assumed 
thermally and electrically independent so that, for the double-circuit, the ampacity 
is 2 x 1.788 A and the electrical model is computable after Section 2.7. 

Once considered the values of ampacity for OHL (foreseen by Italian Standard; 
see Table 4.1), the levels 7, (during cold months) = 2800 A and J, (during warm 


S @® H. © (/K  @ R 


———- > 
UV d) TP d> ie d3 


Figure 4.1 Typical sketch of a mixed line 
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months) = 2100 A have been assumed in (4.1) and (4.2), as specified in each case 
studies. These choices do not seem to be excessively redundant for section @ with 


I, = 2-1788A as shown in Table 4.2: (24288 = 1.28 during cold months and 


2800 
ine = 1.70 during warm months) because it seems a conservative installation 


criterion to have a good power reserve in case of a circuit failure. 

The analysis of the link between nodes S and R, not inserted in the real context 
of a given electric network, constitutes a preliminary assessment and a useful guide 
for further network studies (among which the power flows and network simulations 
are indispensable). 

In analogy with Chapter 3, the regime of this mixed line can be achieved by 
means of a first analysis (Up, and | x| constrained) and successively a second one 
(Ups and |J ;,| constrained), with the aim at individuating the regimes that respect 
the constraints (4.1)—(4.3). 

With regard to the steady state capacitive power Qc absorbed by the UGC con- 
sidered in this chapter (a double circuit of 2500 mm7? cables) it can be roughly es- 
timated (refer to Section 3.2) Oc = 22 Mvar/km at 50 Hz and = 26 Mvar/km at 
60 Hz: it requires almost always a suitable degree of shunt compensation &,, both to 
improve the steady state and transient regimes (see Section 4.7). 


Table 4.1 Positive sequence parameters (50 Hz) of OHL with bundles of 3 ACSR sub-conductors, 
(® 31.5 mm, spacing 400 mm) 


Overhead section © and ® 
Voltage level 380kV 


Resistance at 75 °C (50 Hz) r mQ/km_ 23.10 
Series inductance £ mH/km_ 0.858 
Shunt leakance (50 Hz) g  nS/km 10 
Shunt capacitance c  uF/km 0.0133 
Ampacity (cold months) — I, A 2955 


Italian standard 
Ampacity (warm months) — 2220 
Italian standard 


Table 4.2 (Cp. Table 3.2) Positive sequence parameters (50 Hz) of single-circuit UGC with single- 
core cables directly buried in cross bonding at a depth 1.45 m in the soil with py, = 1.0K-m/W; 
spacing s = 0.35 m; flat laying and soil ambient temperature 3, = 20°C 


Cable section @ 


Voltage level 380 kV 


Phase cross-section mm? 2500 Cu 
Phase resistance at 90°C r mQ/km 13.3 
Series inductance £ mH/km 0.576 
Shunt leakance with tanéd = 0.0007 g = nS/km 52 
Shunt capacitance with €, = 2.3 c  uF/km 0.234 
Steady state ampacity I, A 1788 
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Nevertheless, the value of &, does not modify the methodology for the first and 
second analysis in order to construct the capability charts (presented in Sections 4.4— 
4.6) so that it is more convenient to calculate it definitively in Section 4.7. 

In any case, the hypothesis of uniformly distributed shunt compensation can be 
adopted also in this chapter (analogously to Chapter 3): it allows constructing easily 
the corresponding transmission matrix to introduce in the system analysis. 


4.3 The Transmission Matrices for the System Study 


For each section ©, @ and © and for their cascade composition, the corresponding 
classical matrices can be derived as in Section 2.6. 

It is worth underlining that in the following case studies the matrices M, and 
M_., (OHL) must be computed by considering the corresponding lengths d; and 
d3; the matrix M., must be computed by considering both the length d and the 
compensation degree &, determined as in Section 4.7. 


4.4 First Analysis 


Once fixed at node K the current phasor J x = I, Z0 (on the real axis), and imposed 
at node S the following phasor 


Uys = 230kVZ8 (8 = 0-27), 


a set of regimes of the entire line is univocally achieved, among which only those 
respecting the limit J, at node H must be individuated. In particular (4.4)-(4.11), 
show the subsequent formal passages, which give (4.13), deriving from (4.11). The 
expression (4.13) can be interpreted by the phasorial diagram of Figure 4.2 (with 
some modified proportions for graphical purposes) which, by means of the relations 
shown in the same figure, yields the problem solution: for 6 ranging between 6, 
and 55 (evaluated by means of (4.14) and (4.15)), the voltage phasor Uy, (4) deter- 
mines in this first analysis all the current phasors J 7, (e.g. OL, OL3, OL4, OL2) 
having magnitudes within the limit /, according to (4.2) as long as the magnitudes 
of phasors a, b, c are compatible with the triangular configuration. In particular, the 
phasors OL, and OL, (the former leading and the latter lagging with respect to J x) 
regard currents J ;, with the same magnitude of J, = |J x|. 


Uos| | Ar | Bir | ox 
= ’ (4.4) 
Is Cy, | Dy ||Lr 


Us s(K) My, UK 


4.4 First Analysis 


Figure 4.2 Phasorial diagram of (4.13) 
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O0'=(Aj/A)2) Ik a 

O'H=[( A) C)2/A12)-C)]U os 20 
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rr ee 
cos A = “al z (Briggs) 
(4.5) 
(4.6) 
(4.7) 
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imposed phasors, (4.8) 
1 YsscKy) > YHAK) = Mi: N51 25x. (4.9) 
(4.10) 
_ si 
Aj) 
' (4.11) 
A; 
=-712 
N 1 
(4.12) 
A, i Cy A, 
Ap Yos + 4, —* (4.13) 
6; =a-—fB,—n, +1 (where pp; <0), (4.14) 
69 = 6) + 2B : (4.15) 


4.5 Second Analysis 


In an analogue way, once fixed at node H the current phasor J, = J,Z0 (on the 
real axis), and imposed at node S the following phasors 


Ups = 230kVZ9 (9 = 0-2nr) , 


another set of regimes can be achieved among which only those respecting the limit 
I, at node K are of interest. 

The relations (4.16)-(4.23) show the subsequent formal passages, which give 
rise to (4.25). It can be interpreted by the phasorial diagram of Figure 4.3: by means 
of the relations shown in the same figure, it is possible to see that, for 7 ranging 
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Ly=1.20; 

OO'=(A12/A)) In> a 

O'K=[( 412C1/A))-C12]U os 29 
O'L"=[( A2C)/A))-C)2] Us 29"> C 
OL"=00't[( AnGi/A1)-C12|U,s 28" 
OL'= OO'+[( A72C/A1)-C)2] Us 29" 


|OL"|=|OL'=1. 
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41 


A 
= age in| 
4, 


Figure 4.3 Phasorial diagram of (4.25) 


between 3” and 1 (evaluated by means of (4.26) and (4.27)), the voltage phasor 
Uys (8) determines in this second analysis all the current phasors I x (e.g. OL", 
OL”, OL", OL’) with magnitudes within the limit J,, according to (4.1). 


Uos A, | B, | Gow 
i;| |c,|D,|\ aa (4.16) 
Uss(H) M, vy 


ls (4.17) 


ie (4.18) 
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1 
1 0 
Noy = . (4.19) 
—C,| A 


Dos 1 0 | |\Uos 


Cy} a (4.20) 
fs| | a. | a | [4a 
Uss(H) N50 UsH 
UeKH) =Mn-ssay > UKKGD =Mi-Ns2-¥sH - (4.21) 
UK| Diy —By 1 0 Uos 
— C s 
Ix] Cx) An] a | a | [La ve 
URK(H) My Ns. UsH 
B 
Dox) |Pn- ByG “ae Dos 
a ' (4.23) 
Ix Cyt Ary a aa Ly 
UKK(A) Nx YUsH 
Nxo=MyNs5. (4.24) 
A,,-C A 
Lae = (-Cio+ SS) ta + Ede, (4.25) 
A, A, 


In particular, the phasors OL” and OL’ reproduce again the same relations of mag- 
nitude and angle (with respect to |/ ;,|) detectable in the first analysis for OL; and 
OL. 


0” =n —Bo—12+ 2 (where pt. > 0), (4.26) 
! = 9" +2. (4.27) 


4.6 The Capability Charts 


The procedures presented in the first and second analysis allow individuating, for 
a given mixed line, all the regimes at power frequency compatible with the con- 
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straints (4.1)-(4.3), yielding for each regime the vector (composed of voltage and 
current phasors) at each node: it becomes particularly interesting to evaluate the 
complex (three-phase) power S.; and Sp by using respectively the elements of 
vectors Us s(x) in (4.8) and vs 5(z7) in (4.20) and the elements of vectors UrR(xK) 
in (4.28) and v ppcqy in (4.29). 

For the explicit form of Np, and N p, refer to Section 4.9. 

It is worth noting that (4.28) and (4.29) as well as (4.8) and (4.20) are useful 
not only to compute the end port power but also to verify (as already premised in 
Section 4.2) the respect of ampacity at R and S for OHLs. 


Uor 

Te} EC dys || £2 | | 

AR ©1343) 4. | An | ( (4.28) 
URR(K) M1); Ns, = Usk 

Nri 
Uor 
~ C 1 

IR HC 193) Ags A A, ly (4.29) 
URR(H) M)>, N 52 UsH 


The respect of ampacity along the UGC length appears (in any case-study) as- 
sured a priori, since the ratios between ampacity and J, are 1.28 (cold months) and 
1.7 (warm months). If the ratio were less conservative, it would be advisable to make 
some checks by means of Section 4.10. 

Hence the outlines can be drawn (see the example of Figure 4.4): the outline s 
bounds the “sending-end power area” and the outline r the “receiving-end power 
area”, respecting the fixed constraints. It can be ascertained (see Section 4.9) that 
they are composed of arcs of circumference and ellipse, analytically inferable. 

At regular steps, it is possible to mark the complex power S$, and S p with couple 
of points upon the outlines s and r: their correspondence is easily detectable by the 
mutual proximity (see Figure 4.4). The difference S$. — Sp views with the real 
part the active transmission power losses and with the imaginary part the reactive 
power demand by the mixed transmission system (including reactive compensation 
of UGC, if any). 

The couples of points marked with x represent the particular regimes character- 
ized by |I.x| = |Ly| = Le. 

These capability charts for mixed link give (as already seen for only cable lines) 
immediate visual information on power flow performances. 

As shown in the first example of Figure 4.4, each capability chart is completed 
with the following transmission parameters analogously to Chapter 3: 
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Figure 4.4 A first example of 
capability charts 
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Pr=1850 MW U"GR=242.5 kV 
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d,,do,d3 lengths of sections ©, @ and @; 

I. chosen current constraint at H and K; 

x" 50 Hz subtransient reactance (seen at port S); 

Esh shunt compensation degree (for section @ according to Section 4.7); 

OR RMS value of the subtransient voltage in the no-load mixed line ener- 

gization from the port S for the given degree &,n; 

Ixy RMS value of the current absorbed in S by the mixed line at no load 
in R, in steady state regime with |Up;| = 230kV and with the given 
degree &y; 

ONL Capacitive reactive power absorbed in S by the mixed link at no load in 
R in steady state regime with |U), | = 230kV and with the given &,,; 

Pri maximum active receiving power with cosg = 1; 

AP, active power losses in the mixed line at PR;; 

AQ, whole reactive power required by the mixed line at Pp}. 

The maximum active receiving-end power Pr, (at cos = 1) is generally de- 


tectable in the around of 6, in the first analysis or in the around of # in the second 
analysis and appears very meaningful in order to evaluate the more interesting per- 
formances and power flows of the mixed line, together with the active losses A P| 
and the reactive power AQ, required by the mixed line. 
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4.6.1 Phase Voltage Levels at R 


The voltage levels at node R play a key role in the network operation and can be 
directly visualized in the “receiving-end power area” analogously to Section 3.6: in 
this case the following formula must be applied 


Sr = 3UorlUos — AirxUor)/Biys)” - (4.30) 


4.7 No-Load Energization and De-Energization 


In the mixed lines the energization and de-energization phenomena deserve careful 
consideration also in the cases where the cable length d> is small, since the induc- 
tive reactance of the section OHL, which adds to X”, tends to heighten, in no load 
regimes, the whole capacitive effect. 

The sketch of Figure 4.5 shows the switch-on of the circuit breaker in S in order 
to energize the mixed line. 

As already seen in Section 3.7, the supply in S can be modelled as an equivalent 
generator which is characterized by its electromotive force Up (supposed to be equal 
to 230kV) and the short-circuit subtransient impedance Z” (for simplicity purely 
inductive Z” =jX”). 

The power frequency subtransient voltage Uj}p, due to the closing at S and with 
R open, is completely defined by the following formulae: 


" Uo A, 23 


Gos = ; , (4.31) 
Z" + fu Cos 
” 
sp = ee, (4.32) 
a Alp; 
Uo 


i —_____,_—— , 4.33 
= ae Ajy3 + Z" + C493 eee 


being A,,3/C ,3 the impedance (almost completely capacitive) as seen from S' with 
R open, and expressing with (4.32) the Ferranti’s effect; for the subtransient reac- 


tance X” values ranging in the interval 4-23 Q can be assumed as already seen in 
Section 3.7. 


Z 5 O ug @ ©) R % 2 
di d) d; 
Mp3 ac 


Figure 4.5 Typical sketch for energization of mixed line at no-load at R (subtransient regime) 
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In the hypothesis that, extinguished the transient switching phenomena, the net- 
work voltage regulation restores again at port S the rated value Up = 230kV, the 
no-load steady state (almost entirely capacitive) current which the circuit breaker 
must be able to interrupt is computed by 


ee 
ANL = : 
(Aj93/C€ 493) 
With regard to the energization phenomena, in order to respect the standard switch- 


ing levels (e.g. 1O50kV) with a conservative margin, it is advisable (as in Sec- 
tion 3.7) to impose the following constraint to (4.33) 


(4.34) 


|UGe| < 242.5kV = Un/V3. (4.35) 


As regards the de-energization phenomena, the following constraint to (4.34) must 
be imposed: 


[Zn] < 400 A (4.36) 
which expresses the suggested limit of the already mentioned standard. 
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Therefore, once fixed the parameters d,, d2, d3, X"”, Uo, the values of &, which 
give the respect of both the abovementioned constraints can easily be computed with 
automatic iterations of (4.33) and (4.34). 

In the numerous examples presented in the following Section 4.8, it will be pos- 
sible to distinguish which constraint (4.35) or (4.36) is decisive for the choice of the 
shunt compensation degree, depending on the different configurations. 

Since the determination of &, depends upon numerous parameters, the Fig- 
ures 4.6—4.9 report only some possible situations that are very useful to frame the 
phenomena. 

For example, the curves of Figure 4.6 (where &, = 0) show that the limit length 
dy of UGC are always very close to 10 km and that the constraint expressed by (4.36) 
is almost always decisive. 

The curves of Figure 4.7 (computed for X” = 15 Q) show that (set d3 = 0) 
the section UGC with &, = 0.7 can reach d, = 40km if d,; = 0 (being that the 
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constraint (4.36) is decisive as shown in Figure 4.7b) or dz = 32 kmif d; = 60km 
(being that the constraint (4.35) is slightly decisive as shown in Figure 4.7a). 

It is worth noting that the OHL © prevailingly introduces a longitudinal reac- 
tance (= 0.27 Q/km at 50Hz) which “increases” the short-circuit reactance X” 
worsening the no-load energization subtransient regimes supplied by S. 

Hence, it can be ascertained, by applying iteratively (4.33) and (4.34), that with 
d, = 90km (by holding again &, = 0.7 and d3 = 0) the constraint (4.35) would 
lessen the limit length to d) = 26km. 

Figures 4.8 and 4.9 consider examples where d3 4 0. 

The overlapping of the three curves (grey, dotted, black) in Figures 4.8b and 4.9b 
is due to the particular lengths foreseen for d; and d3 which always give rise to 
the same whole length of OHL (equal to 80 km) and change very little the no-load 
steady state with the UGC of length dp. 

Besides, Figures 4.8a and 4.9a highlight how the subtransient no-load regimes 
are affected by the three configurations in a different way. 
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Figure 4.8 Examples of limit 
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Figure 4.9 Examples of Eon 
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By hypothesising a shunt compensation of UGC with & = 0.7, the length limit 
dy = 32km can be derived graphically (in case of dj = 20km and d3 = 60km), 
due to (Figure 4.9b) the constraint (4.36), whereas (in case of d} = 60km and 
dz; = 20km) it is d2 = 30 km due to (Figure 4.9a) the constraint (4.35). 

For further example, by setting d; = 90km and d; = 20km, the limit length 
dy = 21km, due to the constraint (4.35), can be derived numerically by apply- 
ing (4.33) and (4.34). 

It is evident that the energization of no-load mixed line from the node R must 
be evaluated by considering the lengths of three section in the right order: for ex- 
ample, if the order from S to R of the lengths d, dz, d3 is 60, 30, 20 km, when the 
energization from S is assessed the grey curve in Figure 4.9 must be considered, 
whereas the dark curve must be considered with energization from R. 

If the equivalent generators at both ends were equal (same Up = 230kV and 
same X” = 15Q) it would be better the energization from R, which requires 
a lower compensation degree. 
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On the other hand, it appears very likely that the receiving-end is a very weak 
node (very high X”) so that it results less suitable than S for the energization of 
no-load line: also in these cases, (refer to Section 3.7) the detection of “best-end” 
switching suggested by CEGB [5] is advisable. 

It is obvious that the present procedure is applicable to numerically analyse any 
situation without recurring to specific curves (e.g. change of the ports S and R, with 
magnitude Ups 4 230kV). 

It is trivial to remember the possibility for the TSO of energising the line with 
voltage temporarily reduced. 

These meaningful considerations on the energization and de-energization regimes 
must be followed by an in-depth analysis of power flows and network simulations 
in order to fix definitively the suitable compensation degree for UGC section. 

It is worth remembering, as already stated in Section 3.7, that possible cases, 
where the steady state capacitive power (On, ~ —3UpINL), absorbed by the no 
load mixed line, exceeds the ability of synchronous generators located in the neigh- 
bourhood (dangerous self excitation conditions), must be avoided. 

As already stated, the hypothesis of uniformly distributed shunt compensation 
has been adopted: it has been throughout verified that both the capability charts and 
the line energization and de-energization phenomena remain almost unchanged by 
applying (with the same &,) lumped shunt compensations (see suitable matrices 
M, in Section 2.8.2) on length intervals not greater than 15-20 km. 

Some attention must be paid (in the cases of lumped compensation) for the loca- 
tion of the limit J, at H and K (in the matrix cascade). 


4.8 The Use of Capability Charts as a Guide 


The example of Figure 4.4 already gives an idea of the utility of the capability charts 
for transmission system operators. 

It is worth remembering that the value of &, (evaluated after Section 4.7), is not 
only fundamental for a good line energization and de-energization but also deter- 
mines satisfactory power flows well evident in the same figure. 

Moreover the set of curves parameterized with phase voltages 220-240kV high- 
lights immediately how a good voltage level at R constitutes a further limiting con- 
straint (chiefly for reactive power) within the capability charts for J, = 2800 A. Also 
the values of maximum transmissible active power Pr; = 1850 MW together with 
the power losses A P; and the compensation degree &,, are very useful for economic 
evaluations: in particular, the compensation degree constitutes a direct reference to 
the installation cost of shunt reactors in section @. 

The no-load steady state current Jy, = 391A and the corresponding reactive 
power Qn, = —270 Mvar do not appear problematic. 
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Figure 4.10 Second example 
of capability charts for mixed 
line 
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A second example (see Figure 4.10) highlights how the increased overall length 
of mixed line (90, 30, 20 km) yields now unacceptable voltage drops unless the node 
R has not capacitive load. 

So it is possible to ascertain that when the mixed line length increases, the cri- 
terion of voltage drops and power losses becomes paramount: in correspondence 
of the power flow Pr; = 1767 MW it has |Upr| = 210kV (as graphically in- 


ferable) unsatisfactory for network operation, a ratio AP, /PR, = 3.65% and also 
a reactive power AQ, = 564 Mvar rather high from both technical and economical 
standpoints. 


Figure 4.11 (third example) clearly shows how choosing a current limit of 2100 A 
(just below OHL warm month ampacity) gives regimes more rational than those of 
Figure 4.10. 

Both Figures 4.12 and 4.13, which can represent some examples of mixed line 
for penetration into town centres (d3 = 0), show how the decreased overall length, 
with reference to the previous case (110km rather than 140km), allows (in cold 
months) a higher chosen current limit (2500 A) and good voltage levels at R. 

The comparison between Figures 4.14 and 4.12 underlines the reduction of ca- 
pability charts due to the current limit (2100 A) and the subsequent better voltage 
levels at R. 
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Figure 4.11 Third example 
of capability charts for mixed 
line 


Pr=1390 MW Uo R=242.5 kV 
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Figure 4.13 Fifth example 


of capability charts for mixed Frew. cate EM 


line AP |= 36 MW Tyz=360 A 
AQ = 189 Mvar |Qy;=-249 Mvar 
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Figure 4.14 Sixth example 
of capability charts for mixed Pr=l413 MW [U"R=242.5 kV 
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Figure 4.15 Example of ca- ~ Pees 
pability charts of UGC line ERE MINY, | ae aaa 
(d, = dx = 0); double- ; AP /=15.6 MW. |Iy,=400 A 

circuit with shunt compensa- < AQ |= -76 Mvar_ |Qy=-277 Mvar 
tion 


Ey=0.803 


-2500 —_t tk iii 
2500 -1500 5000-500 1500 2500 
MW 


Figures 4.15 and 4.16 also appear very meaningful because they deal with cases 
where the line between S and R is entirely constituted by double-circuit UGC of 
60 km (having set d; = 0, dz = 60, d3 = 0). 

Figure 4.15 shows the excellent performance (Pr; = 2398 MW) of 60 km com- 
pensated double-circuit cable (due to its very high ampacity limit 3500 A) with good 
voltage levels because of its low longitudinal impedance. 

It is worth noting that the corresponding technical-economical burden in this 
case is rather high since the overall shunt reactive compensation power is (at 50 Hz) 
Ox = 0.804 - (23 Mvar/km) - (60 km) = 1100 Mvar and must be shared at least in 
three sections along the line (see Section 2.8.2). 

The prohibitive conditions which would exist in the case of uncompensated UGC 
(with the same length) are shown in Figure 4.16: the highest steady state no-load 
reactive power On, = —1423 Mvar with corresponding Jy, > 2000A and the 
subtransient voltage level Uj/, = 272 kV during energization. 

It is also of note that the typical shape of capability charts with cusps at points x 
recurs chiefly each time there is a low compensated UGC as in Chapter 3. 

In these cases (chiefly with high d, and low &,y), the differences S 5 — S p appear 
graphically very evident due to the strong reactive (capacitive) power absorbed by 
the cable. 

In the mixed lines (chiefly with high d,, d3 and &,,) the reactive power absorbed 
by the overall link becomes of inductive nature (see Figures 4.4 and 4.10-4.14), 
being prevalent the inductive demand of OHL. 
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Figure 4.16 The prohibitive 7 
conditions (Uvr; In; AQ)) Pr =2037 MW U oR =272.2 kV 
of uncompensated double- AP) =12.6 MW [Jy7= 2060 A 

circuit UGC (60 km) 
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Figure 4.17 Example of ” 9340 kV 
capability charts for 100 km aa 
OHL and 1 km UGC 
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The proposed procedure has been extensively applied to many mixed line con- 
figurations but can give extremely useful suggestions for entirely overhead line as 
well: if, in the configuration of Figure 4.1, the length of UGC section is fixed as 
do < d, + ds, it is possible to highlight with good accuracy the foreseeable be- 
haviour of the considered OHL. 

A meaningful example is given by the capability charts of Figure 4.17 which 
demonstrates very restricted limits for satisfying voltage operations (owing to high 
voltage drops), high active losses (see 54.5 MW!) and a huge quantity of inductive 
power absorbed (see 554 Mvar!). This kind of operation (at the thermal limit) is of- 
ten due to the big hindrances in constructing new OHLs, in several countries, which 
compel to make maximum use of the existing overhead grid. Alternatively, the pro- 
cedures of Chapter 3 can be suitably used by introducing the “positive sequence 
model” of OHL (instead of UGC) with d = 100km. 


4.9 “Receiving Area” and “Sending Area”’ 
as Intersections of Sets 


Once imposed (as in 1’ analysis) the phasors Ix Z0 = I, and Ug (6), if (4.28) is 
made explicit it is possible, being 


Ap, = Dy3 — Bi3€ 19/Ayy > 


Nia Ari| Bri Bry = —Byp3/Ay - 
Nri= , 
Col Dig: Cry = —Cy3 + Ai €y2/Ap ; 
Dry = Ay3/An » 
to obtain (4.37) 
Srl) * 
a = [AriUos + Brilx)-(CriUos + Drilx]) 


= Api CR os’ + BriDailk 
+ Ap, Dai TRU 5 (8) + Bry CR IKU G5 (5) . (4.37) 


Equation 4.37 gives the values of ae) for 6 = 0-2z,, without respecting the limit 
I, at node H. 

Analogously once imposed (as in 2’ analysis) the phasors Jy Z0 = J, and 
Uo 5 (2), if (4.29) is made explicit it is possible, being that 


Ary = Dj3 — By3€/A, 
Bry = —By3/A, » 

Cro = —Cy93 + Aqy3€)/A, ; 
Dp = App3/A, ; 
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to obtain (4.38) 


Sr(2*) ‘ 
———— e [ArUos + Brolw): [CRUos + Dro!A]) 


= AgCR Uns + BeoDinli 
+ ApDrolHUos(3) + BarCRlHUGs(9) - (4.38) 


Equation 4.38 gives the values of Sr) for 3 = 0-27, without respecting for the 
limit J, at node K. 

The simultaneous representation of both curves Sp(1%) and Sp(2%) which 
bound each other, allows drawing (once fixed J.) the bound r of “receiving-end 
power area” (see Figures 4.18 and 4.19). 

With similar procedure, by developing (4.8) and (4.20), it is possible to have 
(4.39) and (4.40) 


S54) " % a 
a: = Ups[Cs\Uos + Dg Ix]* = CSUs? + Le DS Uos(6). (4.39) 


= Ups[C. gU 5 + Dgalu)* = C§\Vosl’ + Ly D5Uos@) . 
(4.40) 


Ss(2*) 
3 


so that, by means of S,(1%) and S,(2%), the bound s of “sending-end power 
area” can be drawn. Equations 4.39 and 4.40 give circumferences, whereas (4.37) 
and (4.38) give ellipses with low eccentricity since beyond the third term with Uys 
there is a forth term (generally slight) with the complex conjugate U5,. 

This kind of approach, which determines the sending-end power area and re- 
ceiving-end power area as intersection of sets, can become very useful to draw the 
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Figure 4.19 Delimiting the 2500 
receiving areas by means of 
(4.37) (curve 1’) and (4.38) 
(curve 2”) 
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capability charts of a line with many and different (for location and magnitude) 
ampacity constraints. 


4.10 Analysis Completion 


It can be useful and suitable to complete the analysis of the system in order to know 
and qualify in detail different regimes concerning complex power more or less far 
from the capability chart area outlines, both those located inside the capability area 
(surely compatible with the current constraints in (4.1) and (4.2)), and those located 
outside (which surely give more or less strong derogations of the same constraints); 
in some cases, it is furthermore interesting to verify how the constraints themselves 
are respected along the line. 

For these investigations, it is extremely useful to take advantage of (analogously 
to Chapter 3) Ossanna’s method and the matrix algorithms. 


4.10.1 Analysis Completion by Means of Ossanna’s Method 
and Matrix Algorithms 


Before proceeding with the use of Ossanna’s formulae reported in Section 3.9.2 
it is necessary and sufficient to explicitly state that in these cases of mixed lines 
A = Aj; and B = B,,, must be always set. 

In this way, once fixed, any complex power inside the receiving-end power area, 
one immediately obtains the corresponding voltage and current phasors at R, which 
allows determining not only the corresponding regime at sending-end port (so that 
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Figure 4.20 Current and voltage magnitudes along the mixed line (see Figure 4.10; d; = 90km; 
dy = 30km; d3 = 20km; &, = 0.79) 


the transmission angle and power losses can be immediately computed), but also 
(by means of suitable matrix algorithms detailed in the Chapter 2) the regimes at 
different locations along the line. 

It is interesting to highlight how the method can be used also to evaluate and lo- 
cate the stresses in the systems during possible low differences in the power outside 
the capability outlines. 

Figure 4.20 shows an immediate application of Ossanna’s method: the current 
and voltage magnitudes along the mixed line of Figure 4.10 are clearly detectable 
for the four shown complex power required at receiving-end. It is worth noting that 
Sp = 1500 MW + j500 Mvar lies slightly outside the receiving capability area of 
Figure 4.10: in fact, the current magnitudes of Figure 4.20 exceed, slightly along the 
line, the ampacity of 2800 A up to 2900 A at port R (with too high voltage drops). 


4.11 Circuital Considerations 


It is worth premising that the considerations which will be developed in this para- 
graph give a useful contribution for the direct interpretation of the matrix relations 
obtained up to here, allow giving the single elements of the matrices themselves 
simpler expressions (even if with exact equivalence): the circuital analysis could 
represent a useful exercise for a student of electrical engineering. 
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4.11.1 The Three Matrices N 44, Ns Np; 


In the following, the matrix relations which summarize in symbolic form the above 
presented equations (4.11), (4.8) and (4.28) are reported: 


(4.11) 


Vos 

Ix (4.8) 
Uos 

IR Cri|Dri| | Le (428) 


Nri 
They are suitable for meaningful interpretations of circuital nature and for the sim- 


ple determination of the parameters, as it will be presented in the Sections 4.11.2- 
4.11.4. 


4.11.2 The Elements of N 7, 


By expanding (4.11), it yields 


Uo =AqiUos + Bailx » (4.41) 


where for the elements A;;,, B71, C771, Dy, precise expressions have already 
been obtained in Section 4.4. 

Now it is possible to deduce new equivalent expressions for the same elements. 
In fact, the first addendum of (4.41) ie. Az;,Up5 = U5, represents the voltage at 
H due to the application of the voltage phasor Uy, when the current x at K is 
zero (see Figure 4.21). 

In this perspective, it is important to note that the regime with J = 0 at K is not 
in contrast with the metallic continuity in K itself, since the port R is supposed to be 
also active. So by adopting the symbols shown in Figure 4.21, it results immediate 
(on the basis of the fundamental transmission equations as in Chapter 2) to write 
(being that 7x = 0) also 


Ups =ApyUSx: Udn = AUix (4.43) 
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from which it yields 
UH = (A2/Ai2)Uos => An = Ad/Arr - (4.44) 


For the determination of unknown Cj,, the first addendum of (4.42) ie. 
Cy\Uos5 = 15, can be interpreted as the current at H due to the application of 
the voltage phasor Uy, when J = 0 is imposed at K so that: 


U 
POs. and Ut = os (4.45) 
2212 
It is immediate to obtain 
‘ : C. (s 
Ty = OoU3n = Us = Cu = (4.46) 
£F12. 4712. 


By following similar considerations and passages (with reference to the scheme and 
symbols of Figure 4.22), it is possible to derive expressions in order to compute also 
the elements B;,, and Dj): 


ry) eo eo eo B, By 
Uy = -Biuls > Le =An-lLys; Von =—7 lk > Bee 
Ai Ai 
(4.47) 
ee ee ee ee A A 
Ie=Apls; Ly =A,-Ly; k= le > Dnge s 
Aj Ai 
(4.48) 


The reader can ascertain that the four elements of N ;,,; computed by (4.44)-(4.48) 
agree perfectly with those computed by (4.11) already presented in Section 4.4, so 
that the phasorial diagram of Figure 4.2 and the references in the same figure hold 
their clear meaning. 


S H K R 
s Ih I=0=1 IR 
> "> > —> 


Figure 4.21 Scheme for the determination of A;,, and C j,, 
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Figure 4.22 Scheme for the determination of B ;,, and D ;,, 
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In particular, it can be stated in general terms that 


ni = angle(C 771). (4.49) 
fy = angle(Dy7)1 x) . (4.50) 


and that the phasor U,),(5) determines, in the interval 5;-d2, current phasor J 
within the limit J,, being obviously that 


5 =a—-Bi-mtm., (4.51) 
69 = 6, + By . (4.52) 


4.11.3 The Matrix N ., 


The two elements of the matrix Ns, can be derived by applying the criteria already 
presented in the previous paragraphs: the fundamental expressions and the succes- 
sive passages do not need explanations but noting that (4.53) regards Figure 4.21 
and (4.54) Figure 4.22. 


e e e e Cc Cc 
Uos =AypUox; Ls =O Vox: Is = 37 Us => Co = 
2212 2°12 
(4.53) 
I A, 1? I? : I => D ! 
Ltr =Aptys ; fs =—7 LK Ysi=7_-: 
pes Ss Ap 1 As 
(4.54) 


4.11.4 The Matrix N p, 


In the following the elements of the matrix N p, are reported, by remembering that 
(4.55) and (4.56) refer to Figure 4.21 whereas (4.57) and (4.58) to Figure 4.22. 


Uos =AwUox: Upp =D3-Uox: 


Ubr = Us => Ap = c : (4.55) 
Uos = AU Ox ; TR = —C3-U 5x : 
LR= ~ ss => CR = = ; (4.56) 
Le = An Ii Ua = Buss; 


B 
Usp =—-— 1x => Bry =-S: (4.57) 
4412 
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TR =AplS ; Ip =Ap LS: 
A A 
a=Sle > DRG 458) 
12 12 


4.11.5 The Matrices N x, Ng» N p> 


In the following, the matrix relations which summarize in symbolic form the above 
presented equations (4.23), (4.20) and (4.29) are reported: 


Uox| |Ax2|Bx2| \Uos 

Ik C x2|D xo| | La G23) 
Nx 

Uos 1 0 | |\Zos 

Is C57|D52|| La Cy) 
N55 

Uor| | Aro| Bro| \Uos 

IR Cro Dro Ty (4.29) 


Also for the elements of these matrices, the simple expressions hereafter reported 
can be obtained (by following the same methodology of circuital analysis). How- 
ever, they have full agreement (via numerical validation) with those previously de- 
rived by matrix algorithms. 


D Cc 
Nee! aaa By, =-—; a ay ee. Dr= 7: 
(4.59) 
C 1 
Ngo: Eso = re : Dgy= rie (4.60) 
D B c A 
No? Argo = vie > Bro= aa Cro = a PR = —s 
221 231 221 pe 
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4.12 Conclusions 


The cable lines together with mixed ones can lessen sensibly many route problems 
by virtue of their adaptability in the territory and have become a recurring topic in 
the technical literature [6]. 

The authors have already presented comprehensive algorithms to deeply analyse 
the transmission performances of UGC in Chapter 3. 

In this chapter, the method has been further developed and extended to take into 
account mixed lines composed of overhead and cable lines. 

Also in these cases, the capability charts can be created, which highlights at quick 
glance the whole field of transmissible complex power compatibly with the assumed 
constraints for both current and voltage and are enhanced with other power trans- 
mission features, which result very useful to evaluate the mixed link operation. 

This gives a meaningful contribution to both the grid planning and operation. 

The cpu-time to compute a capability chart is about four seconds (PC Pentium 
2.8 GHz, RAM 1048 MB). The procedure is versatile as well since it can be easily 
applied to more simple mixed lines (d; = 0 or d3 = 0), to only cable line (d; = 
0 = d3) and to other configurations (e.g. Gas Insulated Lines). 

This chapter is also a powerful guide to choose the shunt compensation degree 
of UGC in order to take under control the transients and steady state regimes of 
energization and de-energization. 

Although the hypothesis of uniformly distributed compensation has been adopted, 
the method is easily applicable to lumped one (at each ends or also at intermediate 
locations). 
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Chapter 5 
Multiconductor Analysis of UGC 


5.1 Introduction 


The cable systems constitute paradigmatic and complex cases of multiconductor 
systems which cannot be studied in detail by means of a simplified single phase 
equivalent circuit. 

In the electrical system research, the multiconductor theory is used in several 
situations. A great attention has always been devoted to the analysis of electromag- 
netic interferences between systems of different kind, e.g. in order to evaluate pos- 
sible harmful effects due to AC power lines on telecommunication lines or metallic 
pipes: it has been well treated in a wide and detailed scientific literature [1-3]. 

Obviously, the modelling used for the above mentioned electromagnetic compat- 
ibility studies can be applied to the multiconductor analysis inside a unique power 
system [4, 5] where besides the overhead lines there are also cable lines with their 
metallic sheaths (UGC) and gas insulated lines with their metallic enclosures (GIL). 

On the other hand, the usual overhead lines at power frequency would also re- 
quire a multiconductor approach if conductors such as the earth wires and the coun- 
terpoises are realistically considered and the configuration asymmetry is not ne- 
glected. 

Moreover, in any system the soil is not a filiform conductor. 

Hence it is easy to understand how, considering the physical reality of the power 
networks, it can be questionable to assume purely three-phase configurations and 
perfectly symmetrical ones, so to use the three sequence modelling; in many cases, 
the multiconductor analysis becomes necessary, since it allows one to achieve (with 
the suitable procedures) great precision results so offering a powerful tool in order 
to validate (or less) approximated and simplified computation methods. 

It is worth premising that in this chapter examples of the multiconductor pro- 
cedures will be applied to a case already considered in Chapter 3 (i.e. cable #b, in 
cross-bonding with phase transpositions, d = 60 km) and in particular to: 
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< Major section > 


‘<——— Minor section —__> 


a b c! 


Figure 5.1 Cross-bonded single-circuit cable with phase transpositions 


* two steady state regimes (i.e. points | and 6 of Figure 3.36) of the capability 
charts (Ups = 230kV) with lumped shunt compensation every 20 km, also 
checking not to exceed the ampacity; 

¢ subtransient no-load energization with and without lumped shunt compensation. 


In the development of the procedures, the multiconductor analysis will individ- 
uate different circuital models, each characterized by their own admittance matrix: 
the suitable final composition of these partial matrices will yield the overall matrix 
which completely characterizes the system and gives the solution algorithms. 

In particular the modelling of 60 km long UGC deals with the cascade connection 
of 32 “major sections” (see Figure 5.1) where each is cross-bonded in order to avoid 
the circulation of induced currents in the sheaths, in spite of the bonding (rg) and 
earthing connections (Rg) at the ends. 


5.2 Multiconductor Cell of Three Single-Core Cables Lines 


In order to gradually develop an exposition of the criteria for the modelling and the 
analysis of the multiconductor system, let us consider (see Figure 5.2) a stretch of 
length Ag < 0.3 km between the two sections S and R of six conductors (three 
phases and three sheaths) parallel to themselves and to the ground surface where 
earth return current flows; if d >> Ag is assumed, the border effects can be ne- 
glected. 

In such a case, the treatment (given by Carson [6], Pollaczek [7] and effec- 
tively expounded and noted in [8]), shows (if the transversal couplings due to the 
phase-to-sheath and sheath-to-earth conductive-capacitive susceptances are treated 
separately) how the longitudinal ohmic-inductive self impedances z; ; and mutual 
impedances z; ; of n conductors (e.g. six as in the present case) can be computed, 
considering also the electromagnetic field inside the earth; once z; ; and z; ; have 
been computed, it is possible to form the matrix Z (6 x 6) and to characterize, by 
means of (5.1), the steady state regime of longitudinal block L of Figure 5.3 (where 
the voltage column vectors u5, Up and the current column vectors is,is;,isr, 


5.2 Multiconductor Cell of Three Single-Core Cables Lines 121 


Figure 5.3 Elementary cell (of length A,) for the three single-core cables 


Lr, iryt-Lpr are shown): 
Us —UR=ZLZyisy . (5.1) 


and by considering the evident relation (5.2) 


ir, =—ts1 > (5.2) 

it yields, (being Z not singular) 
Zpius —Zy'ur=isy. (5.3) 
-Zj'us + Z7'ur =Ir.: (5.4) 


and hence the following matrix relation (5.5) where Y; , (12 x 12) regards the block 
L circuit formed by the six longitudinal links: 


(5.5) 


ina Yra Ua 
(12x1) (12x12) (12x1) 
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In particular, it is important to mark the directions of the currents in correspondence 
to S and R (both towards the circuital block) since the study will be developed by 
means of models identified by nodal admittance matrices. 

Being Ay, sufficiently small, it is possible to lump the uniformly distributed shunt 
admittances at both ends of the cell (transverse blocks 7's; and TR) and to consider 
separately the longitudinal elements in the block L. 

The vectors of the shunt currents at the sending-end i, and at the receiving-end 
Lrr are 


, (5.6) 


ira Yra 
(12x1) (12x12) (121) 


where Y rg (6 x 6) and Y rp (6 Xx 6) are the admittance matrices of each transversal 


block T, which can be computed immediately (see Section 5.2.5) if the susceptances 
distributed on A¢/2 are lumped at each end (see Figure 5.3). 


5.2.1 The Admittance Matrix Y , to Model the Elementary Cell 


The superimposition of (5.5) and (5.6) yields 


(5.7) 


L Up 
(12x1) (12x12) (12x11) 


which completely represents the steady state regime of the elementary cell (of length 
Az). It is important to point out that in the above (5.7) the suitable matrix and vector 
partitions have been introduced: they allow highlighting and methodically grouping 
together the steady state regimes in S and in R; it will be shown (once chosen 
Ag sufficiently small, not greater than 0.3km) how this model of the elementary 
cell with its admittance matrix Y , can be repeatedly used together with the cross- 
bonding modelling to realize cascade connections and to simulate multiconductor 
systems of great length. 
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Figure 5.4a,b Possible con- oe 
figurations of 7 and j conduc- iW (+) j 7 7 
(b) 


tors < dj > 


(a) 


5.2.2 Computation of Z_, by Means of Simplified Carson—Clem 
Formulae 


If djj < 0.135Dca (see Vol. IL, page 154 in [1]), the self and mutual longitudinal 
impedances can be obtained through the Carson—Clem’s formulae (5.8) and (5.9): 


2D 
zig = rit’ -10+. f + j-4n- 1. fn a) [2/km], (5.8) 
I 
D 
Bia 10 sf + jdm 104. fin 7c ) [2/km], (5.9) 
ij 


where r; = resistance per unit length of conductor i (the skin and proximity effects 
can be computed by means of IEC 60287 [10]; it is worth remembering that the 
Milliken conductors, if enamelled, dramatically reduce them); Dca = 660,/p./f 
[m] “Carson Depth”; f = power frequency, 50-60 [Hz]; pg = electrical resistivity 
of the soil [Q m]; a = 2 GMR [m] (Geometrical Mean Radius); usually in EHV 
XLPE-cables the conductor is hollow Milliken type; for GMR, see Vol. II, page 155 
in [1]; dj; = mutual distance [m] between conductors i and j (Figure 5.4a). 

When the conductors i and j are coaxial ones, it yields x;,; = x;,; = x;,; where 
j is the external conductor and i the inner one (see Figure 5.4b). 


5.2.3 Computation of Z , by Means of Complete Carson Formulae 


If dj; > 0.135Dca, self z; ;and mutual z; ; impedances in Z, must be computed 


by means of the following complete Carson’s formulae (which have demonstrated 
its validity also at high frequency for distribution line carrier analysis [11]): 


.{ bo, 2hi 
2,4 = (i + AR) +j (of In Gu a ax;) [S2/km] , (5.10) 
2D3; 
Zz) =Zj, = ARy +i (on a de Axi) [2/km] , (5.11) 
‘ : 20 di; 


where Carson’s correction terms AR;, AX; and AR;;, AX;; are computed by 
means of series [12]; 4; is the height above ground-level of conductor i and has 
a negative value for an underground conductor. Dj;; is the distance between conduc- 
tor i and the image of conductor j or vice versa. 
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Figure 5.5 “Impedance Z1+ Zot Z3-Zy 
equivalent circuit” [9] Conductor 


Zi 


Metallic sheath e {_—i«itgk 
Z5t+ Zot Z,-Ly 
Earth 


5.2.4 Computation of Z_, After Wedepohl 


Another mean to compute Z, is the Wedepohl’s theory which gives approximated 
formulae both of Schelkunhoff loop theory [13] and of Pollaczek one [7] and which 
allows forming the matrix Z ; w differing from Z, only for a different order of the 
phases and sheaths as in (5.12): 


1-4 2.5 3 6 


(5.12) 


In order to calculate each element of the submatrices, the formulae are detailed 
in [9] so that it is immediate to obtain Z, from Zw. In particular, each diagonal 
submatrix (full grey in (5.12)) is constructed from the seven impedances Z,...Z, 
and can be derived from the “impedance equivalent circuit” of Figure 5.5. 


Phase Sheath 
GAGA 4 De FAT OF SZ. | 


Zs+Z,.+Z,-Z, Z,+Z,+Z, Sh. 


The off-diagonal submatrices in (5.12) take account of the mutual inductances be- 
tween cables and may be calculated by suitably approximating the Pollaczek’s for- 
mulae [7]. 


5.2.5 Computation of Y 7, 


The two matrices Y7; = Y 7p can be computed taking into account the shunt 
capacitive and conductive links. 
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The kilometric admittances which must be evaluated (see Figure 5.6) are ex- 
pressed by the following two formulae (see Sections 2.3.3 and 2.3.4) 


y =g,tjoc, = 10° [w- 2m - €,- (tand, +j)]/In (=) (S/km) , (5.13) 
= la 

d 
Y, = 8b +jocp = 10° [w- 2 - ey - (tan dy +\in() (S/km) , (5.14) 
~~ 1b 


in which g, and gy, represent the leakage conductances across inner and outer insu- 
lations, ¢, and ¢, [F/m] the corresponding permittivities. 


For instance, the admittance matrices Y 7, = Y 7p of the T blocks (which are 
easily deduced by the inspection method) are given by: 


1 2 3 4 5 6 
aes = 
2 ya Vi 
y “Xa 
Ay 
7 
eee Bas + yy 
Ya ya Bi Y» 

6 “Ya y te yy 

(5.15) 


When the installation is not completely embedded in soil, a precise computation of 
the self and mutual admittances regarding the sheaths would be problematic but not 
so important since the sheaths are subjected to very low voltages on average along 
the route due to the earthing and the short-circuit links of the cross-bonding: on the 
basis of this consideration it appears licit to neglect I in (5.15). 
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Figure 5.6 Scheme for the computation of kilometric shunt admittances 


5.3 Transposition Joints Modelling: Y , 


The basic idea for representing the transpositions given by the joints at minor section 
ends is to form their corresponding admittance matrices Y , having purely conduc- 
tive elements. 

Figure 5.7 shows the links which perform the transpositions of phases and 
sheaths due to the joints in a, b, c of Figure 5.1 and their corresponding admittance 
“joint matrix” Y,. 

In order to make the model more real as possible, the conductance G imputed to 
the links between S$, and Sy should assume values in accordance with the construc- 
tive characteristics of the joints. 


SS 
° _* 
oe 

phases {! _ > 4 

sheaths {3 ° 
se 
| | 


Figure 5.7 Admittance ma- 
trix Y, for cross-bonding 
joints 


5.4 Earthing of Sheaths and Insertion of Possible Shunt Reactors: Yg; Y p< 127 


Indicatively (manufactures’ data lacking) it can be assumed (for phase and sheath 
links) a conductance G of about 5kS. 


5.4 Earthing of Sheaths and Insertion of Possible Shunt 


Reactors: Y 5; Vue 
In order to complete the model of the cross-bonded circuit it is necessary to form 
the matrix Y,¢ which represents the bonding of the sheaths (between themselves 
by means of a resistance rg and with the earth by means of a cross-bonded box 
resistance Rg) at the ends of each major section (see Figure 5.1). 

There are two procedures in order to obtain the admittance matrix Y g. 

The first direct procedure is based on the inspection method which gives the 
following matrix: 


4 5 6 
re+2Re —Re —Rg 
re(re+3Re) | re(re+3ReE) | re(re+3Re) 


Y —= —RrE re+2RE —Re 
E re(rg+3Re) | re(ret+3ReE) | rete +3. Re) 
—RE —ReE reE+2ReE 
re(re+3Re) | re(e+3Re) | re(e+3.Re) 


The second procedure requires to form first the impedance matrix Z (3 x3) as seen 
at the nodes 4, 5 and 6 (corresponding to the sheaths), in which the self and mutual 
resistances can be immediately expressed as a function of rg (resistance of few mQ 
of the bonding links) and of Rg (earth resistance): 


4 5 6 
Ret+re RE Re 4 
ZE = RE RE + rE RE 5. 


Re+Trg 6 
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Figure 5.8 Admittance matrix Yor for shunt reactors and sheath earthing 


The inversion of Z gives immediately the matrix Y, (a purely conductive one) 
Ye=Z5 


which will be introduced (as it will be shown) in each major section end in the 
modelling of the whole system. 

It is worth noting that the possible presence of lumped shunt reactive compensa- 
tion (insertion between the phases 1, 2, 3 and the earth) by means of three single- 
phase reactors can be easily accounted for by means of the admittance matrix Y p¢ 
to be overlapped (as detailed in the example of Section 5.9) in the right locations. 

Analogously to Yx, there are two methods in order to compute Yy¢: a direct 
calculation by the inspection method as in Figure 5.8. 
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There is also the possibility of computing firstly the matrix Z~¢ (6 x 6), where 


Re Re Re RE Rg 
Z pe Re Re Rg RE RE 
Rp |Z~+ Re Rg RE Re 
Re RE re + Re Re Re 
RE Rg RE re + Re Rg 
RE Re Re Re Rg re + RE 


and then the admittance matrix Yp; = (Z, re) 


5.5 The Multiconductor Supply Model at the Sending-End 


It is worth mentioning that in Chapter 3 (where positive sequence models have al- 
ways been used) the steady regime had been studied by applying a voltage supply 
of Ups = 230kV directly at port S, whereas for the analysis of the no-load ener- 
gization (after Section 3.7) an equivalent subtransient impedance j X” of the supply 
network had been considered: analogously, for the study of multiconductor steady 
state regimes, the following positive sequence voltages must be directly applied 


1 
Uy = a? | -230kV, where a =e 37 , 
a 


and for the study of the multiconductor no-load energization, also a phase model of 
the supply network must be introduced by means of the matrix Z,,, (3 x 3) of the 
self and mutual longitudinal phase impedances. 
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Figure 5.9 Supply modelling (in subtransient regimes) 


Since in the subtransient regimes the equivalent three-phase supply network 
(block A of Figure 5.9) is usually hypothesized as a symmetric structure, it can 
be fully characterized by the corresponding zero, (subtransient) positive and neg- 
ative longitudinal impedances (after Fortescue [14]) whose values are coherently 
computed from the TSO in the power network studies and give rise to the sequence 
matrix Zp. 

The corresponding phase matrix Z ,,, (full matrix composed of the self and mu- 
tual longitudinal impedances of and between the phases) can be obtained (see analo- 
gous examples in [15]) by means of (5.16)—(5.18) which make use of the Fortescue 
transformation F (with Au, and i, consisting of sequence phasors) 


Au, = FAu,, : ip = Fi; Au, = Zrip ; (5.16) 
FAu,, =Zyip > Au, = | igs Ae => Au, = F'Z,F E oh : 
(5.17) 
whence 
Zo = FUL . (5.18) 
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Figure 5.10 Subtransient admittance matrix of supply model (multiconductor) 
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By following the same criteria used for the longitudinal block L of Section 5.2 


it is possible (see Figure 5.10) to achieve for the block A the matrix Y 4, which 
characterizes the power supply in the sub-transient regimes. 


5.6 Equivalent Receiving-End Matrix for Load Modelling 


In order to model a complex power loading at receiving-end, a three-phase balanced 
static load can be considered as in Figure 5.11. 


R 
: 12 3 
g|Aa/All 
Yio Yio —|4|8/4/2 
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Figure 5.11 Model of complex power demand at receiving-end 


It is worth noting that the UGC asymmetry cannot be totally zeroed, even if 
drastically reduced from cross-bonding with phase transpositions: therefore even 
with a perfectly symmetric power supply it has some voltage asymmetry at R. 

Notwithstanding, the value of y,, can be computed after (5.19) with symmetrical 


phase to phase voltages with magnitude U (e.g. U = /3Uor already computed 
with Ossanna’s method for the same hypothesized complex power S, , in the positive 
sequence modelling). 


= ; (5.19) 


5.7 The Cascade Composition of Blocks Modelled by 
‘‘Admittance Partitioned Matrices’: A First, Simple Circuit 


In order to give a clear exposition of how to compose computationally different 
network elements modelled by their partitioned admittance matrices let us consider 
the example of Figure 5.12, where the blocks a, b, c are some elementary cells of 
UGC modelled by the matrices, each with dimensions 12 x 12, Y,, Y,, Y,. given 
by (5.7) in Section 5.2.1. On the basis of (5.20)-(5.22) which lead to the general 
formulation of the current vectors injected in the different sections 1, 2, 3, 4, it can 
be ascertained that the proposed system is expressible by (5.23), where the matrix 


Y is formed of the partial overlapping (“tile style”) of the matrices Y,, Y,, Y.. 
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(5.20) 
(5.21) 
(5.22) 
1 6 24 
| peeeeteseeereeeeeeneel 
: Ys 
6: 
Esl ¥, 
in Y. (5.23) 
yl! 
24: 
i Y u 
(24x1) (24x24) (24x1) 


The classical cascade connection of blocks a, b, c becomes defined by the constrain 
i, =1,=0. 
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5.7.1 The Introduction of Other Blocks in the First Simple Circuit 
and the Steady State Analysis 


The above mentioned simple circuit can now become more realistic if the following 
additions are performed (see Figure 5.13a): 


e the bonding and earthing of the sheaths in section 1 (by means of the model 
represented by the matrix Y, as in Section 5.4) realizing a single-point bonding 
since the cable, in this first simple example with three elementary cells, is short; 

¢ the introduction in the receiving end of a balanced three-phase load modelled by 
the Y,,, as in Section 5.6; 

e the imposition (after Section 5.5) of the subvector Usp (3 x 1) at the terminals 1, 
2, 3 of section S; which causes the injection of the current subvector i py (3 x 1). 


Therefore it is possible to obtain the new relation (5.24) where: 


e the matrix Y has been updated by means of the overlapping, in the suitable posi- 
tions, of the two matrices Y,, and Y, ,; 


(Supply => - = i 
after r 5 am 4 
§ 5.5) ; ; ; + Y,, 
6 6 os 
2 Is 4,5, = 0 i= 0 i3;=0 i= 0 
Ye (§ 5.4) 


(a) 


Up 


(b) 


Figure 5.13 (a) Completion of the model in Figure 5.12. (b) Updated (a) for no-load subtransient 
energization 
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* in the structure of vector 7, the subvectors i, # 0 andi, = 0 (as highlighted 
by the Figure 5.13a). 


(5.24) 


[™ 
ad 


Relation (5.24) simulates the supply of a load (Y,,,) by means of the cascade of 
three cells. 


Y, (5.25) 


Zz 


[™ 
Is 


By introducing suitable partitions in (5.24), as shown in (5.25), (5.26) can be ob- 
tained 


by 


0 = Y3Uy, + Yu, ; (5.26) 
and thus 
u, = —-Y,'Yau (5.27) 


ph > 


where Y, is not singular and u,,, is fixed. 


So knowing the subvector u,, the steady state regime of any cell, and hence of 
the whole multiconductor system, is completely available. 
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5.7.2 The No-Load Subtransient Energization Analysis 


As already mentioned in Section 5.5, in order to study the subtransient regime of en- 
ergization at no-load, the following modifications must be performed on the circuit 
of Figure 5.13a (see Figure 5.13b): 


* zeroing the matrix Y, , in order to model a null load at R; 

* introducing the matrix Y 4,,, (regarding the supply block A) between the section 
S, and the terminals 1, 2, 3 of the section S); 

* imposing the subvector w,,, at the terminals 1, 2, 3 of the section Sy. 


Therefore, it can be derived the formulation of (5.28) where the matrix Y” (which 
includes Y 4,,,) has the dimensions 27 x 27. 


(5.28) 


u” 
haa 


By introducing the suitable partition of the matrix and of the vectors, in complete 
analogy with (5.27), it yields 


fas N-ly" 
u, = —-Y, Yu, 


and in particular the regime of the subtransient voltages at R represented by uw’. 


5.8 The Admittance Matrix Equivalent to k Blocks 
in Cascade Connections 


When a cascade connection of multiconductor blocks is considered, it may be suit- 
able (alternative or in addition to the procedure in Section 5.7) to use also a method 
which allows to synthesize a unique multiconductor block in an equivalent matrix 


Lay without increasing the matrix dimensions (see the matrix Y in (5.23)). 
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Figure 5.14 Cascade connection of two multiconductor blocks 


Let us consider first the connection of two blocks of Figure 5.14 where the admit- 
tance matrices Y,, and Y, (2 x 2n) are partitioned so that the submatrices FE... H, 
P....S, of dimensions n x n, results consistent with the dimensions of the voltage 
and current vectors. 

By observing the current vectors injected into the system, the voltage vectors in 
the different sections and the partitions of the matrices Y, and Y py, it is possible to 
write the following relations 


i,;=Eu,+Fu,, (5.29) 
i',=Gu,+Hu,, 
i, =Ru,+Su,, , 
i, =i,+in, (5.31) 
whose matricial writing is also immediate: 
iy E 0 FF uy 
i, |=; 0 S | R || w (5.32) 
dis G | @ H+Pl a, 
Since in the cascade connection i, = 0, it has 
0=G-u,+Q-u,+(H+ Pu, (5.33) 


or 


u,, = —(H_ + P)'Gu,-(H + P)'Qu, (5.34) 
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which leads to the matrix relation 


iy E-F(H+P)'G| -F(H+P)'0 uy 


(5.35) 


i, -R(H + P)'G |S—R(H+P)'O uy 


Y 


—eq 


Y ,q represents the matrix equivalent to the cascade connection of two multiconduc- 
tor blocks [16]. 

This formulation of equivalent matrix can become very useful when the cascade 
of a great number k of blocks all equal and characterized by the same cell matrix 
Y , must be considered. 

For example, if each block is represented by the same Y , (12 x 12), it is possible 
to compute the whole cascade matrix Y;, (12 x 12) (as seen at the two terminal 
ends) with a very low number of operations even if the block number k is very high. 

If for example k = 1000, the “cascade procedure” (applied nine times) yields at 
first Yoq3s Yun; Ysa -Yieas Ysa; Yoaas Yura; Yoseas Ysiza- 

The representation of k = 1000 in binary code shows that Y jo9)4 might be 
evaluated by applying the cascade algorithm to 


Ysi0a3 Yosoa Viogas Yosa: Yaa Yea; 


so that in total 14 “cascade operations” are needed [17]. 


5.9 Application of Multiconductor Analysis to the System ‘“‘Cable 
#b, 60 km”’ Already Studied in Chapter 3 with Simplified 
Criteria (see Figures 3.21 and 3.36) 


As already mentioned, this paragraph deals with one of the case studies of Chap- 
ter 3 i.e. cable #b (see Figure 5.15 and Table 5.1), with a length of 60km, shunt 
compensated with a compensation degree &, = 0.608. 

With regard of the multiconductor system, Figure 5.16 (together with Figure 5.1) 
results extremely useful in order to visualize both the cable subdivision in elemen- 
tary cells (with indication of minor and major sections) and the lumped shunt com- 
pensation locations. 
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In particular, Figure 5.16 shows that: 


¢ the minor section is suitably set equal to the cable drum length of 625 m; 

¢ each minor section is modelled as a cascade of 5 elementary cells of small length 
Ag = 125m, (after the recommendations of Section 5.2), each modelled by its 
own Y , (after the indications of Section 5.2.1); 

e at the end of each minor section, the model of cross-bonding (characterized, after 
Section 5.3, by the matrix Y ,;) must be introduced (see a, b, c of Figure 5.1); 

¢ by observing Figure 5.1, it can be seen that each major section has a length of 
3 x 625 = 1875 m (since it is modelled by the cascade of three minor sections 
included the three transposition models); at the ends of the major section the 
earthing of the sheaths must be taken into consideration by means of the matrices 
Y¢ suitably located, as depicted in Figure 5.1; 

¢ the model of 60 km long UGC (without shunt compensation) gives consequently 
the simulation of 60/1.875 = 32 major sections; and 

¢ for the lumped compensation (with &, = 0.608) it is foreseen the subdivision of 
the line in three parts constituted by 11, 10, 11, major sections. In the suitable 
locations where shunt reactors are installed, the matrix Vee (or Yi) must be 
computed: it will substitute the matrix Y, after Section 5.4. 


It is worth noting that this method also suit very well the case of major sections with 
different lengths (as it happens in the real installations). 

The model of the system must be completed with the introduction, as shown in 
Figure 5.13, of the balanced load modelling characterized by the matrix Y, , (after 
Section 5.6) by assuming the regime #1 of Figure 3.36 as reference. 

Finally the whole matrix Y (3456 x 3456) arises: after the partition of it similarly 
to the simple example! seen in Section 5.7.1, it is possible to achieve the steady state 
regime of UGC by means of the computation of the unknown vector u, due to the 
application of the supply vector w,,,, so that the multiconductor system steady state 
is completely known. 

Besides the regime #1 of Figure 3.36 many other regimes have been accurately 
investigated: they allow one to ascertain the satisfying correctness of the simplified 
method based on the positive sequence model. 

Figure 5.17 shows, as first example, the phase current magnitudes along the ca- 
ble: they are in whole accordance with those of Figure 3.37 demonstrating that the 
single-circuit positive sequence analysis of lumped compensation (for usual cable 
installations) can be sufficient for engineering purposes. Obviously, the multicon- 
ductor analysis brings to a higher accuracy of system knowledge as the zoom of 
Figure 5.18 (and of other curves) clearly highlights. 

The zoom of Figure 5.18 shows a light exceeding of the ampacity in the first 
kilometres, but it is not problematic if one considers also the conservative limits due 
to the natural fluctuations in the line operation. 

Figure 5.19 shows the phase voltage magnitudes which are almost equal to those 
of Figure 3.38: the scale does not allow separating the three different phases so that 


' The inversion of the submatrix Y, is not problematic since it is highly sparse. 
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Table 5.1 Data assumed in the multiconductor computations 


Multiconductor cell length km 0.125 
Cable drum (minor section) km 0.625 
Cross-bonding section length km 1.875 
(major section) 

Substation earthing Re Q 0.1 
Earth resistivity Pst Gm 100 
Cross-bonded box resistance Re Q 10 
Sheath resistance at 78.2 °C (50 Hz) Ts» mQ/km 70.0 
Link resistance re mQ 1 
Shunt compensation degree sh 0.608 
Number of shunt compensation stations 4 


Figure 5.15 Single-circuit 
cable line (in cross-bonding) 1 350 350mm 


R Ss T 
1 2 3 <—— phases 
4 5 6<——sheaths 


Multiconductor cell lenght=125 m 


Minor section=625 m 


Major section=625 m x 3=1875 m 


Major section=1875 m he 


Figure 5.16 Subdivision of the single-circuit cable line in cross-bonding with indication of lumped 
shunt compensation locations 


the zoom of Figure 5.20 shows the great accuracy and detail of the multiconductor 
procedure and a negligible voltage asymmetry. 

One of the great possibilities offered by the multiconductor analysis is the knowl- 
edge of the sheath electric behaviour. Figure 5.21 shows the sheath voltage and cur- 
rent magnitudes along the major sections of the cable. The zooms in the first three 
major sections in Figure 5.22 show slight differences. The typical behaviour of volt- 
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Figure 5.17 RMS phase current magnitudes along the single-circuit cable #b, (Es, = 0.608 with 
32 cross-bonding major sections) 
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Figure 5.18 Zoom of the RMS phase current magnitudes along the single-circuit cable #b, (E,, = 
0.608 with 32 cross-bonding major sections) 


age magnitudes can be also qualitatively and theoretically foreseen [18]. The sheath 
currents of induced nature are extremely low so demonstrating the effectiveness of 
cross-bonding arrangement: only the capacitive currents (since the multiconductor 
neglects the eddy currents) remain with a very complex composition. 
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Figure 5.19 RMS phase voltage magnitudes along the single-circuit cable #b, (Es, = 0.608 with 
32 cross-bonding major sections) 
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Figure 5.20 Zoom of the RMS phase voltage magnitudes along the single-circuit cable #b, (Eg, = 
0.608 with 32 cross-bonding major sections) 


As is well-known, the use of phase transposition gives a good electrical sym- 
metrization of the flat installation so that the sheath currents result further reduced 
from cross-bonding with phase transpositions to that without them. The phase trans- 
position is unnecessary in case of trefoil arrangement (which has already an intrinsic 
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Figure 5.21 RMS sheath voltage and current magnitudes along the single-circuit cable #b, (E, = 
0.608 with 32 cross-bonding major sections) 


symmetry). An important consequence is the lower power losses in cross-bonding 
rather than only sheaths cross-bonded: this implies a higher ampacity. 

The sheath voltage magnitudes along the line in the minor sections constitute 
rather high touch voltages for operator safety: they could be reduced by increasing 
the number of joints which would allow lessening the length of the major sections 
and consequently the distance between the earthing points. 

The no-load steady-state regime #6 of Figure 3.36 has also been investigated: 
Figure 5.23 shows the voltage and current magnitudes along the cable line. The 
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Figure 5.22 RMS sheath voltage and current magnitudes along the single-circuit cable #b, for the 
first three major sections (€,, = 0.608 with 32 cross-bonding major sections) 


behaviours are very similar to those shown in Figures 3.37 and 3.38. Moreover in 
Figure 5.24 the sheath voltage magnitudes are shown: it is extremely interesting 
to verify that the induced voltages in the sheaths are strictly depended upon the 
inducing phase currents. 

Also the subtransient energization with and without lumped shunt compensation 
has been investigated with the multiconductor approach; in this case for the supply 
network modelling in subtransient regime (see in Section 5.5, Figure 5.9) the fol- 
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Figure 5.23 RMS phase voltage and current magnitudes along the single-circuit cable #b, for 
no-load energization (steady state) (Es, = 0.608 with 32 cross-bonding major sections) 


lowing values of the sequence impedances have been chosen: Zz 4) =2 = j20Q; 
Zp =54+j30Q. 

Table 5.2 shows the voltage magnitudes compared with those computed in Sec- 
tion 3.13: the accordance is fully satisfactory. 

The matrix procedure here applied to the simple case of a single-circuit cable 
system can be enlarged to the study of more complex cable configurations (e.g. 
cross-bonded double-circuit in parallel electrical operation, with any configuration 
of supply phase sequence and any earthing techniques of the sheaths, [19]). 
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Figure 5.24 RMS sheath voltage magnitudes along the single-circuit cable #b, for no load ener- 
gization (€4, = 0.608 with 32 cross-bonding major sections) 
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Table 5.2 Subtransient energization voltage (kV) magnitudes Up with and without shunt reactive 
compensation compared with the results derived in Chapter 3 (for the cable #b, L = 60km, 
Ey = 0.608) 


Multi- Positive 
conductor sequence 
analysis analysis 


240.59 
240.58 240.5 
240.60 


258.92 
258.90 258.8 
258.93 


With shunt reactive compensation 


Without shunt reactive compensation 


YAR AAR 


5.9.1 Comparisons with Other Methods 


The comparison of the results obtained by means of the Carson—Clem’s formulae 
with those obtained by means of the complete Carson’s and Wedepohl’s gives a very 
good agreement: this means that, at power frequency, the use of Carson—Clem’s for- 
mulae is a powerful tool considering their elegant simplicity. Each case has been 
also compared with EMTP cable-constants [12, 20]: the percentage difference is 
a maximum of 1%, which no comparative figure of standard size can clearly high- 
light. 
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The presented multiconductor method is completely general and does not need 
any algebraic manipulation as in the previous contributions [21, 22] or loop current 
analysis [23]. 


5.10 Conclusions 


The elegance and compactness of matrix notation applied to power system (which 
has been highly emphasized by Kron’s Tensor Analysis [24]) find in multiconductor 
cell analysis a powerful and natural application. 

The possibility of choosing in the theoretic settlement the simplified Carson— 
Clem’s equations, the complete Carson’s theory or Wedepohl’s (which give the same 
results at power frequency) allows using the more suitable tool depending upon the 
cases. 

Accuracy of the method has been verified through systematic comparison be- 
tween multiconductor algorithms and the traditional software EMTP cable-constants 
(and FEM validation as well). 

The multiconductor matrix analysis has demonstrated several advantages. The 
possibility of knowing the electric quantities along the line with a detail as better 
as shorter the cell length allows determining the behaviour of all system conductors 
including the earth; it constitutes a sure reference in order to compare the results of 
more simplified methods. 

Other immediate applications of this multiconductor analysis could be: 


¢ Gas Insulated Lines considering the phases plus the enclosures; 

¢ Mixed line composed of cascade connection of overhead line (with one or two 
ground wires) and double or single-circuit UGC; and 

¢ Cables in tunnel installation considering the steel reinforcement of the gallery as 
distributed earthing or other earthing conductors [25]. 


The multiconductor analysis allows: 


¢ Detecting the exact repartition of currents in the phases and in the sheaths of 
any number of UGC circuits and with any bonding configuration (cross-bonding, 
single-point bonding, solid-bonding); 

¢ Knowing the precise steady state electric behaviour of any component of the 
cable line (i.e. cross-bonding boxes); 

¢ Considering possible phase transpositions; 

¢ Including the presence of the metallic grid of the gallery hosting the transmission 
line (if tunnel installation); 

e Studying the subtransient short-circuit regime and the voltages in any line section 
(consequently the touch-voltages); 

e Studying accurately the power loss behaviour along the line and subsequently 
the temperature distribution and ampacity along the line with a suitable thermal 
model. 
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Chapter 6 


A Comparative Procedure for AC OHL 
and UGC Overall Cost 


6.1 Introduction 


For a long period, starting from the XLPE EHV cable market injection, the high 
investment costs of EHV UGC (which can be increased by the shunt compensation), 
were often taken as a sure cause to prefer a priori an overhead line in the network 
widening or reinforcing, without devoting the right consideration to the (sensibly 
different from cable to overhead lines) burdens deriving from territory impact and 
from the energy losses in the whole line lifetime. 

The importance of these costs, at first disregarded, has become evident succes- 
sively [1] due to both the increased ecological-economical needs of the territory and 
the quality of power grid energy efficiency: consequently suitable criteria were in- 
troduced in order to evaluate the costs regarding both the territory impact and the 
energy losses. These ones (burdens of territory impact and energy losses) must be 
added (together with others of less relevance e.g. O&M, dismantling ones) to the 
mere investment cost of a new line so to have a more meaningful “overall cost” 
for a right comparative analysis. In particular, the financial burdens on the territory 
(both private and public property) caused by the line along the whole route have 
generally a not negligible economic weight: they must be computed and introduced 
(as compensation costs) in a complete evaluation of an “overall cost” due to the 
erection of a new line. 

This methodology, which takes into account the technical, environmental and 
social aspects involved in the overall economy of a new line, seems to be in whole 
accordance with the aims of the “Environmental Impact Assessment” recommended 
by the European community since 1985 [2]. 

In this chapter the following costs will be considered: 


* capital costs (I); 

¢ loss energy costs (E); 

¢ burden on territory (T); 

e ashort account on the visual impact; 

* operating and maintenance (O&M) costs (OM); 
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¢ dismantling costs (D); 
¢ UGC shunt compensation investment cost (Al)sn; and 
¢ overall cost (2). 


The costs deriving from the forced outages and the failure reparation are not 
considered as it will be explained in Section 6.7. 

The comparison procedure proposed in this chapter is limited to the power fre- 
quency technologies but in the technical literature there are similar approaches [3, 4] 
to include DC transmission. 

In the following, in order to substantiate the procedure, exemplificative values of 
the parameters are chosen: even if they may vary year by year and country by coun- 
try (especially those of economic nature), the methodology and approach remain 
always valid. 


6.2 OHL and UGC in the Comparative Procedure 


Figures 6.1 and 6.2 show the configurations of OHL (#al) and UGC (2#c1) respec- 
tively (both at 380kV, 50 Hz), which will be considered in the description of the 
comparative procedure. 

Figure 6.1 shows a typical tower used (with different suitable heights) in the EHV 
Italian grid, and reports the necessary data to calculate, by means of the usual for- 
mulae, the positive-sequence line parameters by hypothesizing phase transpositions 
and neglecting the earth wires (see Table 6.1); also in the cases of untransposed 
lines, the errors are slight for the purposes of this chapter. 

The line configuration is important also for computations of the magnetic and 
electric fields and the width of the right-of-way which results, after Italian laws [5- 
7], quite considerable (see Section 6.5) so to weight sensibly on the territory econ- 
omy and therefore in the results of the total cost comparisons. 

A traditional tower only is here considered even if overhead lines could be erected 
nowadays with more innovative design (compact towers, Foster type, etc.). 

The per unit length resistance of OHL has been computed at 75°C, conductor 
temperature when operated at the thermal limit after the Italian Standard CEI 11- 
60 [8]. It is worth underlining that the operation of an untransposed OHL at thermal 
limit current increases the voltage unbalance factor [9] and it is advisable only for 
sufficiently short lines (as those hereby considered). 

For technical evaluations at power frequency, the OHL conductance (g; or g2 
of Table 6.1) can be usually neglected, but not for economical evaluations: corona 
losses and the insulators leakage currents depend strongly [10, 11] upon the weather 
conditions which must be accurately evaluated chiefly where the rainy days often 
occur during the year. 

Figure 6.2 shows the underground installation of the proposed double-circuit 
2#c 1: its fundamental data are reported in Tables 6.2 and 6.3. 
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3 subconductors per phase 
Bundle spacing=0.4 m 
ACSR ®=31.5 mm 


Hy in=7.78 m 
hypothesis a) 


Ay in=11.34 m 
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CSS of Art.2.1.05 
K Ra Stand. CEI 11-4 


Figure 6.1 Typical configuration of an OHL in the 380 kV Italian grid 
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Figure 6.2 Installation of 380 kV cross-bonded, XLPE UGC (2 single-circuits #c1 in electrical 
parallel operation) 


Table 6.3 reports some values of the ampacity depending upon the different soil 
thermal resistivities p,, which are within a wide range foreseen in Clause 3.2.2 of 
TEC 60287-3-1 [12]. 

Owing to the great spacing (5 m), the two single-circuits #c1 are thermally and 
electrically almost independent so that, for double-circuit 2#c1, the ampacities are 
twice the J, values in Table 6.3 and the electrical model can be performed after 
Section 2.7 by means of the single-circuit #c1 parameters. 
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Table 6.1 Positive-sequence parameters of transposed OHL #al (ignoring the earth wires) 


Overhead #al 
Nominal voltage 380 kV 
3 sub-cond. 
Phase bundle ACSR @ = 31.5mm 
spacing = 0.4m 

Resistance at 75 °C (50 Hz) r mQ/km 23.10 
Inductance £ mH/km 0.858 
Capacitance c uF/km 0.0133 
Shunt leakage Fair weather £1 nS/km 10 
(50 Hz) Rainy weather 22 80 
Ampacity (Zone A-Italy) Cold Months I, A 2955 
after CEI 11-60 (Oct.—Apr.) 

Warm Months 2220 

(May-Sep.) 


Table 6.2 Cable geometrical data 


Single-core cable #c1 
Nominal voltage 380 kV 
Cross-sectional area mm 2500 Cu 
Diameter on conductor (Milliken) mm 63.4 
Diameter on XLPE insulation mm 119.9 
Diameter on metallic screen mm 130.1 
Cross-section of screen mm? ~ 500 Al 
Diameter on PE coating mm 141.7 
Mass kg/m 37, 


Table 6.3 Positive-sequence parameters of cross-bonded UGC (single-circuit) 


UGC #c1 

Nominal voltage 380 kV 
Cross-section mm? 2500 Cu 
Resistance at 90 °C (50 Hz) r mQ/km 13.3 
Inductance £ mH/km 0.576 
Shunt leakance (50 Hz) with tand = 0.0007 g nS/km 51.5 
Capacitance with ¢, = 2.3 c uF/km 0.234 
Steady state ampacity 7, with py, = 0.8 K m/W (case a) I, A 1975 

with pn = 1.0K m/W (case a) I, A 1788 
0 = 90°C (XLPE); with Py, = 1.2 K m/W (case a3) I, A 1676 
0, = 20°C with py, = 1.4K m/W (case a4) I, A 1566 

with Pm = 2.0 Km/W (case as) I, A 1327 


Pn = soil thermal resistivity 
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Figure 6.3 Capability charts 
of 2#cl (case @4) 
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sone 


00 ‘ 
-2500 -1500 -500 0 500 MW 1500 2500 


It is absolutely worth outlining that the comparison here proposed (#al versus 
2#c1) is a simple example in order to apply the procedure and does not imply a uni- 
vocal correspondence between OHL and UGC.! 

It rises not only from the great variety of the possible installations within the 
different national rules, but also from the environmental conditions, different zone 
by zone for each country, from the evolution of the different technologies and from 
the criteria more or less conservative used for planning a line. 

In order to confirm a possible operating equivalence of 2#c1 against of #al the 
corresponding capability charts must be compared so that it is necessary also for 
2#cl1 to fix the achievable level /,: it depends not only upon the electrical and ther- 
mal parameters of the cable itself but also upon the thermal parameters of the soil 
(see Table 6.3) where it is laid. It is well known that the different thermal resistivities 
of the soil (chiefly depending upon the moisture content) may strongly vary along 
the route [13] and, unless an extensive measurement campaign is performed, they 
cannot be known with great detail. It is also worth remembering the uncertainties 
of the forecast of the soil temperatures, both for the climate behaviours and for the 
possible presence of external heat sources nearby the route [14]. 

Once fixed, in a more or less conservative way, the J, level compatible with the 
environmental conditions (e.g. case a4) and the length of the line (d = 10km as 
#al), it becomes possible to compute the capability charts of 2#c1 (see Figure 6.3 


' Clause 5 of IEC 60287-3-1 gives the purchaser a useful guide in order select the appropriate type 
of cable together with the manufacturer. 
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with J, = 2+ 1566 = 3132 A) beyond that of #al (see Figure 6.4): by means of 
these charts, it can be verified that the two lines are able to face satisfyingly the 
same operating conditions (see load duration curves in Sections 6.4 and 6.10). 

Figure 6.5 clearly highlights the reduction of the receiving area for #al in “warm 
period” which must be well considered in the operation during the year: it also 
shows the possible advantages of 2#cl in more thermally favourable soils (cases 
Ql|—/3). 

Particular attention must be paid to the different foreseeable levels of availability: 
since the great repairing times of a UGC against those of OHL, the choice of a dou- 
ble circuit UGC (as here considered 2#c1), seems to be suitable for emergency op- 
eration when one of the two circuits is unavailable (e.g. see Section 6.10.2). To this 
aim, it becomes absolutely necessary to consider also the capability charts of #c1. 


6.3 The Capital Costs of OHL and UGC 


In a first approach, the investment costs of the terminal stations are approximately 
equal for OHL and UGC (single or double cable circuits) and so disregarded in the 
comparison; in a detailed analysis, they could be taken into account. 

Beyond the purchasing and installation costs, the capital cost includes the cost 
of other additional works for the lines (e.g. excavation, ducts, tower foundations, 
wayleave). 

By defining 
* a= capital cost of OHL #al per kilometre (M€/km), 

* c =capital cost of UGC 2#cl1 per kilometre (M€/km), 


the following values will be hypothesized, in order to develop the case studies of 
Sections 6.11 and 6.12, 


Cc 
a = 0.35 [M€/km]; — = 10, 


It is worth noting that the values are merely indicative since they can strongly vary 
with market conditions. 


6.4 Energy Losses and Their Actual Cost 


The economical assessment of the energy losses during the operational life of 
a transmission line plays a significant role in the overall cost evaluation. 

The energy losses which have a meaningful economic weight can be computed 
by the steady state power flow studies: the power flow diagrams (and hence the 
power losses) are strictly linked to the line typology inside the electric network 
structure: it is worth remembering for example, in EHV, cross-border interconnec- 
tions (complex power in the first and third quadrant), connections between power 
plant and grid [15-18] (complex power in the first quadrant) or lines inside the 
meshed transmission network. 
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With regard to the latter, the power flow diagrams present generally great fluctu- 
ations both on a daily and monthly basis due to the variations in the load set. For the 
other two line typologies, the power flows can present a constant behaviour on a long 
term (e.g. connections for base service power plants, cross-border interconnections) 
but also strongly variable behaviours (e.g. connections for regulation service power 
plants). 

In the present case studies (see Section 6.10), any power flow diagram (with the 
ratio between reactive and active power Or/PR = 0.2) has been represented by 
means a typical annual duration curve, which efficaciously reports the magnitudes of 
the receiving power Pr surveyed each hour of the year (see, for example, Figure 6.6 
in the hypothesis of perfect availability of 8760 h). 

In order to simplify the computation procedure, it is preferable to approximate the 
continuous curve of Pr with the step curve (see Figure 6.6), where At = 1000h 
(but the last one At = 760h), so that for each time step At it can be assumed 
a constant complex power S, = P;(1 + j0.2). Since the chosen step diagram must 
be the same for both the line typologies, a greater accuracy (e.g. with At = 500h) 
is not necessary. 

By using the Osanna’s method suitably implemented, it is immediate to compute 
the power frequency regime of the line due to each complex power S, at receiving 
end and hence the energy losses in the considered time step. 


Pm 


Figure 6.6 A typical year 
duration curve of receiving 
active power Pr (solid line) Pr 


é : i} 
and approximated step dia- 0 1000 2000 3000 4000 5000 6000 7000 8000 +9000 hours 
gram (dashed line) 
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In such a way, the different foreseeable annual energy losses AEL for OHL and 
UGC (with the same line length d) and subjected to the same duration curve of 
demand power can be evaluated. 

In the following computations, during the n years of the UGC and OHL life, 
the duration curve, the annual discount rate and the loss energy specific cost are 
assumed unchanged: by giving the energy losses a specific cost, 


¢ Loss energy specific cost = Lec, e.g. Lec = 40 [E/MWh]; 
the annual loss cost is given by 
ALC = AEL- Lec ; 


hence, the corresponding actual costs V (Vog_ and Vycc) can be computed by 


means of 
d+i)"-1 
= —— -ALC. 6.1 
idl +i)” sot) 


By assuming: 


¢ Line lifetime = 1 = 35 [years]; 
¢ Real rate of annual interest (discount rate) = i = 0.025 (= 2.5%), 


it yields: 
V = 23.15-ALC. 
In any case, all power line life long, it appears problematic to foresee: 


a) the behaviour of power flow diagrams, generally heavier as years go; 

b) the energy loss specific costs depending upon the market evolution; 

c) the different annual discount rates for the financial calculation of actual costs of 
energy losses. 


However, by implementing some algorithms of financial mathematics in year inter- 
vals with constant parameters, it is possible to develop a numerical program evalu- 
ating the sensibility of (E) and (2) to the foreseeable variations of a), b) and c) 

The values reported in Tables 6.4 and 6.5 regard idealized cases (line length of 
1 km along which the voltage and current remain practically constant) but allow the 
reader to realize the different contribution of P; and Pg, to the total power losses. 

In the OHL the losses due to the insulators and to the corona-effect increase 
severely under foul weather conditions [10, 11] which in some countries may occur 
in several days of the year. 

It is of note that the Joule losses have always a prevailing weight both for OHL 
and UGC and that OHL Joule losses are always greater than those in UGC owing 
to the major longitudinal resistance (T#a1/T24c1 & 3.5). It yields also annual energy 
losses AEL, computed on the basis of the same duration curves, much higher than 
those of UGC and subsequently Von, >> Vuac. 


? IEC 60287-3-2 gives useful comments on these economic aspects. 
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Table 6.4 Power losses [kW/km] of OHL #al and UGC 2#c1 with line length = 1km, J = 
2955 A at U, = 380kV 


I =2955A 

#al Power losses 2#cl Power losses 

kW/km En = 0 kW/km 
P, 605.10 P, 174.2 
Pei 1.45 P, 15 
Poo 11.60 
Py + Pot 606.6 Py + Py 189.2 
Py + Poo 616.7 


P; = Joule losses (UGC and OHL) 

P., = insulator and corona losses in OHL (fair weather) 
P.. = insulator and corona losses in OHL (rainy weather) 
P, = dielectric losses in UGC 


Table 6.5 Power losses [kW/km] of OHL #al and UGC 2#c1 with line length = 1km, J = 
1500 A at U;, = 380kV 


I = 1500A 
#al Power losses 2#c1 Power losses 
kW/km En = 0 kW/km 

P, 156.00 P; 44,90 
Poi 1.45 P, 15.00 
Poo 11.60 

Py + Poi 157.5 Py + Py 59.90 

Py + Poo 167.60 


P; = Joule losses (UGC and OHL) 

P., = insulator and corona losses in OHL (fair weather) 
P.. = insulator and corona losses in OHL (rainy weather) 
P, = dielectric losses in UGC 


Therefore for each kilometre of line length #al the discounted energy loss cost 
is given by: 


= (7B) e/toy, 


and analogously for UGC 2#c1: 


(E). = (==) [M€/kn] . 


6.5 The Burden on Territory 


Nowadays, notwithstanding the numerous researches on magnetic induction expo- 
sure [19-21], there is not yet a worldwide recognized magnetic limit value even 
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if the Guidelines of International Commission on Non-Ionising Radiation Protec- 
tion (ICNIRP) recommends some reference levels for the general public exposure 
to time-varying electric and magnetic fields [22]. 

In Italy for instance, there is a general discipline [5, 6] devoted to the population 
protection from electromagnetic field exposure at 50 Hz. There are different limit 
values B [uT] for short-term and long-term exposure; the short-term exposure is 
100 uT in accordance (at 50 Hz) with ICNIRP [22]. With regard to long-term expo- 
sure, for existing lines there is a value of attention of 10 uT, whereas for new lines, 
a “target of quality” of 3 uT (RMS value). 

It is worth remembering that in Switzerland for new installations the average 
magnetic field long term exposure must not exceed | uT [23]. 

The Italian Law compels to calculate the magnetic fields generated by power 
lines at the ampacity regimes and to individuate [7, 24], along the route, the width 
F of the “corridor” where B > 3 uT; in the “corridor” any extended residence (> 4h 
a day) is interdicted. 

The width F depends upon the phase currents and the line arrangement; for the 
computation of magnetic fields there are several guides [25-27]. 

Figure 6.7 shows the magnetic induction levels generated, in a plane perpendic- 
ular to the line, by the double-circuit UGC 2#cl with Iphase = 2955 A considering 
the sheath currents zeroed by the cross-bonding. In any case, the multiconductor 
analysis (see Chapter 5) would allow a more detailed analysis. It is of note that, for 
UGC, the spacing reduction and/or metallic screening decrease F, but imply a more 
dissipative installation and hence a lower ampacity. 

Figure 6.8 gives an immediate view of the magnetic induction levels generated by 
the OHL #al with phase = 2955 A and span attachments at 32-33 m, neglecting the 
current induced in the earth wires. In many cases, instead of computing the different 
corridor widths on the ground along the span, it is allowable [7, 24] to assume 
conservatively the vertical projection of F, on the ground itself. Also for OHL the 
use of mitigation screening circuits is possible but they are extremely onerous [28]. 
A very slight screening effect is indirectly (see F, = 101.3m in Figure 6.9) given 
by the earth wires, whose currents have been computed by means of multiconductor 
algorithms. 
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Figure 6.7 Magnetic induction levels for the double-circuit cross-bonded UGC (2#c1, Iphase = 
2955 A) and width F¢ of corridor 
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Figure 6.8 Magnetic induction levels for the OHL with J phase = 2955 A (span attachments at 32 
and 33 m) neglecting the presence of earth wires 
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Figure 6.9 Magnetic induction levels for the OHL with J, phase = 2955 A considering the presence 
of currents of earth wires by means of multiconductor algorithms 


The ratio F,/F, = 5.5 is a first impressive index of the different territory im- 
pacts. 

Since in the area where B > 3 uT any extended residential building is not permit- 
ted, it appears necessary to evaluate the loss of value inflicted on the territory itself 
by the new line. To this aim, it is proposed to presume for it an “average building 
parameter ed” similar to that of adjacent areas. 

The value of ed [m3/m7”] is highly variable in the territory and means the average 
ratio, in a given area, between the developed building volume and the surface of area 
itself; it ranges between 3-4 (or much greater) in urban area and 0.8—1 in suburban 
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3125 m* 


Figure 6.10 Visual impression of two areas with the same average building parameter ed = 0.1 
[m?/m?] 


area and can lessen below 0.1 in rural areas. For the sake of a visual impression, see 
Figure 6.10 where two areas have the same average building parameter ed = 0.1 
[m3 /m?], 

As a first approach, it can be considered that every square-meter of land located 
in the “right-of-way” (corridor of width F’) would lose a value: 


wy =k-ed [€/m’], (6.2) 


where k depends strongly upon the property market. A merely indicative value could 
be set equal to k = 400 €/m? so that when ed = 0-0.1 [m?/m7]. 


wy =040 [€/m’). (6.3) 


In fact, even if ed = O, a rural area inside the right-of-way has nevertheless a real 
damage since it will be no more possible to erect any farmer residence. 

In the following comparisons, the reader could ascertain that the variations of wx 
will have a great influence. 

In conclusion, depending upon the foreseeable value of wx, for each kilometre 
of line #a1, it can be supposed, on the territory extent F, - 1000 [m7], the damage: 


(TH). = Fy-10°- w, [€/km] > F,-107-w, (ME/kn] = 0.1- w, 
(with F, = 100m) ; 


analogously for UGC 2#c1: 


(1), = F.-10?- w, [€/km] > F-1073- wy, M€/km] = 0.018- w, 
(with F, = 18m). 
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6.6 The Visual Impact 


The economic evaluation of visual impact is extremely complex owing to its 
strongly subjective nature (as the value of the landscape is something very spe- 
cific and a function of local views and preferences). Notwithstanding the ambiguity, 
when a new line must be installed, this aspect could be evaluated as in Chapter 15 
of [29] and in [30, 31]. 

The economic theory offers many possible approaches but it is worth citing those 
based upon contingent evaluation [32-34]. 

Sometimes, the privilege of UGC on the visual impact is so evident that any 
economical comparison is not worth proposing. 


6.7 Operation and Maintenance (O&M) Costs 


The planned operation and maintenance of a line, during its life, implies some costs, 
which must be considered in the overall cost analysis. 

The annual cost ACogem of O&M is a percentage of investment cost. 

Of course, the annual costs ACogm must be discounted by means of 


Voem = SEE. -ACoam (6.4) 


analogously to (6.1). 

For an OHL, the maintenance costs are between 0.7 and 1% a year (with respect 
to the investment cost) and depend upon the weather and environmental conditions. 
The OHL operation cost ranges between 0.8 and 1%. The OHL O&M (flat instal- 
lation with low salt pollution) can range between 1.5 and 2%. These values (see 
Table 6.6) must be considered as an average indication and can increase in case of 
extraordinary environmental occurrences. 

With regard to UGC, once installed, they do not need particular maintenance due 
to the absence of atmospheric external situations. 

The UGC maintenance can be evaluated as 0.1% of the capital cost (see Ta- 
ble 6.6). 

The procedure does not take into account the different failure repairing times of 
OHL and UGC and their influences on the system costs. 


Table 6.6 O&M annual costs as a percentage of investment costs 


OHL UGC 
Operation 0.8-1.0% 0.1-0.3% 
Maintenance 0.7-1.0% 0.1% 


Operation and maintenance (O&M) 1.5-2.0% 0.2-0.4% 
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To this regard, [35] and [36] give some values for the failure rates of EHV OHL 
and XLPE-insulated UGC respectively. An analogous issue is developed in [37] for 
GIL technology. 

The procedure does not take into account the costs due to undelivered energy nor 
the congestion and re-dispatching costs arising from the differences in UGC and 
OHL failure repair times. Their burdens depend upon the possibilities and adapt- 
ability offered by the network meshing and consequently they are evaluable only by 
knowing the surrounding network; a more detailed analysis regarding the two line 
typologies would be misleading. 

In conclusion, for OHL #al, chosen O&M = 2%, the discounted value is given 
by: 

(OM), = 0.02-0.335- 23.15 = 0.16 [M€/km] , 


and for 2#c1, chosen O&M = 0.4% it has: 


(OM), = 0.004 - 3.5 - 23.15 = 0.32 [M€/km] . 


6.8 Dismantling or Decommissioning Cost 


A comprehensive analysis of a transmission line cannot disregard the end of life i.e. 
the dismantling phase of it. This operation foresees some costs in order to restore 
the place previously occupied by the line itself: it is evident that the considerable 
delay with respect to the investment date gives a consequent low discounted value. 

A possible order of magnitude is evaluated in [15] which gives a discounted value 
of the OHL dismantling cost: 


17.500 


= Gap 


= 0.007 [M€/km] (i = 2.5% andn = 35) . 


For UGC 2#c1, the same reference [15] reports a value of 5% of the capital cost so 
that: 


0.05 3.5 


6.9 The Cost of UGC Shunt Reactive Compensation 


The circuital modelling of UGC shunt compensation (uniformly distributed or 
lumped) has already been throughout described in Chapters 3 and 4; in this para- 
graph there are some suggestions to introduce the relevant costs (capital costs and 
energy losses). 
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Figure 6.11 Capability charts 
of UGC 2#cl by assuming an 
idealized distributed compen- 
sation 


Figure 6.12 Capability charts 
of UGC 2#cl by assuming 
lumped compensation at both 
ends 
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In order to consider a real case, let us consider the two lines #al and 2#c1 (with 
abovementioned parameters and ampacities) with the same length d = 25km. 
For such a length, the UGC 2#cl requires the compensation degree & = 0.53 
which gives the capability charts of Figure 6.11 (uniformly distributed compensa- 
tion) and of Figure 6.12 (lumped compensation at both ends considering also the 
reactor power losses). 

Figure 6.13 shows the capability charts of the uncompensated #al. 

Therefore, it is possible to verify that the capability chart of Figure 6.12 (lumped 
compensation) gives a really good agreement with that of Figure 6.11 (distributed 
compensation) and presents conservative margins (grey zones in the Figures 6.11 
and 6.12 and better levels Upr) compared with the capability chart (OHL #al) of 
Figure 6.13, so that, in abovementioned grey zones, any complex power Sp com- 
patible with #al are also compatible with 2#c1. 

The computation of annual energy losses of 2#c1 (followed by the discounted 
cost evaluation) must be performed by applying Ossanna’s method (as already sug- 
gested in Section 6.4) on the matricial model of the lumped compensated circuit 
(see Chapter 2). 

With regard to the investment costs, the increase A,, must be added to the UGC 
investment cost (see Section 6.3): it is due to the costs of 1) purchase and installation 
of the shunt reactors at both ends (usually three-phase with unchained magnetic 
flux), of ii) the disconnecting and breaking devices and the cost of iii) the territory 
occupation. 
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It is worth remembering that, in order to continue the operation also in emergency 
situation with a unique UGC circuit, it becomes necessary to foresee suitable taps 
on the three-phase reactors to reduce the compensation power, so avoiding the more 
expensive solution with two shunt reactors. 


6.10 Two Case Studies: #a1 vs. 2#c1 with d = 10km 


In the first case study, the duration curve of Figure 6.14a is assumed: it is suitable 
for a line belonging to the meshed network; in the second one the duration curve of 
Figure 6.14b, typical for a base power plant line, is considered. 

As already stated in Section 6.4, the receiving-end complex power is Sp = 
Pr(i +j - 0.2); it can be ascertained that also when Pr = 1800 MW, the corre- 
sponding S'p is compatible with the capability charts of Figures 6.3 and 6.4 (if, for 
OHL #al, the regimes at Pr = 1800 MW are hypothesized in the “cold months”). 

The energy losses and the corresponding discounted costs have been computed 
in accordance with Section 6.4; for the greater losses in #al due to foul weather (see 
Table 6.1), the hypothesis of 90 days a year with four rainy hours is assumed. 


Pr[MW] 7 Pr[MW] 


760 8760 hours 6600 8760 hours 


(a) (b) 


Figure 6.14a,b Schematic year duration curves 


6.10.1 First Case Study with Duration Curve of Figure 6.14a 


Table 6.7 reports the different cost components so that it is immediate to ascertain 
their different weights and immediately suggests to plot the overall costs (X’), and 
(*’), as a function of wx (see Figures 6.15). It is worth noting that the term (T) 
comes from the results of Section 6.5. 

Obviously, the abscissa w, = 24.3 €/m? of the intersection point (point of eco- 
nomical indifference) represents the value for which (X'), = (2), and clearly de- 
limits the intervals of convenience. 

It easy to verify that generally OHLs, having low investment costs but wide corri- 
dors F,, are convenient in lands potentially suitable for building with moderate w,. 
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Table 6.7 Composition of the overall cost [M€/km] 


#al 2#c1 
(I) 0.35 3:5 
(E) 2.05 0.70 
(OM) 0.16 0.32 
(D) 0.007 0.037 
(T) 0.1: wx 0.018 - wx 
(2) = () + (BE) + (OM) + (D) + (1) 2.57 + 0.1: wy 4.55 + 0.018 - wy 
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Figure 6.15 Overall costs 0 
per kilometre (with target of 0 10 20 24.330 40 50 60 
quality = 3 uT) wy [€/m*] 


6.10.2 Second Case Study with Duration Curve of Figure 6.14b 


The duration curve of Figure 6.14b can match the annual diagram of the complex 
power Sr = PrRj0.2) which must be delivered at node R by a thermoelectric power 
plant (three groups rated 320 MW) located at node S (d = 10km); the reduced 
value 640 MW for 2160 hours corresponds to the unavailability of each group for 
the planned maintenance of 30 days. 

Table 6.8 shows the composition of overall costs where the weights of the en- 
ergy losses are lower than those of Section 6.10.1 so that the point of economical 
indifference (see Figure 6.16) is wy = 33.3 €/m”. 

It is worth noting that the maximum power Pr is compatible with the capability 
chart of a single circuit #c1 (with ampacity of case a4) highlighting very important 
features for generation delivery continuity. 


Table 6.8 Composition of overall cost [M€/km] 


#al 2#c1 
(I) 0.35 3.5 
(E) 1.04 0.41 
(OM) 0.16 0.32 
(D) 0.007 0.037 
(T) 0.1: wy 0.018 - wy 


(Z) = () + (E) + (OM) + (D) + (1) 1.56 + 0.1- wx 4.27 + 0.018 + wx 
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Figure 6.16 Overall costs (2) [ME/km] 
per kilometre (with target of 
quality 3 uT) 


6.10.3 Sensitivity to the Principal Parameters 


Figure 6.17 shows (2’) as a function of different exposure limits set by countries for 
a location with wx = 24.3 €/m? which could occur in a rural area. 

Once again, the UGC provides a rather favourable solution when restrictive limits 
of B are compelled (in Figure 6.17 for 0.2uT < B < 3.0uT). 

Figure 6.18 highlights the sensitivity of (’) to the ratio of the investment costs 
(l)-/(l)a ranging between 8 and 10 and allows verifying how the variations of UGC 
investment cost (due to possible trader allowances, scale effect) play a key role in 
the comparative evaluations. 
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6.11 Case Study of Section 6.9 with Duration Curve 
of Figure 6.14a 


The compensation power Omyar, necessary on the whole for the double-circuit line 
2#c1 (with &, = 0.53), can be estimated about 


Omvar = 9.53» (11+ 2)pivar/km] * 25 fm] = 291.5 [Mvar] . 


Detailed evaluations of the whole system at both ends (where the cost of i) is the 
prevailing one) have brought to As, & 4.96 [M€] corresponding to a kilometric 
investment cost (Al)sn & Ag = 0.198 [M€/km]. 

Table 6.9 shows the composition of overall costs where the weight of (Al)sn 
is clearly detectable; Figure 6.19 highlights the point of economical indifference 
Wx = 27.4€/m’ slightly greater than that of Figure 6.15. 

It is worth remembering that for d > 25 km intermediate shunt compensations 
and greater &, are needed as explained in Chapter 3. 


Table 6.9 Composition of overall cost [M€/km] 


#al 2#c1 
(I) 0.35 3.5 
(E) 2.096 0.801 
(OM) 0.16 0.32 
(D) 0.007 0.037 
(AD sh 0 0.198 
(7) 0.1- wx 0.018 - wx 


(5) = () + (CE) + (OM) + (D) + (ADsn + (1) 2.61 + 0.1- wy 4.86 + 0.018- wy 


(2) [MEé/km] 
0 
FF, = 100m 


Figure 6.19 Overall costs per 
kilometre with compensated 0 
UGC (d = 25km) and year 0 10 20 27.430 40 50 60 
duration curve of Figure 6.14a wr [€/m*] 
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6.12 Conclusion 


In this chapter, the importance of an overall cost evaluation between UGC and OHL 
has been clearly highlighted. 

Ensuring reliable and economic connections and respecting the environment is 
a crucial task, often requiring innovative solutions. 

Whereas overhead lines have been the selected solution, extra high voltage ca- 
bles have made a breakthrough with the introduction of XLPE insulation reducing 
dramatically the losses whilst maintaining an excellent performance level. 

In spite of the low investment cost of overhead lines, UGCs have other important 
tangible benefits as well as some advantages which are less tangible. 

This chapter makes a detailed technical and economic assessment of these two 
different technological solutions, considering the best social and economic solution 
in different conditions and shows that the evaluation must be carried out on a case 
by case basis. 

The general procedure and approach have been reported by presenting a particu- 
lar Italian situation as an example. The model however can be applied to any country 
with its own specific transmission standards, rules and/or laws. 

In general, overhead lines and cables have been debated as competitors in a com- 
parison without precise rules. It is apparent that where land has already been devel- 
oped for residential use or where development potential is very high, underground 
cables are the preferred option having less environmental impact even if there are 
higher capital costs. Conversely, in areas of low value overhead line transmission 
having a greater environmental impact and lower capital costs remains the preferred 
option. 

The authors are aware that a “comprehensive verdict” in the technical, econom- 
ical, social and environmental fields is highly difficult and complex; nevertheless, 
they believe that the present schematic procedure can be deepened and enlarged to 
other aspects chiefly if it is performed during the planning phase when the differ- 
ent constructive, environmental, social and financial parameters can be accurately 
assessed. 

More and more in the future the overall cost comparison between OHL and UGC 
will have to be evaluated: the authors hope to have contributed in constructing the 
engineering basis of this comparison. 
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